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INTRODUCTION 


The  Naval  Research  Laboratory  (NRL)  Is  conducting  a  Special 
Focus  Program  in  Advanced  Accelerators,  l.e..  Modified  Betatron  Design. 
This  program  has  been  stimulated  by  the  rapid  advance  In  Free  Electron 
Laser  (FEL)  research  over  the  past  several  years.  Great  progress  has 
been  made  In  both  the  theoretical  understand! ng  of  the  mechanism  and  the 
experimental  demonstration  of  the  validity  of  many  features  predicted  by 
theory.  The  FEL  has  attracted  this  Interest  because  of  Its  promise  as  a 
tunable,  powerful,  efficient  source  of  radiation  over  a  very  large  fre¬ 
quency  range.  The  strong  dependence  of  output  frequency  on  electron  en¬ 
ergy  (  oj « y2)  makes  scaling  attractive.  Furthermore,  the  availability  of 
very  high  power  electron  beams  (>>  GW)  makes  possible  the  generation  of 
high  power  radiation  In  regions  of  the  spectrum  when  other  practical 
high  power  sources  do  not  exist. 

JAYCOR  has  supported  this  program  through  Investigations  In 
the  design  and  evaluation  of  several  subsystems  Including  the  coll  de¬ 
sign,  the  application  of  the  long  pulse  induction  llnac,  and  most  Impor¬ 
tant,  the  Injection  scheme  for  the  accelerator.  The  results  of  these 
Investigations  are  presented  In  the  following  publications  produced  dur¬ 
ing  this  effort. 

The  author  wants  to  acknowledge  the  support  of  the  NRL  staff 
and  co-contractors.  Special  appreciation  Is  given  to  Dr.  J.  A.  Pasour, 
Dr.  J.  Golden  and  Dr.  C.  W.  Roberson  of  the  NRL  staff  for  their  guidance 
and  encouragement. 
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A  FREE  ELECTRON  LASER  DRIVEN  BY  A 
LONG  PULSE  INDUCTION  LINAC 

X.  Introduction 

In  the  last  five  years,  free  electron  laser  [FED  research  has  expanded 

rapidly.  Great  progress  has  been  made  both  in  the  theoretical  understanding 

of  the  mechanism  and  in  the  experimental  demonstration  of  the  validity  of  many 

features  predicted  by  theory.  The  FEL  has  attracted  this  interest  because  of 

its  promise  as  a  tunable,  powerful,  and  efficient  source  of  radiation  over  a 

very  large  frequency  range.  The  strong  dependence  of  output  frequency  on 

2 

electron  energy  (u  *  y  )  makes  scaling  attractive.  Furthermore,  the 
availability  of  very  high  power  electron  beams  (>>1  GW)  makes  possible  the 
generation  of  high  power  radiation  in  regions  of  the  spectrum  where  other 
practical  high  power  sources  do  not  exist. 

In  this  paper  we  will  present  an  overview  of  FEL  research  and  then 
describe  in  some  detail  an  ongoing  experiment  designed  to  produce  mm-wave 
output  from  a  long  pulse  induction  linac  (t  •  2  usee).  In  Section  II  we 
will  describe  the  FEL  mechanism,  and  a  brief  history  of  FEL  experimental 
research  will  be  presented  in  Section  III.  Section  IV  will  describe  the 
various  FEL  operating  regimes  and  discuss  the  beam  requirements  for  operation 
in  the  Raman  regime.  Accelerator  development  and  diode  characteristics  will 
be  discussed  in  Section  V  and  the  long  pulse  induction  linac  described.  Our 
FEL  experiment  will  be  described  in  Section  VI  along  with  a  comparison  of 
various  wigglers.  The  experimental  results  will  be  presented  in  Section  VII 
an/  future  '  .actions  for  the  experiment  and  conclusions  in  Section  VIII. 

M ana*'  A  approved  December  23,  1982. 
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II.  PEL  Mechanism 


Tto  understand  the  PEL,  let  us  consider  a  very  low  current  electron  beam  so 
that  space  charge  forces  are  not  important.  We  assume  the  initial  velocity  of 
the  electron  beam  is  in  the  axial  direction  and  we  have  a  simple  linear 
wiggler  in  the  y  direction. 


Let 


and 


Then 


3  -  B  sin(k  2)  y  • 

w  w  w 


v  »  v  cos  (k  z)  x  . 
w  w  w 


(1) 


(2) 


(3) 


Here,  vfa  is  the  beam  velocity,  Bw  the  wiggler  magnetic  field,  and  vw  is  the 
wiggle  velocity  that  results  from  the  Lorentz  force  acting  on  the  particle  as 
it  passes  through  the  wiggler  field.  We  now  assume  the  presence  of  a  linearly 
polarized  radiation  field 

A  A 

E  ♦  B  »  E  cos  (kz  -  ut)  x  +  B  sin  (kz  -  wt)  y  ... 


This  radiation  field  will  exist  as  part  of  the  noise  spectrum  in  the  case  of 
the  oscillator,  or  will  be  supplied  externally  in  the  case  of  an  amplifier. 


2 


As  a  result  of  the  wiggler  and  the  radiation  field,  a  pondermotive  force 
in  the  z  direction  develops 


F  =3(v  x  B  +  v  x  3  )  z  «  sin:  (k  -  :<  }  z  -  mt]  z 

—  c  w  —  —  w  —  W  (3) 


This  force  arises  from  the  interaction  between  v^  and  the  magnetic  component 
of  the  radiation  field  and  from  the  perturbed  velocity  v_  (which  results 
from  v^  x  3)  interacting  with  the  wiggler  field. 

This  pondermotive  force  drives  a  current,  From  the  continuity 

equation  q  3Sn/?t  a  n  ,  we  obtain  a  density  modulation  of  the  form 

<$n  «  cos  [  (k  +  k  )  z  -  at]  .  (6) 

—  w 

This  density  modulation  of  the  electron  beam,  driven  by  the  pondermotive 
force,  results  in  a  current  which  can  cause  radiation  in  phase  with  the 
radiation  field: 


ij  *  q  in  v  «  cos  (kz  -  mt)  x  .  (7) 

—  w 

Therefore  the  pondermotive  force  causes  the  beam  to  bunch  and  the 
electrons  to  radiate  coherently  and  in  phase  with  the  existing  radiation. 
Thus  the  original  radiation  grows,  which  increases  5n,  increasing  5J  and 
so  on  until  saturation  occurs. 
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If  <5n[(k  +  k  ) ,  <u]  is  at  a  frequency  and  wave  number  that  is  also  a  beam 
w 

mode,  then  we  have  a  collective  free  electron  laser.  Then  to  determine  the 
wavelength  scaling,  we  assume  the  phase  velocity  of  the  pondermotive  wave  is 
near  the  beam  velocity 


v  m  -  *•  v 

ph  k+k  b  * 

w 


(8) 


Then  using  a)  =  ck ,  we  obtain 


,  1-8  .  w 

\  M  — —  X  *  - - 

3  w  „  2 

2y 


(9) 


This  is  the  wavelength  scaling  law  that  is  such  an  important  feature  of  the 
PEL.  If  the  beam  current  is  high,  the  pondermotive  force  will  drive  a 
collective  space  charge  wave.  This  enhances  the  interaction  and  increases  the 
intrinsic  efficiency.  However  the  higher  current  introduces  an  additional 
complication,  that  of  the  beam  rapidly  expanding  due  to  its  own  space 
charge.  The  simplest  configuration  to  confine  the  beam  is  an  axial  magnetic 
field.  This  introduces  the  cyclotron  modes  on  the  beam,  which  can  result  in 
competing  interactions. 


III.  FEL  Experimental  Research 


As  is  typical  in  a  rapidly  developing,  complex  field  such  as  free  electron 
lasers,  experiments  have  lagged  behind  theory  to  a  large  degree.  In  the  last 
2-3  years,  however,  scores  of  experiments  have  been  proposed  and  undertaken  to 
explore  the  many  aspects  of  FEL  operation  predicted  by  theory.  Although  only 
preliminary  results  have  been  reported  from  most  of  these  experiments, 
indications  are  that  the  FEL  mechanism  is  a  viable,  tunable  source  of 
radiation  from  millimeter  waves  to  the  infrared.  Experiments  now  underway  are 
attempting  to  extend  the  range  of  operation  into  the  visible  and  even  beyond. 

Free  electron  laser  experiments  have  evolved  along  two  distinct  paths 
depending  on  the  type  of  electron  accelerator  used.  High  current,  relatively 
low  particle  energy  electron  beams  typically  from  Marx  type  generators  have 
been  used  for  FEL  experiments  in  the  collective  regime.  More  conventional 
accelerators  such  as  RF  linacs,  which  produce  much  higher  energy  electrons  but 
at  much  lower  current  levels,  have  proved  useful  for  FEL  research  in  the 
Compton  or  single  particle  regime.  Madey  and  coworkers1  at  Stanford  have 
pioneered  the  Con?>ton  regime  research  using  a  superconducting  linac  together 
with  a  helical  magnetic  wiggler  field  (Elias  et  al. ,  1976;  Deacon  et  al. , 
1977).  They  initially  performed  amplifier  experiments  at  X  ■  10.6  pm  and 
later  demonstrated  laser  oscillator  operation  at  X  *  3.4  um  using  a  43  MeV 
beam  having  a  peak  micropulse  current  of  2.6  A.  T3ie  efficiency  (laser 
energy/beam  energy)  of  the  oscillator  was  less  than  10  ,  but  it  was  suggested 

that  a  much  higher  efficiency  could  be  achieved  by  providing  for  multiple 
passes  of  the  electron  beam  through  the  FEL  interaction  region.  However,  the 
increasing  beam  energy  spread  with  each  pass  becomes  a  limiting  factor. 
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Mora  recently ,  the  Stanford  group  has  collaborated  with  LURE  on  an 
experiment  using  the  ACO  storage  ring  in  Orsay,  France  (Bazin  et  al. ,  1982). 

In  this  experiment  a  superconducting  undulator  (B^  »  4  kG,  L  ~  1  m)  was 
placed  on  a  straight  section  of  the  storage  ring,  which  was  typically  operated 
at  150  MeV  with  a  bunch  current  of  ~  10  raA.  An  argon 

laser  (X  *  4880  A  or  5145  A)  having  a  power  density  of  -  1.6  kW/cm  was 
amplified  in  the  interaction  region.  The  maximum  gain  reported  was  ~  4  x  10~4 
per  pass,  hence  oscillator  experiments  would  be  extremely  difficult  with  the 
original  apparatus.  Consequently,  the  group  is  attempting  to  increase  the 
gain  by  substituting  an  optical  klystron  for  the  undulator. 

Somewhat  similar  experiments  are  underway  at  Frascati  using  the  Adone 
storage  ring  (Barbini  and  Vignola,  1982)  and  at  Srookhaven  using  the  VUV 
storage  ring  (Luccio,  1982)  .  Another  storage  ring  FEL  experiment  has  been 
proposed  for  BESSY  in  Berlin  (Gaupp,  1982). 

Another  set  of  similar  experiments  employing  rf  linacs  in  conjunction  with 
C02  lasers  and  planar,  permanent  magnet  wigglers  have  recently  been 
undertaken.  These  experiments  were  initiated  at  Los  Alamos,  TRW,  and  Math 
Sciences  Northwest.  The  permanent  magnet  wigglers  used  in  these  experiments 
are  convenient  for  studying  variable  wiggler  efficiency  enhancement  schemes  by 
changing  either  the  wiggler  period  or  amplitude  with  axial  distance.  In 
initial  experiments,  the  groups  have  concentrated  on  measuring  the  energy  loss 
of  the  electrons  passing  through  the  interaction  region,  because  the  gains  or 
efficiencies  are  so  low  that  it  is  very  difficult  to  accurately  measure  the 
amplified  output  signal  in  the  presence  of  the  large  input  laser  signal  (20  - 
1000  MW) , 


6 


Based  on  measurements  of  electron  energy  loss,  the  TRW  group  has 
calculated  an  efficiency  of  0.07%  and  a  gain  of  2.7%  using  a  wiggler  with  an 
amplitude  taper  of  2.25%  (Boehmer  at.  al. ,  1982).  This  efficiency  is  an  order 
of  magnitude  higher  than  the  theoretical  value  with  no  taper.  By  using  a  more 
powerful  laser  than  the  TRW  group  ( 1  GW  vs  2  MW)  the  Los  Alamos  group  can 
achieve  high  efficiency  (~  2%)  but  lower  gain  (Warren  et.  al. ,  1982).  (Mote 
that  the  input  laser  power  is  not  included  in  the  efficency  calculation) .  The 
Los  Alamos  experiment  is  now  being  converted  to  an  oscillator  mode. 

Additional  low-current-regime  experiments  are  being  performed  using 
microtron  accelerators  at  Bell  Labs  (Shaw  and  Patel,  1982)  and  at  Frascati 
(Bizzarri  et  al. ,  1982).  These  accelerators  operate  at  up  to  -20  MeV  with 
peak  current  of  -  5  A,  and  the  Frascati  microtron  is  being  upgraded  to  >  30 
MeV.  Bell  Labs  is  using  a  10  m  long  helical  wiggler  with  a  20  cm  period  to 
generate  100  -  400  urn  radiation  by  tuning  the  beam  energy  from  10-20  MeV. 
Frascati  uses  a  2.25  m  long  permanent  magnet  wiggler  with  a  period  of  5  cm  for 
output  at  -  15-20  urn. 

A  different,  but  quits  interesting  approach  is  being  pursued  at  the 
University  of  California  -  Santa  Barbara  (Elias  and  Ramian,  1982).  The  UCSB 
group  is  using  a  recirculating  electrostatic  accelerator  to  achieve  high 
average  power,  good  beam  quality,  and  high  overall  FEL  efficiency.  If  all  the 
beam  could  be  recovered,  the  device  could  essentially  operate  dc.  The  group 
hopes  to  use  the  2A,  3  MeV  electron  beam  to  generate  FIR  to  submillimeter 
radiation  in  a  single  stage  FEL,  and  then  to  perform  a  two  stage  experiment  to 
generate  visible  to  IR  output. 


High  current,  relatively  low  partiele  energy  FEL  experiments  have  achieved 
impressively  high  gains  and  radiation  power  levels  at  much  lower  frequencies 
than  the  high  particle  energy  experiments.  However,  these  experiments  are 
characterized  by  a  different  set  of  problems.  High  current  beams  are 
typically  produced  by  Marx- type  generators  and  in  general  have  higher  energy 
spreads  than  low  current  beams.  Self  field  effects  become  important  in  the 
propagation  of  these  beams,  and  the  axial  magnetic  field  requried  to  confine 
the  beam  can  result  in  large  amounts  of  cyclotron  type  emission. 

The  first  operation  of  a  device  employing  what  we  now  call  the  FEL 
mechanism  wee  reported  over  20  years  ago  by  Phillips  (1960).  He  called  the 
device  a  Ubitron  and  built  several  microwave  tube-type  versions  which  behaved 
impressively  in  S-band  (~  2-3  GHz).  The  wiggler  typically  consisted  of 

alternately  wound  co-axial  coils  and  the  interaction  length  was  -  50  cm.  Peak 

microwave  output  power  in  excess  of  1  MW  was  achieved  using  a  beam  of 

-  150  kV  and  <  100  A,  corresponding  to  an  efficiency  of  -  10%.  Scaling  this 

device  into  the  mm  wavelength  range  was  restricted  at  the  time  because  of  the 
limited  voltage  of  available  electron  guns.  However,  the  introduction  of  high 
current,  higher  voltage  Marx  type  generators  made  possible  a  renewed  effort  in 
the  generation  of  high  power  mm  waves  using  this  mechanism. 

High  current  (>  1  kA)  FEL  experiments  were  originally  referred  to  as 
stimulated  scattering,  and  their  origins  can  be  traced  to  early  experiments  at 
Cornell  Ohiversity  and  the  Naval  Research  Lab  in  which  large  radiated  powers 
were  observed  when  an  intense  beam  was  modulated  by  a  slow  wave  structure  or  a 
rippled  magnetic  field  (Nation,  1970;  Friedman  and  Herndon,  1972).  These  and 


similar  experiments  were  interpreted  using  one  of  two  different  theories,  the 

cyclotron  maser  instability  (Sprangle  and  Manheimer,  1975)  and  stimulated 

Raman  scatter ing( Sprangle  et.  al,  1975),  which  predicted  large  growth  rates 

and  high  efficiencies.  The  cyclotron  maser  instability  is  a  purely 

relativistic  effect  associated  with  relativistically  gyrating  electrons.  This 

mechanism  has  been  successfully  used  in  the  development  of  gyrotrons, 

relatively  compact  tube-type  devices  which  operate  at  the  fundamental  or 

harmonics  of  the  cyclotron  frequency  (See,  e.g.,  Granatstein,  ed. ,  1981). 

Stimulated  scattering  experiments,  on  the  other  hand,  have  continued  to  employ 

Marx  generators  at  higher  voltage  and  current  levels  in  order  to  take 

2 

advantage  of  the  y  frequency  scaling  and  the  growth  rate  scaling  with 
current.  After  the  Stanford  group  popularized  the  term  "free  electron  laser" 
and  the  equivalence  of  the  FEL  and  stimulated  scattering  mechanisms  was  shown 
(Kroll  and  McMullin,  1978;  Sprangle  et.  al. ,  1979),  the  high  current 
experiments  also  began  to  be  referred  to  as  FEL's. 

The  first  high  current,  relativistic  beam  experiment  to  be  interpreted  as 
stimulated  scattering  was  performed  at  the  Naval  Research  Laboratory 
(Granatstein  et  al. ,  1974).  In  experiments  designed  to  produce  microwaves  at 
\  ■  2  cm  via  the  cyclotron  maser  interaction,  strong  submillimeter  radiation 
was  also  observed.  A  quantitatively  consistent  explanation  was  that  some  of 
the  2  cm  radiation  reflected  off  the  output  window  and  subsequently  interacted 
with  a  sufficiently  cold  part  of  the  1.5  MeV  beam  to  produce  the  high 
frequency  scattered  radiation.  In  a  subsequent  experiment  designed  to 
optimize  this  effect,  400  um  radiation  was  generated  at  a  power  of  —  1  MW 
(Granatstein  et  al. ,  1977). 
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The  next  important  step  in  the  development  of  high  current  FEt,' s  came  in 
an  experiment  at  Columbia  university  in  which  the  high  power  electromagnetic 
pump  was  replaced  by  a  periodic  magnetic  wiggler,  or  quasi-wave  (Mross  at  al. , 
1976;  Efthimion  and  Schlesinger,  1977).  A  750  keV,  5  kA  beam  was  used  to 
generate  >  10-100  kW  of  mm  or  cm  radiation.  Subsequent  experiments  at 
Columbia  increased  the  millimeter  wave  power  level  to  more  than  1  MW  (Marshall 
et  al.,  1977a).  %is  group  also  performed  mode  structure,  spectral,  and 
growth  rate  measurements  which  indicated  that  the  mechanism  responsible  for 
the  observed  output  was  weak  pump  Raman  scattering  (Marshall  et  al.,  1977b; 
Gilgenbach  et  al.  ,  19  79)  . 

These  early  experiments  at  MRL  and  Columbia  led  to  a  collaborative  effort 
on  SRL's  VESA  generator  in  which  a  quasioptical  cavity  was  used  to  provide 
feedback  (McDermott  et  al. ,  1978).  Approximately  1  MW  of  400  urn  radiation  was 
coupled  out  of  the  cavity,  and  the  output  wavelength  agreed  with  theoretical 
predictions  for  the  y  *  3.4  beam.  An  important  result  of  this  experiment  was 
the  demonstration  of  line  narrowing  when  feedback  was  introduced.  The  line 
width  decreased  from  Al/X  >  10%  to  AX/\  •  2%. 

me  efficiency  (<  0.03%)  of  these  high  current  experiments  was  strongly 
limited  by  the  poor  beam  quality.  Also,  researchers  at  MIT  have  shown  that 
the  output  from  such  experiments  can  easily  be  dominated  by  cyclotron  emission 
(Shefer  and  Bekefi,  1982).  Consequently,  the  NRI,  VEBA  group  (Parker  et  al.  , 
1982)  designed  an  improved  diode  in  which  90%  of  the  diode  current  was  removed 
from  the  beam  with  a  collimator.  The  collimated  1.5  kA,  1.4  MeV  beam  had  an 
instantaneous  axial  velocity  spread  of  -  0.1%.  The  beam  was  passed 
through  a  helical  wiggler  and  generated  -  35  MW  of  radiation  at 
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X  »  4  mm.  This  power  was  ~  2.5%  of  the  propagated  beam  power.  A  strong 
frequency  dependence  of  the  output  on  the  axial  guide  field  was  observed,  with 
the  highest  output  occuring  at  a  field  slightly  above  the  cyclotron 
resonance. 

Another  recent  collective  regime  experiment  at  Columbia  University  has 

been  reported  in  which  the  electrons  were  given  an  initial  transverse  velocity 

before  their  injection  into  the  wiggler  (Grossman  et.  al. ,  1983) .  Output  was 

observed  at  a  frequency  equal  to  the  sum  of  the  doppler  shifted  cyclotron  and 

2 

usual  FEL  frequencies;  i.e.,  ■  2^  (fl^/y  k^v^)  .  In  this  experiment, 

<  1  MW  of  power  at  X  ■  1.5  mm  was  generated.  This  mechanism  was 
independently  proposed  by  the  MIT  group,  which  called  it  a  Lowbitron  (McMullin 
and  Bekefi,  1981),  and  experiments  to  study  the  effect  are  also  underway 
there. 

Finally,  an  important  trend  that  is  occurring  in  FEL  research  is  an 
attempt  to  bridge  the  gap  between  the  two  current  regimes.  High  current 
experiments  are  being  scaled  to  higher  voltages  in  an  attesq>t  to  produce 
higher  output  frequencies  and/or  higher  powers.  In  an  experiment  being 
undertaken  on  the  ETA  induction  linac  (V  ~  4.5  MV)  at  Livermore,  a  high 
current  beam  (  ~  1  kA)  is  being  used  in  an  attespt  to  efficiently  generate 
very  high  FEL  output  powers  at  f  ~  100  GHz  (Prosnitz  and  Sessler,  1982) .  At 
the  Naval  Research  Laboratory,  a  program  is  underway  to  develop  contact  high 
current,  and  high  voltage  accelerators,  which  could  be  used  as  FEL  drivers 
(Roberson  et  al. ,  1982).  The  long  pulse  duration  experiment  reported  here  is 
part  of  that  program  and  is  sealable  to  high  energies. 
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IV.  FEL  Operating  Regimes 


There  are  three  regimes  of  free  electron  laser  operation.  The  operating 
regimes  can  be  characterized,  in  part,  by  the  velocity  distribution  relative 
to  the  phase  velocity  of  the  pondermotive  wave.  A  phase  space  plot  is  shown 
in  Figure  1. 

In  the  Raman  regime  the  beam  appears  cold  to  the  pondermotivs  wave.  There 
are  no  particles  in  resonance  with  the  wave.  The  amplitude  of  the  wave  is  too 
small  to  trap  any  beam  particles  at  z«0.  As  the  wave  grows  in  space,  the 
amplitude  becomes  large  enough  to  trap  the  beam  particles,  thus  terminating 
the  linear  growth  phase. 

In  the  kinetic  Compton  regime,  the  beam  appears  warm  to  the  pondermotive 
wave.  There  are  particles  in  resonance  with  the  pondermotive  wave.  The 
linear  growth  is  proportional  to  the  slope  of  the  distribution  function.  This 
is  a  Landau  type  growth  mechanism.  Kinetic  theory  is  required  to  describe  the 
FEL  process  in  this  regime.  The  pondermotive  wave  is  resonant  with  a  range  of 
particles,  as  indicated  by  the  dashed  lines  in  the  figure.  The  growth 
mechanism  of  the  FEL  in  this  regime  has  a  quantum  mechanical  analogy  with 
Conqoton  scattering.  It  is  not  a  single  particle  effect,  since  it  requires  a 
kinetic  description,  and  so  it  is  referred  to  as  kinetic  Compton.  The  phase 
space  plot  shows  the  growth  of  the  wave  in  space.  Initially  the  wave 
aiqplitude  is  small  conpared  to  the  width  of  the  resonance.  The  pondermotive 
wave  will  grow  until  the  amplitude  becomes  large  enough  to  trap  the  resonant 
particles,  thus  terminating  the  linear  growth  phase. 
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Figure  1.  Velocity  distributions  and  phase  space  plots  of  the  three  FEL 


operating  regimes . 


In  the  Compton  regime,  the  beam  energy  is  high  and  the  current  low. 

Single  particle  effects  dominate  over  collective  effects.  'Bie  phase  velocity 
of  the  pondermotive  wave  is  less  than  the  beam  velocity.  There  are  two  modes 
in  which  a  single  particle  Compton  FEL  may  be  operated.  If  the  amplitude  of 
the  pondermotive  wave  is  too  small  initially  to  trap  the  beam  particles,  as  in 
the  Stanford  experiment,  the  gain  is  due  to  the  perturbation  of  beam  particles 
by  the  pondermotive  wave.  The  interchange  between  wave  and  particle  energy  is 
oscillatory  at  the  same  period  as  the  pondermotive  wave.  These  particles  are 
represented  by  a  shaded  region  above  the  pondermotive  wave  ellipses  in  Fig.  1. 

If  the  wave  amplitude  is  sufficiently  large,  as  in  the  Los  Alamos 
experiment  (Warren  et.  al. ,  1982) ,  some  or  all  of  the  beam  particles  may  be 
trapped  in  the  potential  ’veil  of  the  pondermotive  wave.  The  gain  of  the  FEL 
in  this  region  is  due  to  the  loss  of  particle  energy  as  the  particles  rotate 
in  the  potential  troughs  of  the  pondermotive  wave.  This  occurs  at  the  trapped 
particle  bounce  frequency.  The  maximum  gain  is  obtained  when  he  beam 
particles  are  at  the  bottom  of  the  ellipse  in  phase  3pace. 

The  efficiency  of  the  free  electron  laser  can  be  increased  in  all  of  these 
schemes  by  decreasing  the  phase  velocity  of  the  pondermotive  wave  after  it  has 
trapped  the  beam  particles.  The  decrease  in  particle  energy  appears  as  an 
increase  in  wave  energy. 

To  operate  a  high  current  beam  in  the  Saman  regime  requires 


where  u>.  is  the  beam  plasma  frequency  (Sprangle  et  al.  ,  1979).  Combining 
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This  condition  can  be  shown  to  be  equivalent  to  requiring  X"  >  X'  ,  where 

X'  and  X'  are  the  wavelength  and  Debye  length  in  the  beam  frame  (Hasegawa, 

1978) .  The  energy  spread  Ay/y  which  results  from  the  normalized 

emittance  e  of  a  beam  with  radius  rK  is  (Neil,  1977) 
n  d 
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The  Lawson-Penner  relation  which  relates  the  emittance  to  the  beam  current  is 


e  2  =  52I(kA)  (cm  -  rad)  (13) 

n 

where  S  is  a  scale  factor  that  is  typically  0.1  -  0.3  for  existing 
accelerators . 


Combining  these  equations  we  find  for  a  Raman  interaction,  that 


X(cm) 


>  5.7 


2 

( kA/cm  )  . 


(14) 


At  a  current  density  of  1  kA/cm^  and  S  =  0.3  the  wavelength  must  be 
greater  than  0.3  cm  to  operate  in  the  Raman  regime  when  y  =  3.  To  operate  a 
1  micron  FEL  in  the  Raman  regime  for  such  a  beam  would  require  y  *  2600,  or 
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1.3  GeV,  and  a  higgler  with  a  wavelength  of  135  meters.  Consequently,  the 
Raman  FEL  is  not  viable  in  the  IR  region  unless  S  can  be  greatly  reduced. 

It  is  clear  from  Sq.  14  that  there  are  three  parameters  one  can  improve 
upon  to  extend  the  Raman  region.  T5ie  kinetic  energy  of  the  beam  can  be 
increased,  the  current  density  decreased  or  the  scale  factor  in  the  eroittance 
relation  decreased.  However,  there  is  a  limit  to  the  extent  one  can  increase 
the  energy,  decrease  the  current  density  and  maintain  a  collective 
interaction.  When  the  number  of  particle  in  the  debye  sphere  is  anall,  the 
system  will  not  support  collective  oscillations. 

The  wavelength  condition  of  Eq.  14  is  most  sensitive  to  S,  so  it  is 
important  to  analyze  the  degree  to  which  S  can  be  reduced.  The  Lawson- Penner 
relation  is  not  derived  from  first  principles  but  is  a  phenomenological 
relation  which  has  been  shown  to  hold  in  many  accelerators.  The  success  of 
the  Lawson-Penner  relation  is  related  to  the  fact  that  the  reliable  current 
density  from  oxide  thermionic  cathodes  is  generally  <  13  A/cm2.  With  higher 
current  density  cathodes,  the  scale  factor,  S,  can  be  reduced  considerably. 

In  a  recent  experiment  using  cold  graphite  cathodes,  an  S  parameter  of  0.12 
•was  obtained  with  a  14  kA  beam  by  aperturing  the  beam  to  about  30  percent  of 
the  diode  current  (Sloan  et  al,  1932).  An  S  parameter  of  .13  at  4  kA  was 
obtained  by  aperturing  to  about  25%  of  the  current  from  a  spiral  knife-edge 
cold  cathode  in  a  20  MeV  induction  li.nac  (Xulke  et  al.  ,  1981).  Finally,  a 
much  reduced  value  of  S  -  0.04  at  1.5  kA  -was  calculated  for  a  recent 
experiment  in  which  only  -  5-10%  of  the  current  from  a  graphite  cold  cathode 
was  extracted  from  the  diode  (Parker  et  al.,  1982). 
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Wie  S  parameter  is  best  understood  in  terms  of  the  brightness,  another 


invariant  that  is  related  to  the  emittancet 
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Thus,  the  brightness  is  a  measure  of  our  ability  to  make  a  high  current 
density,  cold  beam.  To  the  extent  that  the  emittance  is  an  invariant,  the 
brightness  is  determined  at  the  diode;  hence  if  we  compress  the  beam  at  a 
later  stage  we  increase  Ay/Y-  Therefore,  for  high  power  operation  it  is 
desirable  to  have  a  high  current  density  cathode.  From  the  above  three 
experiments  we  can  conclude  that  the  brightness  of  cold  cathodes  with 
appropriate  aperturing  is  a  factor  of  six  to  sixty  higher  than  that  of  hot 
cathodes.  This  will  extend  the  Raman  wavelength  region  by  a  corresponding 
amount.  However,  the  price  one  pays  for  cold  cathode  operation  is  diode 
closure.  This  can  limit  the  pulse  length  to  a  few  hundred  nanoseconds.  In 
addition,  the  efficiency  of  the  system  is  reduced  by  a  factor  of  5  to  10  when 
a  large  fraction  of  the  beam  is  terminated  by  aperturing. 

The  emittance  is  often  the  dominant  contributor  to  energy  spread. 

However,  the  self  potential  and  the  wiggler  gradient  can  also  contribute  to 
the  beam  energy  spread,  which  in  general  is  given  by 
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where  the  second  term  is  the  contribution  from  the  self  potential  and  the 


third  is  from  the  wiggler  gradient.  From  this  we  can  find  an  optimum  radius 
of 
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For  high  current  beams  propagated  in  a  strong  axial  magnetic  field,  the 
equilibrium  radius  (Neil,  197  7)  should  be  compatible  -with  the  optimum  radius. 

The  variation  in  y  as  determined  by  the  emittance,  wiggler  gradients  and 
self  field  of  the  beam  results  in  a  velocity  spread  of  the  beam  that  is 
independent  of  time.  If  the  diode  voltage  is  time  dependent  (as  for  example 
due  to  diode  closure) ,  this  will  shift  the  output  wavelength.  If  the  gain  is 
sufficiently  high  in  an  oscillator  the  time  variation  •'ill  a  broad  band 

output.  The  time  variation  in  an  amplifier  experiment  can  cause  a  loss  of 
resonance  with  the  input  signal.  The  range  of  unstable  wave  numbers  is 
(Sprangle  et.  al. ,  1979) 


Ak  *  3  k.2/k  /  (IB) 

i 

where  k^  is  the  spatial  growth  rate.  This  can  be  related  to  the  variation  in 

2  iv  2  2 

Y  through  X  *  \w/2Y  *  Then  — <  4  A  is  required.  Substituting 
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the  maximum  growth  rate  value  for  k^  we  get 
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Therefore,  the  requirements  on  the  time  variation  of  the  beam  energy  to 

keep  the  growth  of  the  output  wavelength  in  the  unstable  region  is 
3  1 

proportional  to  X  .  For  a  1  kA/cm  ,  1  MeV  beam  in  a  1  kG  wiggler  field  an 
output  wavelength  of  0.7  mm  requires  a  A y/y  time  variation  of  less  than  1%. 

V.  Accelerator  Development 

During  the  past  fifty  years  accelerators  have  developed  primarily  along 
two  directions:  (1)  low  current,  high  voltage  and  (2)  high  current,  low 
voltage  devices.  Consequently  FEL  research  has  followed  similar  lines.  The 
history  of  traditional  accelerators  is  quite  well  known.  lhe  average  current 
in  these  accelerators  is  generally  less  than  an  ampere  with  micro  bunch 
currents  reaching  tens  of  amperes.  The  energies  acheived  in  electron 
linacs  exceeds  10 10  eV  and  the  beam  energy  in  a  proton  synchrotron  is 
approaching  10  *  eV.  High  current  accelerators  have  been  developed  primarily 
as  Marx  generators  driving  a  pulse  forming  line,  which  is  connected  to  a 
diode.  The  energy  of  these  devices  is  typically  around  1  MeV,  although  10  MeV 
devices  have  been  built.  The  pulse  current  from  a  pulse-power  line  diode  has 
reached  10  MA.  These  devices  do  not  scale  well  in  energy  because  the  total 
energy  is  achieved  in  a  single  gap,  hence  the  stress  on  the  insulator  becomes 
excessive  at  high  voltages.  A  free  electron  laser  requires  high  voltage  for 
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short  wavelength  operation  and  high  current  for  efficiency  and  high  gain. 

Only  the  induction  linac  has  been  a  serious  candidate  as  a  FEL  driver  in  this 
parameter  range.  Experiments  at  Lawrence  Livermore  Laboratory  on  the 
Experimental  Test  Accelerator  (ETA)  have  acheived  currents  approaching  10  kA 
and  beam  energies  of  5  MeV.  A  scaled  up  version  of  this  accelerator  is 
designed  for  50  MeV,  10  kA,  30  ns  operation.  The  current  from  this  device  is 
probably  excessive  for  efficient  FEL  operation  and  the  pulse  length  too  short. 

A  long  pulse  induction  linac  (t  *  2  usee)  was  built  at  the  National 
Bureau  of  Standards  (Leiss,  et  al* ,  1980)  and  is  currently  in  operation  at  the 
Naval  Research  Laboratory.  Figure  2  is  a  schematic  of  the  Long  Pulse 
Induction  Linac.  The  accelerator  consists  of  two  major  components;  (1)  an 
injector  and  (2)  an  induction  accelerator  module. 

The  electron  gun  in  the  injector  has  a  16.5  cm  diameter  tungsten  dispenser 
thermionic  cathode.  Electrons  are  accelerated  in  the  gun  through  a  series  of 
12  annular  electrodes,  spaced  by  ceramic  insulator  rings.  The  last  electrode 
supports  a  95%  transmission  tungsten  mesh  at  ground  potential.  The  gun  has 
also  been  operated  with  a  graphite  brush  cold  cathode  (Ramirez  and  Cook,  1980; 
Prohaska  and  Fisher,  1982) . 

The  electron  gun  is  immersed  in  an  oil  filled  tank.  The  gun  voltage  is 
fed  from  a  pulse  line  driving  a  12:1  step  up  transformer.  ftae  injector 
typically  produces  a  0.8  kA  beam  pulse  of  400  keV.  The  electron  beam  is 
transported  to  the  induction  accelerator  module  by  a  series  of  focusing  coils. 
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INDUCTION  LINAC 


Figure  2.  Schematic  of  Long  Pulse  Induction  Linear  Accelerator. 


The  induction  accelerator  module  consists  of  two  core  sets.  One  of  the 
core  sets  give  a  4  to  1  step  up  voltage  and  the  other  a  5  to  1.  The  cores  are 
wound  with  0.001  inch  mild  steel  foil,  separated  by  0.00025  inch  mylar 
sheets.  Each  accelerating  gap  produces  about  200  kV  of  acceleration. 

Typically,  the  output  energy  of  the  electron  beam  generated  by  the  linac 
is  0.3  MeV,  the  current  approximately  0.8  kA  and  the  pulse  length  2  ysec. 
However,  for  the  data  reported  in  this  work  one  of  the  accelerating  gaps  was 
not  operating,  so  the  beam  voltage  was  0.55  MeV.  Also,  a  graphite  brush  cold 
cathode  was  used  and  the  0.3  kA  diode  current  was  reduced  to  0.2  kA  due  to 
losses  in  beam  transport.  The  temporal  variation  in  the  voltage  with  the  hot 
cathode  is  less  than  3  percent  over  1.6  ysec.  This  is  the  longest  pulse 
induction  linac  in  existance,  with  a  pulse  duration  more  than  an  order  of 
magntiude  longer  than  other  induction  linacs.  The  pulse  length  becomes  an 
important  consideration  for  free  electron  laser  experiemnts  where  one  wishes 
to  study  the  nonlinear  dynamics  of  the  beam  or  efficiency  enhancement 
schemes.  Applications  that  require  a  significant  amount  of  energy  in  the 
radiation  field  also  require  long  pulses  to  avoid  the  problems  encountered 
with  excessive  electric  field  strengths  at  short  pulses. 

Figure  3a  shows  the  voltage  and  current  traces  of  the  injector  for  the 
long  pulse  induction  linac.  The  voltage  remains  constant  to  within  a  few 
percent  over  the  2  microsecond  pulse  length.  In  contrast  Figure  3b  shows  the 
voltage  and  current  from  the  graphite  cathode  diode  of  the  Febetron,  a 
relatively  long  pulse  Marx  type  generator.  The  oscillations  on  the  voltage 
and  current  are  due  to  the  finite  element  pulse  line  of  the  Febetron  Marx. 

The  diode  voltage  decreases  25  percent  in  250  nsec  as  a  result  of  diode 


closure.  This  closure- induced  voltage  collapse  is  typical  of  cold,  plasma 
cathodes,  and  in  light  of  the  condition  of  Sq.  19,  is  a  key  problem  for  long 
pulse  FEL  operation.  Efforts  to  prevent  diode  closure,  such  as  magnetic 
insulation,  generally  result  in  higher  beam  emittance. 

The  relationship  between  the  cathode  plasma  velocity  and  the  diode 
parameters  can  be  determined  by  considering  the  simple  circuit  shown  in  Fig. 
4.  There  is  a  voltage  source  Vc,  an  internal  line  impedance  Z^,  and  a  time 
dependent  diode  impedance  ZD(t).  The  current  i3 


I  -  VC/(ZL  +  z0) . 


(20) 


The  diode  voltage  is 


VD  "  IZD  *  W1Zl  +  ZD> 


(21) 


and  the  diode  inqpedance  is 


,  _  ad2(t) 

ZD  1/2  ' 

D 


(22) 


where  a  is  a  constant  and  d  the  anode  cathode  spacing.  Me  will  assrne 


d-d  -  at  (23) 

o 

where  a  is  the  plasma  closure  velocity.  The  closure  velocity  for  high  current 
cold  cathodes  is  typically  -  2.5  cm/us.  Now  if  we  differentiate  Eq.  (21) 
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with  repsect  to  time  and  asssume  the  diode  impedance  is  matched  to  the  line  at 
t  ■  o,  then  vc  *  2  VQ  and 


AV_ 


AZ 

1  _2 

2  Z_ 
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Then 


and 
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4  Ad 
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VQ(t)  -  V  (o)  (r) 
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4/5 


VQ(0)  (l  -  at/dj 


4/5 


(26) 


(27) 


Hence  for  a  diode  that  is  initially  matched  we  can  obtain  a  simple  expression 
for  the  time  dependence  as  a  function  of  diode  spacing. 

Figure  5  is  a  plot  of  the  diode  voltage  as  a  function  of  time  for  Z.  *  ZQ, 
Z»  »  10  Z0  and  Z^  *  0.1  ZQ  when  t  *  o.  When  the  line  impedance  is  ten  times 
the  diode  impedance,  the  device  tries  to  behave  as  a  constant  current  source 
and  the  voltage  decays  faster  than  in  the  matched  load  case.  When  the  line 
inpedance  is  small  compared  to  the  diode  impedance,  the  device  tries  to  act  as 
a  constant  voltage  source  at  the  expense  of  an  increased  current.  However, 
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T*at/d0 

Figure  5.  Calculated  diode  voltage  and  current  for  three  different  diode 

impedances.  The  line  impedance  is  and  ZQ  is  the  diode  impedance 
at  t  *  0. 
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for  high  currant  beams,  the  increasing  self  potential  of  the  beam  with  current 
will  cause  the  beam's  kinetic  energy  to  decrease  faster  than  the  diode 
voltage. 

Although  diode  closure  in  high  current  cold,  cathode  machines  is  a 
limiting  factor  for  long  pul3e  operation,  there  are  many  operational 
advantages  to  cold  cathodes.  Thus,  we  have  tenporarily  replaced  the  hot 
cathode  in  the  induction  linac  with  a  graphite  brush  cold  cathode  of  the  same 
area.  The  average  electric  field  at  the  cathode  surface  is  about  30  kV/cm. 

The  cathode  rapidly  turns  on,  and  the  closure  problem  does  not  seem  to  be  as 
severe  as  in  the  Febetron.  Perhaps  the  large  area,  low  current  density 
reduces  the  closure  velocity.  The  emittance  is  worse  than  with  the  hot 
cathode,  and  more  of  the  beam  is  lost  in  transport.  However,  the  emittance  of 
the  beam  transported  to  the  wiggler  appears  to  be  about  the  same  as  with  the 
hot  cathode. 

VI.  The  FEL  Apparatus 

Figure  6  is  a  sketch  of  the  free  electron  laser  experiment.  There  is  a 
uniform  axial  field  of  120  cm  length  which  is  varied  from  1  to  5  kG,  with  a 
typical  operating  field  of  2  kG.  'Hie  beam  from  the  induction  linac  is  focused 
into  the  solenoid.  The  experiment  has  not  yet  been  run  with  the  mirrors  in 
place  as  shown,  but  operated  in  a  superradiant  mode. 

Two  wiggler  configurations  have  been  investigated:  (1)  a  pulsed  linear 
wiggler  and  (2)  a  radially  symmetric  diffusive  wiggler.  The  amplitude  of  the 
linear  wiggler  field  can  be  varied  from  0.1  kG  to  1  kG.  The  wiggler 


wavelength  in  3.0  cm,  and  the  overall  length  of  the  wiggler  is  120  cm.  The 
wiggler  amplitude  rises  adiabatical ly  in  30  cm,  has  a  uniform  straight  section 
of  50  cm  and  decays  adiabaticallv  in  the  last  30  cm. 

Two  diffusive  wigglers  have  been  used:  a  6  cm  period  aluminum  wiggler  and 
a  4.5  cm  copper  wiggler.  The  perturbed  axial  component  of  the  on-axis  field 
is  4  percent  of  the  axial  field  for  the  6  cm  wiggler  and  6  percent  for  the  4.5 
cm  wiggler. 

The  output  radiation  is  extracted  through  a  large  area  window,  and  the 
power  and  spectrum  are  determined  by  gas  breakdown  thresholds  or  with  high 
pass  filters  and  calibrated  detectors. 

A.  Linear  Wigglers 

There  are  several  advantages  to  linear  wigglers,  including  ease  of 

assembly,  changing  the  periodicity  to  operate  at  different  wavelengths,  and 

tapering  the  period  and/or  field  amplitude  for  efficiency  enhancement 

( Sprangle  et  al. ,  1979).  As  a  result  of  our  experiments  and  analysis,  we  find 

that  a  linear  wiggler  in  an  axial  guide  field  produces  an  elliptical 

polarization.  The  ratio  of  the  major  to  minor  axis  of  the  ellipse  is  k  v./f2  , 

w  z  o 

where  k  is  the  wiggler  wave  nanber,  v  the  axial  beam  velocity,  and  2  the 

*  2  o 

cyclotron  frequency  in  the  axial  guide  field.  As  a  result  of  the  asymmetry  of 
the  wiggler,  there  are  no  focusing  forces  in  one  of  the  directions 
perpendicular  to  the  beam  propagation  and  an  electron  drift  results.  Although 
this  drift  can  be  very  small  in  some  parameter  regimes,  it  is  always  present 
in  linear  wigglers  immersed  in  an  axial  guide  field  -unless  additional  focusing 
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forces  are  provided.  This  problem  was  avoided  in  a  high  current  device  using 
a  linear  wiggier  by  eliminating  the  solenoidal  field  (Phillips,  1360:. 
Focusing  was  provided  by  increasing  the  wiggier  field  at  the  edges. 

To  understand  how  the  elliptical  polarization  and  perpendicular  drift 
arise,  we  will  solve  the  equations  of  motion  in  the  combined  wiggier  and  axia 
fields.  The  linear  wiggier  field  components  are 


3  =  b(z)  3  cosh  k  v  cos  k  z  (23) 

wy  w  w'  w 


3  =  -  b(z)  3  sinh  '<  v  sin  k  z 

wz  w  w*  w 

where  3W  is  the  peak  wiggier  field  on  axis  and  b(z!  is  an  adiabatic  taper  m 
the  wiggier  amplitude  at  the  entrance  and  exit  of  the  wiggier.  The  amplitude 
increases  to  its  full  v< lue  in  ten  periods  in  our  experiment. 


We  define  Cl  =  qB  /ym  and  1  =  aS  /ym.  To  simplifv  the  analvsis  we  will 

o  o  v  *  w 

make  the  following  assumptions 
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3  ' 
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V 
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k  y  <  <  1 
w 


(29) 


and  v 

z 


v,  a  constant, 
o 


The  single  particle  equation  of  motion  then  gives 
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V  = 
X 


(ft  v  -ft  cosh  k  y  cos  k  z  v  ) 
o  y  w  w  w  z 


V  =  -V  0, 

y  x  o 


(  30) 


v  =  v  ft  cosh  k  y  cos  k  z 
z  x  w  w  w 


where  we  have  neglected  the  correction  to  3Z  from  the  wiggler  since  3W/3Q 


<< 


and  v  <<  v  . 

y  z 


We  take  the  derivative  of  each  of  the  components  with  respect  to  time  and 
use  the  assumptions  of  Sq.  29  to  get 


"2  2 
v  +  ft  v  =  ft  k  cosh  k  y  sin  k  z  v 
X  ox  w  w  w  w  z 


Similarly 


v  +  ft  v  =  ft  ft  cosh  k  y  cos  k  z  v 
y  o  y  o  w  w  w  z 


We  have  the  following  solutions 


ft  k  cosh  k  v  sin  k  z  v 
W  W"  W  z 


_  2  .2  2 
ft  -  k  v 

o  w  z 


ft  ft  cosh  k  v  cos  k  z  v 
o  w  v*  w  z 


-  2  .2  2 
ft  -  k  v 
o  w  z 


ft  v  cosh  k  y  cos  k  z 


(31) 


(32) 


(33) 
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Q  ft  cosh  k  y  sin  k  z 
o  w _  w  w 


w 
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The  displacement  in  the  x  direction  is  proportional  to  the  wiggler  field 
and  axial  velocity.  TSie  displacement  in  the  y  direction  occurs  only  if  there 
is  a  nonzero  solenoidal  magnetic  field.  It  is  proportional  to  both  the 
solenoidal  and  higgler  field.  This  displacement  is  a  result  of  the  Lorentz 
force  arising  from  the  x-coraponent  of  vw  and  BQ. 


The  ratio  of  the  maximum  excursion  in  the  x-direction  to  that  in  the  y- 
direction  is  seen  from  Sq.  33  to  be 


Ax/Ay  *  k  v  /Q  .  (34) 

w  z  o 

Thus,  this  ratio  is  a  measure  of  how  close  S3  is  to  the  cyclotron  resonance 

o 

which  occurs  at  =  kw*/2.  Since  k^vz  is  approximately  constant  in  the 

experiment,  the  ratio  of  the  maximum  displacements  Ax/Ay  should  vary 

as  1/fi  . 
o 

Figure  7  is  a  series  of  exposures  made  from  x-rays  produced  when  the  beam 
strikes  a  target  placed  beyond  the  wiggler.  These  photos  give  us  the  shape  of 
the  beam  versus  solenoidal  field.  In  Fig.  7(a),  the  wiggler  field  is  zero  and 
the  beam  is  approximately  circular.  In  Fig.  7(b),  with  a  500  Gauss  wiggler 
field  the  beam  becomes  elongated  in  the  x-direction.  The  linear  wiggler  field 
is  in  the  y  direction,  consistent  with  the  above  analysis.  As  the  solenoidal 
field  13  increased  with  the  wiggler  field  constant,  the  shape  goes  from 
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Bq  s  2.6  Byy  s  0.5 

(0 


B0  *  3.48  Bw  *  0.5 

(d) 


gure  7. 


Time  integrated  x-ray  exposure  of  beam  striking  a  target  placed 
beyond  linear  wiggler  (wiggler  field  is  vertical)  for  various 
values  of  solenoidal  field. 
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elliptical  to  circular,  as  shown  in  Figs.  7(b)  -  7(d).  Hence,  we  would  expect 
a  nearly  linearly  polarized  output  for  (b)  and  a  nearly  circular  polarization 
for  (d) . 

Figure  8  is  a  plot  of  the  ratio  Ax/Ay  as  a  function  of  the  solenoidal 
field.  The  results  are  plotted  three  ways.  The  solid  line  is  the  theory,  the 
small  dots  were  computed  numerically  and  the  circled  x  points  are 
experimental.  At  small  values  of  magnetic  field  the  approximations  break  down 
and  there  is  a  large  variation  in  the  results. 

We  have  demonstrated  that  the  polarization  of  a  linear  wiggler  in  a 
solenoidal  magnetic  field  may  be  varied  by  changing  the  magnitude  of  the 
solenoidal  field.  In  addition,  we  see  from  Eqs.  34  that  a  change  in  the 
amplitude  of  the  wiggler  field  alone  does  not  change  the  polarization.  We 
have  confirmed  this  by  propagating  a  beam  the  entire  length  of  the  wiggler  in 
a  solenoidal  magnetic  field  of  2  kG.  The  wiggler  field  is  varied  up  to  1  kG 
without  changing  the  shape  or  intensity  of  the  beam. 

As  a  result  of  the  asymmetry  of  the  wiggler,  there  are  no  focusing  forces 
in  the  x-direction.  This  results  in  a  net  drift  of  off  axis  particles  (Pasour 
et  al.,  1982).  This  drift  is  shown  clearly  in  Fig.  9,  which  consists  of 
electon  trajectories  in  the  x-y  plane.  The  trajectories  are  calculated  by 
numerically  solving  the  equations  of  motion,  including  self-fields.  Initially 
the  electrons  S  x  B  drift  azimuthally  as  they  are  injected  into  to  the 
adiabatically  increasing  wiggler  field.  However,  when  they  reach  the  constant 
amplitude  portion  of  the  wiggler,  the  electrons  drift  rapidly  in  a  direction 
perpendicular  to  both  the  axial  and  wiggler  fields. 
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Figure  3.  Plot  of  bean  elongation  (dx/Ay)  vs.  axial  magnetic  field  in  a 
linear  wiggler  (B^  in  y-direction) .  Solid  line  is  theoretical 
curve  (Eq.  34) ,  dots  are  values  obtained  by  numerical  integration 
of  orbit  equations,  and  crosses  are  experimental  results. 


Figure  9.  Electron  trajectories  in  a  linear  wiggler  and  an  axial  magnetic 


field.  In  this  case,  B  =  1  kG,  B  *  250  A,  X  =  3  cm,  and 

WO  w 

Y  -  2.2. 
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An  approximate,  empirical  expression  for  the  drift  has  been  found  which  is 
in  good  agreement  with  results  from  computer  calculations  of  the  electron 
orbits  when  <  kv^.  In  terms  of  the  wiggler  gradient,  the  expression 

is 


-  i  Ll 

vd  *  2  3 


B 


7(B)' 

w  o 

+  2  x  B 

2(B  )  o 

w 


(35) 


2  2  "> 

where  v,  vx  v  •  This  expression  has  the  same  form  as  the  usual  grad  B 
drift  and  can  be  thought  of  as  arising  from  the  crossed  wiggler  gradient  and 
axial  field.  It  can  also  be  written  as 


,  x  v 

1  W  2 

v  JI  —  v  — - - 

d  2  z  2 

o 


ffl  k  v  )2  cosh  k  y  sinh  k  y 
v  w  w  z  w^  w 


(36) 


ffl  2-  !c  2v  21 

v  O  W  Z  ' 


This  expression  is  reasonably  valid  for  k^  £  0.8,  Bwy  >>  Bwz,  and 
v,  <<  v^.  Mote  that  the  drift  is  very  small  near  the  axis  but  increases 
exponentially  with  '^y. 

Table  1  compares  drift  velocities  for  various  cases  as  calculated  by  the 
computer  trajectory  code  to  those  calculated  from  the  Be.  36.  In  general  the 
agreement  is  very  good.  One  consequence  of  the  large  drift  near  resonance  is 
that  it  limits  the  degree  of  gain  enhancement  achievable  through  the  magneto¬ 
resonance  effect. 
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Comparison  of  drift 

velocities  computed  from 

the  code 

with  those 

from  Eq.  9  for 

various  cases. 

v  . 

v  . 

--  (Eq-  9) 

B0(kG)  Sw(kC) 

0 

z 

X(cra) 

—  (code) 
c 

2 

1 

2.2 

0.29 

0.4 

0  .047 

0  .050 

4 

1 

2.2 

0.57 

0.4 

0  .042 

0  .046 

2 

0.5 

2.2 

0.29 

0.4 

0  .011 

0  .013 

2 

1 

2.2 

0.29 

0.2 

0.019 

0  .022 

2 

1 

3.0 

0.21 

0.4 

0.029 

0  .034 

4 

1 

3.0 

0.42 

0.4 

0.019 

0  .022 

4 

0.5 

3.0 

0.42 

0.4 

0.0051 

0  .0056 

4 

1 

10.0 

0.11 

0.4 

0.0043 

0  .0049 

LO 

5 

10.0 

0.28 

0.4 

0.057 

0  .056 
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3.  Helical  Wiggler 


Because  the  linear  wiggler  drift  is  so  sensitive  to  slight  off-centering 
or  large  diameter  of  the  beam,  a  helical  wiggler  would  seem  to  be  somewhat 
better  suited  to  the  experiment,  at  least  initially  when  beam  focusing  is  not 
optimized.  If  the  trajectory  calculation  is  repeated  -with  the  3ame  parameters 
as  those  in  Fig.  9  but  with  a  helical  wiggler  substituted  for  the  linear  one, 
it  is  found  that  the  electron  drifts  azimuthally  and  is  well  confined,  as 
shown  in  Fig.  10  (Pasour  et  al. ,  1982}.  This  behavior  results  from  the 
radially  increasing  wiggler  field,  which  for  k  r  <  0.8  can  be  written  as 
( 3lewett  and  Chasman,  1977) 


3  »  B  (l  +  -^k2r2)  sin  (9  -  k  z) 

r  w  8  w  w 


B-  »  3  (1  +  x  k  2r2)  cos  (9  -  k  z)  (37). 

0  w  8  w  w 


1  2  2 

B  «  -  krB  (L  +  —  k  r  )  cos  (9  -  k  z)  . 
z  w  3  w  w 


Although  a  helical  wiggler  is  in  general  more  difficult  to  construct  than 
a  linear  one  and  is  not  nearly  as  amenable  to  tapering,  the  orbital  stability 
that  it  provides  is  a  major  advantage  when  an  axial  guide  field  is  used. 

Also,  it  is  possible  to  taper  the  amplitude  of  the  wiggler  field  in  the 
helical  wiggler  by  carefully  varying  the  winding  radius  as  a  function  of  z. 
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Figure  10.  Electron  trajectories  in  a  helical  wiggler  and  an  axial  magnetic 

l  field.  Parameters  are  identical  to  those  in  Figure  9. 

1) 
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C.  Diffusive  Wiggler 


An  axially  symmetric  wiggler  can  be  constructed  by  placing  conducting 
rings  in  a  pulsed  solenoidal  magnetic  field  !  Jacobs  et  al.  ,  1980).  .As  the 
solenoidal  field  diff fuses  in,  it  induces  eddy  currents  in  the  conducting 
rings  which  generate  an  opposing  field.  This  results  in  a  magnetic  field  with 
a  periodic  radial  component.  'Si is  process  is  illustrated  in  the  computer 
plots  of  Fig.  11. 

The  magnetic  field  at  a  distance  r  from  the  the  axis  is  approximately 


3  =  z  f3 
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1  o  w 
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where  IQ  and  1^  are  modified  Bessel  functions  and  is  the  anqplitude  of  the 
axial  field  modulation  on  axis. 

Figure  12  is  a  computer  plot  of  the  axial  component  of  the  magnetic  field 
of  the  diffusive  wiggler  used  in  the  system.  When  the  wiggler  is  not  present, 
the  solenoidal  field  is  2  kG.  Hence,  the  wiggler  not  only  modulates  the 
solenoidal  field,  it  reduces  the  average  value.  To  make  a  smooth  transition 
into  the  wiggler  field,  the  conducting  rings  extend  all  the  way  to  the  end  of 
the  solenoid  where  the  field  drops  to  one  half  the  peak  value. 

Since  the  diffusive  wiggler  is  axially  symmetric,  it  does  not  produce  the 
kind  of  radial  particle  drift  which  is  characteristic  of  linear  wiggiers  in  an 
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wiggler.  Plots  show  field  lines  at  progressively  later 
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Axial  Magnetic  Field  (Gauss 


Axial  Position 

Figure  12.  Plot  of  B,(r=0)  vs.  2  in  the  diffusive  wiggl 
is  at  1  =  0. 
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axial  field.  However,  the  perpendicular  component  of  the  wiggler  vanishes 
along  the  axis  of  the  diffusive  wiggler,  so  that  only  electrons  that  are 
significantly  off  axis  can  participate  in  the  type  of  FEL  interaction 
described  previously.  It  has  been  pointed  out,  however,  that  an  interaction 
that  relies  on  the  transverse  energy  in  the  bean  can  be  exploited  to  obtain 
radiation  from  such  a  wiggler.  This  device  has  been  called  the  Lowbitron 
(McMullin  and  3ekef i,  1981).  The  radiation  frequency  is  approximately 

u)  =»  (1  +  3  2)  y  2  (v  k  +  2  /y)  .  (39) 

2  Z  2  W  O 

This  radiation  is  a  result  of  electrons  interacting  with  the  axial 
component  of  the  wiggler  field.  The  Lowbitron  interaction  requires 
k  c3,Y /2  <  1.  The  electron  gyroradius  must  be  sufficiently  small  that  the 

transverse  field  modulation  felt  by  the  electron  can  be  neglected  co-oared 
with  the  longitudinal  modulation. 

VII.  Radiation  Measurements 

We  have  carried  out  a  series  of  measurements  on  the  radiation  using  a 
diffusive  wiggler.  Radiation  measur ements  have  primarily  consisted  of  a) 
spectral  analysis  using  cylindrical  cut-off  filters  and/or  a  gas  breakdown 
spectrometer  and  b)  approximate  power  measurements  using  either  pressure 
thresholds  for  gas  breakdown  or  calibrated  crystal  detectors  and 
attenuators.  Careful  consideration  has  been  given  to  the  identification  of 
the  interaction  modes,  both  cyclotron  and  FEL,  and  scaling  measurements  have 
been  performed  to  verify  these  modes. 
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A.  Gas  Breakdown  Spectrometer 


The  radiated  power  can  be  estimated  from  the  pressure  at  which  gas 
breakdown  occurs.  An  empirical  formula  for  the  rms  breakdown  electric  field 

%  is 

E,  =  A?fl  +  (a/PX)2]1/2,  (4  0) 

D 

where  A  =  3000  Vm-~Torr-'''  and  a  =  0.9  Torr-tn  for  air  or  nitrogen,  P  is  the  gas 
pressure  in  Torr,  and  X  is  the  wavelength  in  meters.  Equation  (40)  has  been 
checked  against  McDonald's  gas  breakdown  data  (McDonald,  1966)  and  agrees  in 
the  'worst  case  to  within  30%.  This  formula  is  valid  on  the  high  pressure  side 
of  the  breakdown  curve  and  for  X  <  30  cm,  pulse  lengths  >1  us,  low  repetition 
rate  (<  100  pps) ,  and  electron  diffusion  lengths  small  compared  to  the  chamber 
dimensions. 

If  the  radiation  is  reflected  from  a  boundary,  a  standing  wave  is  set  up 
with  an  amplitude  given  by 

E2  =  S  2  fl  +  T2  -  2T  cos  (2  k  •  x] ;  (41) 

where  E^  is  the  magnitude  of  the  incident  wave,  T  the  reflection  coefficient, 
and  k  the  wavenumber.  The  distance  d  between  the  peaks  in  the  standing  wave 

■>  4 

is  determined  by  setting  2  k  •  x  «*  2t ,  so  that 


d 


X 

2  cos  9 


(42) 
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where  X  is  the  wavelength  of  the  radiation  and  3  the  angle  of  incidence  with 


respect  to  the  surface  normal.  3y  adjusting  the  gas  pressure  so  that  the 
standing  wave  field  amplitude  is  slightly  higher  than  the  breakdown  field, 
localized  gas  breakdown  will  occur  at  the  standing  wave  peaks.  Then  X  can  be 
determined  by  measuring  the  distance  between  adjacent  breakdown  spots,  and  the 
power  can  be  estimated  from  the  pressure  required  to  initiate  observable 
breakdown. 

Figure  13  shows  a  schematic  of  the  spectrometer.  Radiation  is  collected 
in  the  horn  on  the  left.  It  then  travels  through  a  high  pass  filter  and 
expands  through  the  horn  on  the  right  to  a  collimating  or  weakly  focusing 
lens.  The  microwaves  are  reflected  from  a  metallic  boundary  located  inside 
the  gas  filled  chamber. 

Figure  14  is  a  time  integrated  photograph  of  the  resulting  interference 
pattern.  Microwave  radiation  has  entered  from  the  left  of  the  photograph 
where  the  lens  is  located  (but  cannot  be  seen)  and  is  reflected  off  a  copper 
plate  located  at  the  right  of  the  photograph  (outlined  by  a  light  ring) .  All 
the  white  spots  are  due  to  light  produced  when  the  gas  breaks  down.  Type  57 
(ASA  3000)  Polaroid  film  was  used  with  a  Graphflex  camera  (f/4.5)  to  obtain 
these  results.  The  gas  density  was  selected  so  that  enough  light  was 
available  for  photographing  the  spots,  but  was  kept  low  enough  to  prevent 
microwave  reflections  from  the  plasma.  Thus  for  the  measured  pressure  (25 
Torr)  and  wavelength  (4.5  cm),  the  electric  field  is  determined  from  Eq.  40  to 
be  about  1  kV/cm.  From  the  spot  diameter  of  2.5  cm,  the  power  is  estimated 
for  a  plane  wave  to  be  3  kw. 
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Figure  14.  Open-shutter  photograph  of  gas  breakdown  in  25  Torr  nitrogen. 

From  the  2.25  cm  spot  spacing,  the  radiation  wavelength  is  4.5  cm. 
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This  technique  underestimates  the  power  because  there  is  considerable 
energy  outside  the  radius  of  the  observed  spot.  For  example,  in  the  case 
above  the  spot  diameter  is  nearly  a  factor  of  2  smaller  than  the  wavelength, 
so  clearly  the  radiation  envelope  is  larger  than  the  measured  spot  size. 

We  have  identified  this  radiation  as  a  TEq^  mode  that  is  excited  at  the 
2nd  harmonic  of  the  electron  cyclotron  beam  interaction.  With  similar 
measurements  we  have  observed  more  than  100  kW  of  power  at  a  wavelength  of  7. 
cm.  This  mode  was  identified  as  a  TE^  mode  excited  at  the  fundamental 
cyclotron  beam  interaction.  The  electron  beam  current  in  these  experiemnts 
was  130  A,  the  voltage  550  keV  for  2  usee,  and  the  wiggler  period  was  6  cm. 
The  gas  breakdown  spectrometer  has  been  described  in  more  detail  elsewhere 
(Mako  et  al.  ,  1982)  . 

3.  Mode  Analysis  and  Experimental  Results 

The  strongest  interactions  are  expected  when  the  phase  velocities  of  the 
beam  modes  and  waveguide  modes  are  equal.  The  dispersion  relation  for  the 
waveguide  mode  is 


2 


to 


2 

0) 

CO 


.  2  2 
+  .<  c 


(43) 


where  u»  is  the  cutoff  frequency.  The  dispersion  relation  for  the  free 
electron  laser  is 


(k 


k  )  v  -  u  /y 
w  z  p 


3/2 


(44) 
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The  last  term  is  small  at  the  current  densities  of  this  experiment  and  is 
neglected.  The  cyclotron  mode  dispersion  relation  is  given  by 

a)  =  k  v  +  n  H  /y,  (45) 

z  o 

where  a  eSc/!n  an<*  'f  have  been  explicitly  written  to  make  the  energy 
scaling  clear. 

Then  the  intersections  of  the  FEL  and  cyclotron  modes  with  a  waveguide 
mode  are  found  from  Sqs .  43-45  to  be 


<jd 


FEL 


k  v  V 

w  z  z 


1  + 


ui 

eye 


nfl  y 
o  z 


1  +  3 
—  z 


(46) 


(47) 


where  w  ■  2+c/aX  ,  a  is  the  waveguide  radius  and  x__  *  3.41  for  TE, , ,  2.61 
co  nm  nm  11' 

for  TM^,  2.06  for  TE23.,  1.14  for  TM1 1  and  1.14  for  TMq2  modes. 

A  sketch  of  three  dispersion  relations  is  shown  in  Fig.  15.  The  cyclotron 

mode  will  be  above  the  FEL  mode  if  nH  /y  >  k  v  .  There  are  typically  both 

o  w  z  c  1 

high  frequency  and  low  frequency  intersections  of  the  FEL  and  cyclotron  modes 
with  the  waveguide  mode.  However,  the  high  frequency  FEL  intersection  is  the 
one  which  results  in  the  usual  FEL  interaction.  Also,  the  gyrotron  typically 
operates  at  the  low  frequency  cyclotron  intersection,  but  a  cyclotron 
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Radiation  growth  is  expected  at  the  intersection 


interaction  i3  also  possible  at  the  high  frequency  intersection.  Figure  16  is 
a  plot  of  the  waveguide  beam  mode  interactions  that  are  possible  for  the 
parameters  of  the  experiment.  Both  the  first  and  second  harmonics  of  the 
electron  cyclotron  interaction  are  shown,  assuming  y^  3  1.6.  The 
justification  for  this  assumption  will  be  given  below.  The  arrows  in  the 
figure  indicate  the  frequencies  at  which  the  strongest  interactions  were 
observed  experimentally. 

To  positively  identify  the  FEL  mode  at  f  *  16  GHz,  we  have  varied  the 

magnetic  field  and  wiggler  wavelength.  It  turns  out  that  the  cyclotron  mode 

is  very  insensitive  to  Bq  in  our  parameter  range.  3ecause  the  magnetic  moment 

is  approximately  conserved,  the  leading  term  for  the  cyclotron  wave  2  y  ^  is 

o  z 

nearly  a  constant  for  small  changes  in  the  magnetic  field.  Thus,  the 

frequencies  of  the  electron  cyclotron  modes  increase  only  slightly  with 

increasing  magnetic  field.  However,  the  FEL  frequency  increases  with  a 

2 

decreasing  magnetic  field  because  y^  varies  as  the  inverse  of  Bq. 

Figure  17  contains  the  results  of  the  mode  identification  experiments. 

The  measurements  were  made  using  high  pass  cutoff  filters,  so  each  bar  is  an 
indication  of  the  resolution.  This  technique  integrates  the  total  power  above 
the  cut  off  frequency  of  the  particular  filter.  Thus,  the  power  in  a 
particular  band  is  just  the  power  measured  with  that  filter  minus  the  power 
measured  with  the  next  smaller  diameter  filter. 

The  arrow  is  the  calculated  FEL  frequency  with  y  *  1.6.  At  \  *  7  cm 

z  w 

and  B^  »  2.45  kG  the  calculated  frequency  is  11.6  GHz  and  this  is  where  we 
see  a  factor  of  two  increase  in  the  signal  amplitude.  When  the  magnetic  field 
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Frequency  (GHz) 

Figure  16.  Beam-waveguide  mode  interactions  which  are  possible  with  the 


experimental  parameters.  The  arrows  denote  the  measured  output 
frequencies . 
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Frequency  (  GHz  ) 

Figure  17.  Frequency  scaling  of  one  particular  mode  with  magnetic  field 
and  wiggler  period,  as  determined  by  using  a  series  of  high- 
pass  cut-off  filters.  Except  as  limited  by  the  detector  band¬ 
width,  the  signal  amplitude  shown  in  each  frequency  interval  is 
proportional  to  the  integrated  power  at  all  frequencies  above 
the  cut-off  of  the  corresponding  filter.  Thus,  power  in  a 
particular  interval  is  indicated  by  a  decrease  in  signal  amplitude 
at  the  next  higher  interval.  The  arrows  locate  the  theoretical 
frequency  of  the  TE^  FEL  mode,  which  agrees  well  with  the 
observed  spectra. 


is  reduced  to  2.17  kG,  the  FEL  frequency  should  increase  to  13.2  GHz,  and  we 


observe  a  factor  of  two  increase  in  that  cutoff  filter.  At  X  »  6  ca  and  3  * 

w  o 

2.17  kG,  the  FEL  interaction  should  occur  at  16.5  GHz  and  we  observe  a  factor 

of  10  increase  in  that  cutoff  filter.  We  have  also  observed  that  the  electron 

cyclotron  mode  at  21  GHz  increases  with  magnetic  field.  We  have  looked  for 

2 

radiation  from  the  Lowbitron  interaction  at  u  »  2y  (2  /y  +  k  v  ) .  This 

z  o  w  z 

interaction  requires  perpendicular  velocity  and  interacts  with  variations  in 
the  z  conponent  of  the  wiggler  field.  We  did  not  see  Lowbitron  radiation  that 
was  comparable  in  amplitude  with  the  FEL  and  2nd  harmonic  cyclotron  radiation. 

In  summary,  we  have  observed  free  electron  laser  radiaton  from  the  long 
pulse  induction  linac.  Tile  FEL  'was  operated  in  the  superradiant  mode.  The 
power  radiated  in  the  FEL  mode  was  approximately  13  kW  with  the  beam 
interacting  with  the  TE^  waveguide  mode.  Second  harmonic  electron  cyclotron 
radiation  of  conparable  amplitude  was  observed  in  the  TE^.  mode.  Most  of  the 
power  resulted  from  the  low  frequency  cyclotron  interaction.  More  than  100  kW 
of  radiated  power  at  the  fundamental  electron  cyclotron  mode  was  observed  in 
the  TE^  mode.  In  all  cases  the  duration  of  the  radiation  was  2  usee.  The 
thermal  spread  of  the  beam  in  the  region  requires  FEL  operation  in  the  kinetic 
Compton  regime. 

The  radiation  spectrum  is  very  sensitive  to  the  average  value  of  y  , 
which  was  determined  by  fitting  the  observed  spectrum  with  interaction 
frequencies  calculated  using  an  assumed  value  of  <yz>.  Then  small  variations 
in  the  magnetic  and  wiggler  period  were  made  to  determine  if  the  mode  'was  an 
FEL  or  cyclotron  mode.  In  this  manner  we  were  able  to  determine  that  <y^> 
had  to  be  1.6  +  .2  . 
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1 


Another  way  of  estimating  y is  by  measuring  the  beam  radius  in  the 
uniform  field  of  the  solenoid.  Since  the  cathode  in  this  experiment  is  in  a 
magnetic  field  free  region,  the  canonical  angular  momentum  ?a  is  zero.  Then 


ymr  v?  =  erAg, 


(48) 


where  Ag  is  the  magnetic  vector  potential  which  for  a  uniform  magnetic  field 
is  Ag  =  r  3^/2  .  Then 


,  eB  r 
3  -  1  — i- 
3  2  ymc 


(49) 


Experimentally,  the  rms  beam  radius  was  measured  from  time  averaged  x-ray 
scintillation  pictures  to  be  1.3  cm  when  3z  =  2.17  kG  and  y  *  2.1.  Thus 
conservation  of  alone  gives  3  ■  0.4  at  this  rms  beam  radius.  Then 


1/2 

f  =  Y/fl  +  B  2  Y2 1  =  1-6 


(50) 


In  addition,  computer  simulations  of  the  induction  li.nac  free  electron 

laser  configuration  have  been  carried  out  (Thompson  et.  al. ,  1982).  The  mean 

value  of  Y  from  the  simulations  was  1.61.  Most  of  the  energy  soread  comes 
z 

at  the  transition  between  the  induction  li.nac  transport  system  and  the 
solenoidal  field.  Figure  18  is  a  plot  of  the  induction  linac  FSL  magnetic 
field  profile.  The  diode  is  in  a  field  free  region  at  z  *  0  and  the  last 
focusing  coil  of  the  linac  is  at  350  cm.  The  edge  of  the  solenoid  is  at 
approxima tely  375  cm.  The  2  kG  magnetic  field  i3  required  for  a  beam 
equilibrium  radius  that  is  consistent  with  a  practical  wiggler  iiameter. 
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Figure  19  is  a  plot:  of  results  from  a  static  computer  simulation  of  the  beam 


transport  from  the  cathode  into  the  solenoid.  The  last  accelerating  gap  of 
the  induction  module  is  not  energized,  so  y  =  2.1  as  the  beam  enters  the 
last  focusing  coil.  When  the  beam  has  reached  the  uniform  field  region  <v  > 
is  1.6. 

This  problem  of  beam  transport  and  matching  represents  one  of  the  ma^or 

design  concerns  for  high  current  free  electron  lasers.  The  equilibrium  of 

high  current  beams  is  sensitive  to  self  electric  and  magnetic  fields  of  the 

beam.  Abrupt  changes  in  the  wall  diameter  or  magnetic  field  profile  induces 

betatron,  or  zero  frequency  cvclotron  oscillations  at  u  =  0  =  !<3c  -  ~  /y ,  so 

o 

the  oenod  is  X  =  2i%cy/1  . 

o 

A  spread  in  3  across  the  beam  radius  will  cause  these  oscillations  to  phase 
mix  after  a  few  periods.  The  perturbation  then  appears  as  an  effective 
temperature. 

In  the  present  simulation  the  initial  normalized  emittance  was 

s  =  180  t  mrad-cm,  whereas  the  effective  final  emittance  was 
n 

z  =  yr3,  =  940  t  mrad-cm.  Almost  ail  of  this  increase  in  emittance 

n  1  rms 

came  in  the  transition  to  the  solenoid.  This  is  why  the  final  emittance  for 
the  hot  and  cold  cathode  is  about  the  same.  These  problems  can  and  have  been 
avoided  in  Marx-pulse  line  beam  generators  by  immersing  the  diode  in  a  very 
strong  magnetic  field.  This  minimizes  the  radial  excursion,  and  effective 
emittance  growth.  This  approach  becomes  difficult  when  long  pulse  times  are 
required  and  one  must  use  large  area  hot  cathodes  to  get  kiloampere  electron 
beams.  One  can  conclude  from  this  that  the  design  of  the  beam  transport 
system  for  an  induction  linac  :E1  is  at  least  as  important  as  the  diode 
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Figure  19.  Variation  of  y,  with  z  for  electrons  in  a  computer  simulation 

of  the  experiment.  The  large  spread  in  y,  occurs  as  the  beam  enters 


the  solenoidal  field. 


1 

VIII.  Conclusions  and  Future  Directions 

Improvements  in  the  beaut  transport  are  clearly  needed  to  operate  in  the 
Raman  regime.  Cne  of  the  advantages  of  operating  in  the  cold  beam  limit  is 
that  the  efficiency  of  the  FEL  is  higher  and  can  be  enhanced  by  tapering  the 
wiggler  period.  Once  the  beam  is  trapped  in  the  pondermotive  wave,  the  phase 
velocity  of  the  wave  can  be  reduced  gradually  to  extract  energy  from  the 
beam.-  Figure  20  is  a  plot  of  efficiency  vs  distance  for  an  untapered  and 
tapered  wiggler.  These  results  are  from  a  computer  simulation  of  the  present 
experiment.  Whereas  the  intrinsic  efficiency  is  about  5  percent,  the  tapered 
wiggler  efficiency  is  in  excess  of  20  percent. 


To  go  to  shorter  wavelengths  the  FEL  can  be  operated  as  a  two  stage 

2 

device.  The  output  from  stage  one  at  X,=  \yiy  can  be  used  as  the  output 
for  the  second  stage  to  give 


.  2  a  4 

4Y  8y 


(51) 


To  obtain  a  high  voltage  high  current  accelerator  and  achieve  a  high  gain, 
high  power  short  wavelength  FEL,  the  long  pulse  induction  module  can  be 
converted  into  a  racetrack  accelerator  (Roberson,  1981’,  Mondelli  and  Roberson, 
1932).  This  takes  advantage  of  the  long  pulse  to  give  a  voltage 


V  =  V  T/T 
9 


where  V  is  the  voltage  aain  of  the  module,  T  the  time  the  module  is  on,  and 

g 

T  the  time  it  takes  the  beam  to  go  around  the  racetrack.  In  the  present 
case  a  30  nsec  path  length  would  result  in  a  26  MeV  beam  with  the  2  js  ,0.4 


MV  module. 


Figure  20.  Theoretical  FEL  efficiency  vs.  interaction  length  with  an  without 
a  tapered  wiggler.  The  parameters  used  for  the  calculation  are 
those  of  the  experiment. 
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VF  and  TF  Coil  Design  Objectives 


1)  Maintain  Field  Uniformity  and  Low  Field  Errors 


2)  Minimize  Energy  Losses  (<10%) 


3)  Minimize  Deflections  (  <  1mm) 


4)  Establish  a  Reliable  and  Accessible  System 
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List  of  Svbols 


rv  :  RADIUS  OF  VERTICAL  BUS  CONDUCTOR 

:  LENGTH  OF  VERTICAL  BUS  CONDUCTOR 

V 

Z  :  AXIAL  POSITION  OF  THE  BOTTOM  OF  THE  EUS 

B3us  :  IS  THE  VF  (AVERAGE) IN  THE  REGION  WHERE  THE  BUS 

IS  GOING  TO  BE  LOCATED 

b:  SIDE  LENGTH  OF  SQUARE  CROSS  SECTION  COIL 

<3  :  ELECTRICAL  CONDUCTIVITY 

T  !  TIME  TO  REACH  PEAK  CURRENT 

S  :  TOTAL  CONDUCTOR  LENGTH 

L  :  COIL  OR  CONDUCTOR  INDUCTANCE 

^LOSS^COIL  1  FRACTIONAL  ENERGY  LOSS  (MAGNETIC  AND 
RESISTIVE  UNLESS  SPECIFIED  OTHERWISE) 

F  :  FORCE  PER  UNIT  LENGTH 
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List  of  Symbols,  Con't. 


d :  LENGTH  OF  COIL  BETWEEN  FIXED  POINTS  (SPAN) 

a  :  DEFLECTION  DISTANCE 

Y:  Young's  modulus  of  elasticity 

Tdf:  diffusion  time  through  a  thin  wall  cylinder 

t  :  WALL  thickness  of  coaxial  bus, 
y0  :  permeability  of  free  space 

Rt  :  total  AC  coil,  bus  and  joint  resistance 

* 

Lt  :  total  coil  and*  feus  inductance 
LBlie  :  bus  inductance 

ouS 

Leo i l  :  C0IL  INDUCTANCE 

Eco I L  :  MAGNETIC  field  energy  from  coils  only 

C  :  CAPACITANCE 
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Consider  a  sguare  cross  section  coil  with  side  of  length,^ 


9  Clearance  requires:  b  <  7cm 
o  Low  energy  loss  requires  the  conductivity,  a  > 

o  Shall  deflection  (a)  requires  the  modulus,  y  >  -E_  /AW 

~  32a  \  t)  / 

For  T  =  3,2  ms,  S  =  170  m,  L  =  590  JJH,  EL0SS/EC01L  =  .1 

7 

(resistive  only) ,  d  =  70cm,  a  =  ,1cm,  F  =  10'  dYNEs/cm 
(Hoop  Stress)  and  b  =  5cm  then 

~  4  x  10^  (Q-m  )~1  ,  Y  ~  loll  dYNEs/cm  2 

=> Aluminum  or  Copper  can  satisfy  these  Requirements. 

Aluminum  is  selected  when  cost  and  weight  are  considered. 

-so¬ 


ts 


'COIL 


(_b2  l^LOSS 


Consider  a  sguare  cross  section  coil  with  sides  of 


Field  Errors  from  Eus  Eddy  Currents  Induced  by  VF 
For  no  diffusion  through  a  vertical  bus,  the  VF  error  $ 


IS  GIVEN  UNDER  THE  BUS  ON  THE  MIQPLANE  BY 


For  rv=  5cm  ,  &v=  100cm  and  located  at  r=  90cm  , 


z  =  33cm  then  B  ^us  =  1.7  Bz  and  at  r=  90  cm,  z  =  0 


0  ABZ/BZ  ~1  %  (Worst  Case,  inner  VF  bus) 

=>ThE  INNER  VF  BUS  MUST  BE  MADE  MAGNETICALLY  TRANSPARENT 


and/or  HAVE  A  SMALL  MAGNETIC  MOMENT, 

-=>ThE  TF  BUS  CAN  BE  LOCATED  IN  A  REGION  OF  LOW  VF. 


==C>BUS  CONDUCTORS  ABOVE  AND  BELOW  THE  DECK  GIVE  SMALL 
(  < .1%)  VF  ERRORS. 
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VF  and  TF  Bus  Selection 


Consider  coaxial  bus  with  outer  radius  a 

VF  Bus 

o  Coil  compatible  and  clearance  require:  a<5  cm 

2tdf 

o  Magnetic  diffusion  requires:ct< — ^ 

At/6 

TF  Bus 

o  Coil  compatible  and  clearance  requires:  a  <  10  cm 

•  Low  loss  requires:  cr  >  — 

L27TAt 

VF,  rDF=  ,5  ms,  a  =  2,5  cM,t  =  1  cm 
cr <  3  x  IQ6  Ca-M)-1 

Stainless  Steel  or  Equiv.  is  required. 

IE,r=3.6  ms,  s  =  24  m,  L  =  3yu.H,A  =  1’cm,  a  =  5  cm 

q-^9  x  106  (fl-M)'1 

Aluminum  or  Copper  will  work. 
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E 

COIL 

^LOSS, 


T.  F.  COAXIAL  BUS  TO  COIL 
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Joint  Problem 


•  Multilam-a  louvered  interface  that  provides  many 

POINTS  OF  CONTACT 


e  Joints  are  demountable  and  don't  require  plastic 

DEFORMATION  OF  JOINT  CONDUCTORS 

•  LOW  RESISTANCE  (~  4]jQ/jOINT  WITH  210  LOUVERS  FOR  TF) 
AND  HIGH  CURRENT  0210  «A  FOR  TF,  1  KA/LOUVER 
TESTED  >1  SEC)  JOINTS  ARE  POSSIBLE 
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COIL  TO  COIL  V.  F.  BUS 


jRC.-AL  .  .EL 


ERROR 


CAT  I NG  BUS  WILL  REDUCE  ERROR/  BUT  FORCE  IS  STILL  DIFFICULT  TO  HOLD 
K  LBS/FT)/  COAXIAL  BUS  APPEARS  TO  BE  THE  BEST  SOLUTION. 


TOTAL  ENERGY  AND  CIRCUIT  PARAMETERS 


After  conductor  selection  the  total  energy  can 

BE  CALCULATED, 

0  ^T  ^COIL 

From  Table  V, 

ETF  =2.12  MJ  and  Evp  =  .66  MJ 

The  capacitance  is  determined  from 

rs  4LT/RT2 

•  x  =  tan  ( - 1  x),  WHERE  C  =  1  ~ 

2LT  x2  +  1 

See  table  V. 


r  R_  Lallc 
1  +  T  +  J1AS. 

T  COIL 
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I  Al»7  H7  THE  MODIFIED  BETATRON  ACCELERATOR  EXPERIMENT < U )  dAYCOR 
ALEXANDRIA  VA  F  M  MAKO  06  JAN  84 
dAYCOR -0206-84-001/6207  NOOO 1 4 - 8 1 -C-2095 


NCLASSI F I  ED 


F/G  20/7 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURIAU  01  STANDARDS  196i- A 


PARAMETERS  AND  SPECIFICATIONS 

TF  SYSTEM 

VF  SYSTEM 

COIL  AC  RESISTANCE 

0.65 

2.9 

BUS  AC  RESISTANCE  (M$ 

0.52 

12.4  . 

JOINT  DC  RESISTANCE  (l^ 

0.4 

2 

COIL  INDUCTANCE  (fA\) 

85 

590 

BUS  INDUCTANCE  (f*M) 

4.1 

8.7 

TIME  TO  REACH  PEAK  CURRENT  (MS) 

3.6 

3.2 

PEAK  CURRENT  (kA) 

210 

45 

BANK  ENERGY  (MJ)  (NOMINAL) 

2.12 

0.66 

BANK  CAPACITANCE  (mF) 

60.8 

7.2 

BANK  VOLTAGE  (kV) 

8.5 

13.6  • 

1 

COIL  CROSS  SECTION  AND 

MATERIAL 

4"x4'W  WALL 
6061-T6  alum. 

3"x3"x*l"  WALL 
6061-T6  alum. 

VERTICAL  BUS  TO  COIL,  TYPE, 

SIZE  AND  MATERIAL 

CIRCULAR  COAX., 

4"  0,D.,3.36"I.D. 

OUTER  CONDUCTOR 
1,5"  DIAM.  INNER  ‘ 
cond.,  6061 -T6 

ALUMINUM 

CIRCULAR  COAX., 

3.5"  0.D.,3"  I.D. 

OUTER  CONDUCTOR 

304  S.S.,  1"  DIAM., 
INNER  COND.,  6061 -T 
ALUMINUM 

COIL  TO  COIL  BUS,  TOP  AND/OR 
BOTTOM  DECK,  TYPE,  SIZE  AND 

MATERIAL 

RECTANGULAR  COAX. 
4"x8"x*i"  WALL 
OUTER  COND.,  l"x 
5"  INNER  COND. 
6061-T6  ALUMINUM 

SAME  AS  TF  BUT  AT  V 
AND  BOTTOM  OF  DECKS 
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J 

Return  Strap 


Joint 

of  I  F  Coil  To  Coil  Connection 


Summary  and  Conclusions 


Coil  and  busing  design  issues  are  resolved 


Energy  losses  have  been  minimized 


Low  field  errors  can  be  attained 


Deflections  can  be  kept  small 


Detailed  bus  work  design  is  in  progress 


External  Injection  into  a  High  Current 
Modified  Betatron  Accelerator 

F.  Mako,  W.  Manhelmer,  C.  A.  Kapetanakos 
Naval  Research  Laboratory,  Washington,  DC  20375 
0.  Chernln 

Berkeley  Research  Associates,  Springfield,  VA  22150 
and  F.  Sandel 

Sachs  Freeman  Associates,  Bowie,  MO  20715 

NRL  Memorandum  Report  5196 
September  9,  1983 


Abstract 

By  axially  tapering  the  current  density  at  each  end  of  an  axial  pinch  a 
scheme  is  developed  for  injecting  electrons  across  magnetic  field  lines  into  a 
modified  betatron  accelerator.  This  scheme  produces  only  a  minimal 
perturbation  to  the  circulating  electron  beam.  The  axial  pinch  provides  beam 
equilibrium  in  the  transverse  field  and  the  tapering  reduces  the  non-adiabatic 
growth  of  the  perpendicular  particle  velocity  from  0.5  c  to  .05  c  for  the 
parameters  of  the  NRL  modified  betatron. 


Introduction 


Presently,  there  is  substantial  interest  In  increasing  the  current  limits 

of  the  conventional  betatron.  4  Several  of  these  new  approaches  require  at 

3-6 

least  a  large  toroidal  magnetic  field  in  addition  to  the  conventional 
betatron  field.  The  added  toroidal  field,  however,  makes  the  injection  into 
the  accelerator  considerably  more  involved,  because  the  injected  beam  must  be 
transported  across  field  lines  before  entering  the  torus. 

In  the  modified  betatron,  i.e.,  an  accelerator  that  combines  a  betatron 

and  torodial  field,  this  difficulty  was  avoided  by  locating  the  injector 

inside  the  torus. ®  Such  an  injection  scheme  has  several  advantages  but  also 
three  short  comings.  First,  it  requires  a  large  opening  on  the  torus,  which 
introduces  large  field  perturbations.  Second,  the  debris  of  the  diode  can 
have  an  adverse  effect  on  the  quality  of  the  vacuum  system  and  third,  it 
requires  an  injector  with  a  short  fall  time  to  avoid  perturbing  the  beam  after 
the  first  revolution. 

To  avoid  the  difficulties  of  internal  injection,  we  have  developed  a 
scheme  to  propagate  the  beam  across  magnetic  field  lines  and  thus  to  locate 

the  injector  outside  the  torus.  In  the  proposed  scheme,  an  axial  pinch 

provides  an  equilibrium  for  the  beam  to  cross  field  lines.  Tapering  of  the 
current  density  at  both  ends  of  the  axial-pinch  is  necessary  to  minimize  the 
increase  of  the  perpendicular  particle  velocity  which  results  from  the  axial- 
pinch. 

I.  Electron  Beam  Motion 

The  proposed  external  Injection  scheme  for  the  high  current  modified 
betatron  is  shown  in  Fig.  1.  The  anode  of  the  injector  is  located  at  y  ■  o 

i 

and  the  intense  beam  propagates  in  a  combined  transverse  and  axial  external 
magnetic  field  until  it  reaches  the  torus.  The  components  of  the  magnetic 
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field  for  a  uniform  currenC  density  channel  are  given  byt 


B  (z)  -  B*  *  B 


a 


xo' 


B,  (x) 


-B  i 
%  a* 


yo 


(O 

(2) 

(3) 


where  is  the  magnitude  of  the  plasma  current  magnetic  field  at  the  plasma 


radius  a  and  Bxo,  ByQ  are  the  cartesian  components  of  the  toroidal  magnetic 

\ 

field  (Bn)  at  the  beam  trajectory.  The  fields  (B  ,  B  ,  B  )  are  taken  to  be 
□  *  y  xo  yo 


constants  and  the  vertical  magnetic  field  (Bzq)  is  ignored  since  its  effect  is 


the  same  as  that  of  In  addition,  the  self  fields  are  ignored,  because 


the  axial  pinch  plasma  density  is  5-6  orders  of  magnitude  greater  than  a 
typical  intense  electron  beam  (~  IKA/cm^)  and  thus  charge  and  current 
neutralization  is  assumed. 


After  a  change  of  variables  (from  time  t  to  axial  position  y)  using  the 


transformation 


dt 


°y  dy  an<*  the  ParaxiaI  ray  approximation 


2  2 

[  ( x "  +  z'  «  1) ,  where  the  prime  is  a  derivative  with  respect  to  yj  for  the 


transverse  motion,  the  equations  of  motion  become 

ft 


X  *  -z 


<5^ 


<*) 


a^.  n  q 


u  a  "  ft  '  u 

o  i>  o 


(5) 


dt 


where  ft. 


MZe,  a 

ay  y 


|e|  B 


4 


1  +  x'2  +  z'2 


(6) 


,  n„ 


e  B 
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,  and  e  ,  m,  y  and  u  are  the  i 

o  +  , 


i  vt  " 


♦  my  ’  y  my  '  x  my  '  i  «'  ■  o 

electronic  charge,  mass,  relativistic  factor  and  the  magnitude  of  the (velocity 

vector  respectively*  Note  also,  the  paraxial  ray  approximation  implies; 


flx  «  ft^  <  fty  and  small  displacements. 
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According  to  Eq.  (5)  the  B  field  displaces  the  particle  orbits  along  the 

ax 

z-axis  by,  -  -g-  a.  This  is  a  general  result  and  not  a  consequence  of  the 
♦ 

paraxial  ray  approximation.  Beam  containment  requires  therefore  that, 
ft  «  ft  that  is,  the  peak  field  of  the  plasma  current  must  be  large 

X  y 

compared  to  the  crossed  magnetic  field. 

The  solution  of  Eqs.  (4)  and  (5)  is  given  by 


~a  j-  nx  1 

z  (y)  -  -$T  a  +  xo  sin  kvy  +  (zQ  +  ft"  a)  C0S  kyy  COS  ky 

ib  l-  *  ib  J  *s 


to 


*  y  1  gin  IcV 

-  xoky)  cos  kyy  +  <xo  +  2oky  +  hy  Q-  *)  sin  kyyj  — jp-2,  (7) 


x  (y)  =>  £xq  cos  k^y  -  (zq  +  ■—  a)  sin  cos  ky 

+  ['  (zo  “  xo  ky}  sln  kyy  +  (ky  zo  +  xo  +  ky  ft^  a)  cos  kyy]  “ 1 T^’  (8) 

and  indicate  oscillatory  motion  around  the  displaced  axis.  In  Eqs.  (7)  and 
/  2  *  * 

( 8 )  k  »  — k  =‘\k  +  — and  x  ,  z  ,  x  ,  and  z  are  the  initial  positions 

2u  y  au  o  o  o  o 

o  o  ■ 

and  velocities  respectively.  Particle  motion  is  characterized  by  two 
frequencies,  a  fast  cyclotron  motion  at  Qy  and  a  slow  rotation  about  the  y- 
axis  with  a  frequency  of  ft-  -  u  ft  /a  ft  ,  where  u  is  the  magnitude  of  the 
initial  axial  velocity. 

Although  the  axial  pinch  can  provide  the  means  for  crossed  field  transport 
it  also  is  the  source  for  the  growth  of  perpendicular  energy.  From  Eqs.  (7)- 

4 

(8),  with  xQ 


zq  "  o,  the  perpendicular  energy  Is  found  to  be. 


eA  (y)  -  2ra  m 


yo 


“;v 


[iTil 


Since  »  Qx  most  of  Che  perpendicular  energy  comes  from  the  pinch, 
except  for  particles  near  Che  axis*  However,  in  practice,  most  electron  beams 
have  more  particles  at  a  large  radius. 

An  accurate  calculation  of  the  perpendicular  velocity  cannot  be  based  on 
the  paraxial  equations.  However,  the  more  general  equations  are  not  easily 
tractable  analytically.  For  this  reason,  it  is  assumed  that  *  o.  Under 
this  assumption,  it  can  be  shown  from  the  conservation  of  axial  Py  and 


canonical  angular  momentum  that, 

dg 


%  ^2a 

and  — — -  is  the  small  expansion  parameter. 

/Y® 

From  (9)  the  maximum  radial  velocity  is. 


(9) 


S  9hJ  t-o 

where  pQ  is  the  initial  value  of pi  For,  By  ■  2KG,  *  2KG,  pQ  =»  a  -  2  cm, 
and  Y  ■  7,  the  ratio  _£max  , 

which  is  too  large  to  be  acceptable. 

uyo 

By  numerically  integrating  the  exact  relativistic  particle  equations  we 
obtain  Upmax/c  *  0.45  (where  c  *  speed  of  light).  The  nu  ^rically  calculated 
radial  velocity  is  shown  in  Fig.  2a.  The  results  of  Fig.  2a  were  obtained  for 
the  general  case,  i.e.,  when  a  transverse  field  is  present,  in  addition  to  the 
axial  and  plasma  current  magnetic  fields.  The  values^of  the  various 
parameters  are  listed  in  the  figure  captions.  Shown  also  in  Fig.  2a  are  the 
axial  profiles  of  the  fields. 
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A  subtantlal  reduction  of  the  perpendicular  velocity  can  be  achieved  Sy 
axially  tapering  the  current  density  at  both  ends  of  the  pinch.  Although  « 
pinch  field  is  the  primary  source  of  u^,  non-adlabatic  effects  contribute 
significantly  to  u^  but  can  be  removed  by  an  adiabatic  taper.  Results  are 
shown  in  Fig.  2b.  For  a  large  taper  length  a  perturbation  calculation 
predicts  a  u^  that  is  within  20%  of  the  numerical  results.  To  substantially 
reduce  the  perpendicular  velocity  the  taper  length  must  be  several  times  the 
betatron  wavelength.  An  expression  for  the  betatron  wavelength  can  be 
calculated  from  Eqs.  (6),  (7),  (8)  and  is  given  by, 

X  -  2*  /£  [i  +  4  ?°rl/2 
y  y 

Figure  3  shows  the  reduction  of  the  perpendicular  velocity  as  a  function  of 
the  ratio  of  taper  length  to  betatron  wavelength,  this  is  shown  with  and 
without  transverse  fields.  For  By  »  2KG,  *  2KG,  a  =  2  cm,  y  ■  7  the 
betatron  wavelength  is  approximately,  X  *  10.4  cm.  To  reduce  u^/c 
substantially  would  require  a  taper  length  of  20-30  cm. 

The  feasibility  of  tapering  the  current  density  has  been  tested 
experimentally.  Presently,  a  taper  length  of  25  cm  has  been  achieved 
experimentally.  The  measured  azimuthal  magnetic  field  B^  of  the  pinch  with 
and  without  tapering  is  shown  in  Fig.  4.  This  current  density  profile  was 
achieved  by  geometrically  tapering  the  vacuum  chamber.  A  cylindrical  resistor 
taper  is  presently  being  examined  for  making  a  radially  compact  taper. 

II  Plasma  Interaction  - 

Calculations  indicate  that  the  axial  pinch  does  not  substantially  interact 
with  the  electron  beam  nor  with  the  external  fields.  However,  there  are 

% 

limits  imposed  on  the  plasma  density.  To  avoid  an  unacceptably  large 
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displacement  of  the  plasma  column  from  the  J  x  B  force  and  enhancement  of 

y  * 

the  beam  emittance,  the  plasma  density  should  satisfy  the  Inequalities  (in 


2  fl2  B,  B  t 

_ £  B _ >  n  >  *  xo 

8r  DZ  (Z  +  1)  lx1  p2ln  [a  Y&2/2Z4/3r  ]  ^  AzAmpa 

e  o  o  e 


where  e  is  the  maximum  beam  emittance,  l  is  the  length  of  the  plasma  column, 

Z,  A  are  the  atomic  number  and  weight  of  the  gas,  Az  is  the  maximum  plasma 

e2  X  2 

channel  displacement,  mp  is  the  proton  mass,  rg  =  — j,  aQ  = **— j,  D  =  1  for 
0  me  me 

monotoraic*3  2  for  diatomic  gas.  Also,  T  is  the  time  that  has  elapsed  from  the 

beginning  of  the  pinch  to  when  the  electron  beam  is  injected.  Instabilities 

will  be  considered  elsewhere. 

Summary 

Using  a  tapered  axial  pinch  of  suitable  plasma  density  and  current  appears 
so  far  to  provide  a  low  perturbation  scheme  for  beam  injection  Into  a  modified 
betatron.  Extensive  analytical  and  computational  work  indicates  that  the 
transverse  particle  velocity  at  the  exit  of  the  plasma  column  is  substantially 
smaller  than  in  an  untapered  pinch.  Preliminary  experimental  results 
indicates  that  a  pinch  with  tailored  axial  current  density  is  achievable. 
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X 


CURRENT  DENSITY  PROFILE 

Figure  1.  Schematic  of  tapered  axial-pinch  injector. 
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B  (KG) 


FIELDS  VS.  y 


y  (cm ) 


Transverse  particle  velocity  vs  axial  positions.  The  particle 

parameters  are:  y  ■  7,  y  •  0  and  P.  ■  a  »  2  cm.  Taper  length 
Y  ±  v 

from  10-90  percent  amplitude  is  2.2  cm.  C  »  speed  of  light. 


Figure  3.  Maximum  transverse  particle  velocity  vs  taper  length.  The  same  parameters  as  in  Figure  2a  but 
with  variable  taper  length  and  crossed  fields. 
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Electron  trajectories  through  a  linear  magnetic  wiggler  in  an  axial  guide  field  are  calculated 
numerically.  Off-axis  electrons  are  shown  to  drift  in  a  direction  perpendicular  to  the  wiggler  field 
because  of  the  wiggler  gradient.  Effects  of  self-fields  and  initial  conditions  are  analyzed,  and  the 
results  compared  to  those  obtained  with  a  helical  wiggler.  An  empirical  analytic  expression  for  the 
linear  wiggler  drift  is  presented,  and  means  of  stabilizing  the  drift  are  discussed. 

PACS  numbers:  41.70.  + 1,  41.80.Dd 


Recently,  there  has  been  much  interest  in  free  electron 
lasers  (FEL'si  which  use  spatially  varying  magnetic  fields  to 
modulate  a  relativistic  electron  beam  ;R£B).'"!  Three  types 
of  magnetic  wigglers  have  been  widely  used:  helical  (usually 
produced  by  a  bifilar  helical  current  winding4-7),  linear  ipro- 
duced  by  an  array  of  permanent  magnets*  or  by  linearly  al¬ 
ternating  windings’),  and  radial'0  I  produced  by  a  series  of 
spaced  conducting  or  ferromagnetic  washers"  immersed  in 
an  axial  field  or  by  a  series  of  alternating  coils3).  Of  these, 
only  the  first  two  produce  a  perpendicular  field  on  axis  and 
are  therefore  suitable  for  use  with  a  small  diameter  beam  of 
solid  cross  section.  In  this  paper,  we  will  analyze  the  proaga- 
tion  of  a  solid,  high  current  REB  through  linear  and  helical 
wigglers  with  a  superimposed  axial  guide  field. 

FEL  experiments  fall  into  two  categories  depending  on 
the  beam  current  Compton  regime  FEL’sl-JJ  have  used 
high  energy  I  tens  of  MeV),  low  current  I  5  1  A)  beams  while 
Raman  FEL’s3-4  have  used  lower  energy  I  —  1  MeV),  high 
current  ( 2  1  kA)  beams.  Helical  wigglers  have  been  used  in 
both  current  regimes,  but  until  recently9  linear  wigglers  were 
used  only  in  the  Compton  regime.  In  this  case,  the  linear 
wigglers  consist  of  permanent  magnets  and  there  is  no  guide 
field.  However,  for  high  current  beams,  a  guide  field  is  re¬ 
quired  to  contain  the  beam.  Although  there  are  advantages 
to  a  linear  wiggler  from  the  standpoint  of  ease  of  assembly 
and  versatility  le.g.,  changing  the  periodicity  or  tapering  the 
period  and/or  field  amplitude  for  efficiency  enhancement12) 
we  will  show  that  there  is  no  equilibrium  for  off-axis  parti¬ 
cles  when  a  beam  is  propagated  through  a  linear  wiggler  in 
an  axial  field  The  particles  drift  out  of  the  wiggler  unless 
additional  focusing  forces  are  provided 

Various  authors  have  considered  particular  cases  of 
electron  motion  through  magnetic  wigglers.  These  have 
typically  been  single  particle  (low  current)  calculations. 
Blewett  and  Chasman7  considered  motion  of  high  energy 
electrons  (—24  MeV)  through  a  helical  wiggler  and  found 
stable  helical  orbits  with  superimposed  betatron  oscillations. 
Friedland11  has  treated  the  case  of  lower  energy  electrons 
(—300  keV)  in  an  idealized  radially  uniform  helical  wiggler 
with  a  superimposed  axial  guide  field  He  showed  that  var¬ 
ious  stable  trajectories  were  possible  and  derived  stability 
criteria  relating  the  allowed  wiggler  and  guide  field  strengths 
to  beam  energy  and  wiggler  period  .  These  “stable”  re¬ 
gions  are  given  by 

“JAYCOR,  AJtuafr*.  VA  21304 


and 


(2a) 


iiQ>ck, 


(2b) 


where  30  is  the  axial  guide  field  J2,  *  eBg/ym  is  the  corre¬ 
sponding  cyclotron  frequency,  <>  is  the  electron  velocity,  and 
k  =»  2tr/Aw .  Physically,  these  conditions  stem  from  the  reso¬ 
nance  in  the  perpendicular  velocity  of  an  electron  in  the 
wiggler  and  guide  fields4: 


»,J20 
o,k  —  J20 


(3) 


Freund  and  Drobot14  have  considered  this  case  further  and 
also  conclude  that  stable  trajectories  with  nearly  constant 
axial  velocities  and  relatively  large  wiggler  amplitudes  are 
possible  when  fi0<ck.  This  is  consistent  with  condition  |2a). 

The  present  analysis  employs  a  relatively  simple  com¬ 
puter  code  which  solves  the  equations  of  motion  of  an  elec¬ 
tron  in  any  electric  and  magnetic  field  configuration  using  a 
fourth-order  Runge-Kutta  method.  Self-electric  and  mag¬ 
netic  fields  are  calculated  by  assuming  that  the  electron  is  at 
the  edge  of  an  azimuthaUy  symmetric  beam  of  current  I,  so 
that  the  self-fields  can  be  written  as  imks) 


E. 


-/ 

2jre0/v, ' 


B , 


2irr 


(4) 


where  tit )  is  the  electron  radius.  The  axial  self-fields  are  ne¬ 
glected. 

For  the  particular  cases  considered  here,  the  external 
magnetic  field  consists  of  a  soienoidal  field  produced  by  a 
15.3-cm-i.d.  x  2-m-iong  solenoid  together  with  a  wiggler 
that  begins  in  the  uniform  portion  of  the  soienoidal  field.  The 
wiggler  field,  which  may  be  either  helical  or  linear,  rises 
adiabaticaily  over  ten  periods  and  then  oscillates  with  con¬ 
stant  amplitude.  The  envelope  enclosing  the  wiggler  ampli¬ 
tude  is  given  by 


0<z<lCW. 

*>10*,.. 

(5) 
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Here,  z  is  the  axial  distance  from  the  beginning  of  the 
wiggler.  This  variatioa  closely  fits  the  linear  wiggler  used  in 
the  NRL  induction  linac  FEL  The  field  on  axis  from  this 
wiggler  is  plotted  in  Fig.  1  along  with  the  envelope  equation 

(5). 

The  linear  wiggler  field  components  are 
B‘,=b  [z\Bw  cosh  far  cos  kz, 

B‘,-0,  16) 

Si  —  —  b  (z)S„  sinh  far  sin  fcr, 

where  S.  is  the  peak  wiggler  field  on  axis.  For  the  NRL 
linear  wiggler.4  which  has  A w  —  3  cm  and  winding  layers 
spaced  by  3.2  cm,  these  expressions  are  valid  for  r  S  1  cm. 

We  will  be  primarily  concerned  with  relatively  low  en¬ 
ergy  (~  1  MeV),  high  current  1—kA)  beams.  Of  particular 
interest  are  the  effect  of  self-fields,  guide  field  and  wiggler 
amplitude,  and  initial  beam  conditions  corresponding  to  a 
field  immersed  {v,0  =  0)  or  a  shielded  source  \Pt0  —  01, 
where  u,0  is  the  initial  azimuthal  electron  velocity  and  P,0 
the  initial  canonical  angular  momentum.  If  the  axial  self- 
magnetic  field  is  neglected,  P,a  —  0  implies  ut0  —  erB,/ 
(2 ym)  in  a  uniform  field  B,. 

First,  we  consider  a  beam  with  y  -  2.2,  I  —  230  A, 
B0  —  2  kG,  Bm  —  1  kG,  and  *  3  cm.  Note  that  Eq.  (2a) 
would  require  S„  <  3  kG  for  stability  in  this  case  with  a 
helical  wriggler.  Thex-y  trajectories  of  electrons  injected  with 
i\o  —  0  and  r0  —  0.4  cm  into  a  linear  wiggler  are  shown  in 
Fig.  2.  Each  particle  initially  begins  to  E  x  B  drift  in  the  self- 
fields,  but  when  it  reaches  a  region  of  large  wiggler  field  it 
begins  drifting  in  they  direction.  In  each  case  the  electron 
reaches  an  imaginary  wall,  located  at  r  —  1.2  cm  only  par¬ 
tially  transversing  the  wriggler.  It  should  be  noted  that  the 
assumptions  used  to  calculate  the  self-fields  become  invalid 
as  the  beam  distorts.  Consequently,  to  determine  if  these 
self-fields  are  responsible  for  this  drift,  the  calculations  are 
repeated  with  /  —  0  (Fig.  3|.  A  third  trajectory  is  also  plotted 
for  Bw  —  0.3  kG-  The  electron  that  originates  on  they  axis  is 
well  confined,  but  electrons  off  the  y  axis  again  drift  to  the 
wall,  with  a  velocity  much  higher  for  Bm  —  1  kG  ((»*> 
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FlO.  1.  Plot  of  niinrwH  t'/X  from  NTU.  linaar  wiggkr  lOfatber  with 
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FIG.  2  Electron  tnyocton—ia  linear  wiggler  with  a  »  -  j  cm.  B,  -  2  kG. 
r  -  12  r,  -  0  4  cm.  /  -  230  A.  and  B.  -  1  kG. 

=0.047  c)  than  for  Bm  -  0.5  kG  |  (vg  >=0.01 1  c). 

This  behavior  can  be  explained  by  the  increasing  gradi¬ 
ent  in  the  wiggler  field,  and  consequently  the  emergence  of  a 
significant  axial  field  component,  as  the  distance  from  the  y 
axis  increases.  Although  the  present  drift  arises  from  the 
gradient  in  B„  and  is  in  the  direction  ofVBw  xBo.  ltis  quan¬ 
titatively  different  from  the  usual  guiding  center  approxima¬ 
tion  because  the  field  variation  over  one  gyroperiod  is  so 
extreme.  For  example,  when  far  —  0.8  ( x=0.4  cm  in  the 
present  case),  B',  varies  from  ■+•  Bm  to  —  S„  over  one  peri¬ 
od.  However,  the  physical  mechanism  is  the  same  as  for  the 
usual  gradient  drift;  i.e.,  the  gyroradius  in  the  part  of  the 


FIO.  3.  Electron  tnjectones  m  linear  wiggler  with  .4.  -  3  cm.  t,  «  2kG, 

r  •  2.2  m  0.4  as.  and  /  -  0.  B.  w  1  kG  tot  parades  s  sad  b  seed  Bm 

m  0.3  kG  Tor  parade  C. 
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TABU  I.  Companion  of  drift  velocities  computed  from  the  code  with  those  from  Eq.  l9l  for  vinous  cases 


JolkGl 

i„,k Gl 

V 

Tlcmi 

Vt/e  |co0*i 

9"«/e(Eq. 

2 

1 

2.2 

0.29 

0.4 

0.047 

0.030 

X 

1 

12 

0.37 

0.4 

0.042 

0.04« 

2 

0.5 

12 

029 

0.4 

0.011 

0.013 

-> 

l 

12 

0.29 

0.2 

0.019 

0.022 

Z 

i 

3.0 

0.21 

0.4 

0.029 

0,034 

X 

1 

3.0 

0.42 

0.4 

0.019 

0.022 

X 

0.3 

3.0 

0.42 

0.4 

0.0031 

0.0039 

X 

l 

10.0 

0.11 

0.4 

0.0043 

0.0049 

10 

J 

10.0 

0.2S 

0.4 

0.037 

0.039 

1  orbit  where  B,  is  a  minimum,  or  where  ||  xi|  is  a  minimum  in 

(  Fig.  3,  is  larger  than  where  3.  is  a  maximum.  The  addition  of 

'  seif-fields  merely  imposes  an  additional  ExB  rotation  on 

i  the  orbit,  so  that  an  electron  that  originates  on  the y  axis  i  and 

is  therefore  confined  when  /  =  0)  begins  to  drift  into  a  region 
of  increasing  3 ' .  Therefore,  it  actually  is  “lost”  sooner  than 
an  electron  on  the  x  axis  which  initially  Ex  B  drifts  into  a 
region  of  smaller  3 ' . 

An  approximate,  empirical  expression  for  the  drift  has 
been  found  which  is  in  quite  good  agreement  with  the  code 
results  when  <  f30  <  leu,  (below  the  cyclotron  reso¬ 
nance  I.  In  terms  of  the  wiggier  gradient,  the  expression  is 


_  1  W)1  Bp 

2  a>  2un2  V 


(7) 


which  has  the  same  form  aa  the  usual  grad  B  drift  but  is 
quantitatively  different.  The  single  panicle  equations  of  mo¬ 
tion  are  ti,  *  —  o,/l„  u,  —  —  ot/2,  +  u,il„  tad  li, 
mo,/],.  Assuming  v,  ~  constant  and  that  the  gyration  velo¬ 
cities  are  v,  m  sin  kx  and  u,  m  Vu  cos  kx,  it  follows  that 


FIG.  4.  Electron  memorise  in  timer  wigglcr  *tth  1 
kO,  and  r,  •  0.4  cm.  Iti  a,  -  I  kO.  ,bi  B,  m  10  kO. 


bit/lm  cosh  kx 

n\-kWt 


!lcf)w  cosh  kx 

ai-kV,  ‘ 


18) 


for  3  '  >3 ‘t .  These  assumptions  are  reasonably  valid  for 
kx  50.8  and  for  B0  sufficiently  far  from  resonance  and  Bm 
small  enough  that  o.  <vt.  Furthermore,  to  insure  that 
cosh  kx ~ constant  over  the  orbit,  we  restnet  v._,  <v_t  =u, 
lie.,  J?0  <  kvt).  Then  usmg  Eqs.  (81  and  (6)  in  Eq.  i7),  we  ob¬ 
tain 


»(S 


kv,  {/)  mko.f 

ir0  &i-kV'f 


cosh  kx  srnh  kx. 


'91 


Table  I  compares  drift  velocities  for  varous  cases  from  the 
code  with  /  »  0  to  those  from  Eq.  19).  In  general,  the  agree¬ 
ment  is  very  good. 

Clearly,  the  gradient  drift  is  more  severe  when  Bw  is  a 
considerable  fraction  of  Bg.  However,  if  Bu  is  held  constant 
at  1  kG  and  B0  increased  to  4  kG,  a4\y  m  2.2)  is  not  signifi¬ 
cantly  reduced.  This  relative  insensitivity  to  B0  is  due  to  the 
increased  ol:  and  hence  a  larger  u4,  as  B0  approaches  the 
cyclotron  resonance  field  B,.  One  consequence  of  this  large 
drift  near  resonance  is  that  it  limits  the  degree  of  gain  en¬ 
hancement  achievable  through  the  magnetoresonance  ef¬ 
fect.* 

If 30 sufficiently  exceeds.#,,  the  drift  can  be  quite  small. 
For  y  m  2.2  and  Xw  m  3  cm.  the  resonant  magnetic  field  B, 
s7  kG.  With  Ba  m  g  kG  and  3,  «  1  kG.  the  particle  still 
drifts  to  the  wall  as  shown  in  Fig.  4ial.  The  dnft  is  now  in  the 
opposite  direction  to  that  below  resonance  since  ut  changes 
sign  with  B,  >  B,.  This  has  the  effect  of  changing  the  phase 
of  u,  oscillations  with  respect  to  those  of  B, ,  i.e..  v,  is  positive 
(for  x>  01  when  B,  is  a  minimum,  so  that  the  particle  drifts  in 
the  j-  y  direction.  When  Ba  is  increased  to  10  kG  for  this 
case,  the  drift  is  small  enough  that  the  electron  remains  con¬ 
fined  for  >  30  periods  as  shown  in  Fig.  4tbi.  The  confinement 
remains  very  good  when  the  self-fields  of  a  300- A  beam  are 

In  principle,  operation  of  FEL  experiments  with 
B0> B.  is  possible  and  has  been  demonstrated  with  a  high 
current.  y=3.5  beam  in  a  helical  wtggier.*  However,  com¬ 
peting  processes  such  as  the  cyclotron  maser1*  interaction 
can  in  general  produce  large  radiated  powers  at  frequencies 
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close  to  those  of  the  FEL  interaction  for  bean  energy  —  1 
MeV  and  B,  2  10-20  kG.  Consequently,  analysis  of  experi¬ 
mental  results  is  more  difficult.  Also,  arbitrarily  large  guide 
fields  are  not  possible  with  a  magnetically  shielded  diode 
simply  because  the  electrons  will  be  mirrored  by  the  field. 
The  equilibrium  beam  radius'6  in  such  a  case  depends  only 
on  y,  /.  2?o,  and  the  beam  emit  tan  ce  e,  so  for  particular  beam 
parameters  suitable  values  of  5,  are  limned.  The  required 
field17  in  kG  for  a  matched  beam  radius  R  in  cm  can  be 
written 


1.36/  _  H.56e; 
R  l&y  R* 


( 10) 


where  /  is  in  kA  and  f,  is  the  normalized  emittance  in 
rad  cm.  For  example,  if  f,  =0.14  rad  cm  (about  the  lowest 
value  expected  for  a  thermionic  cathode  beam  with  /~750- 
1000  Ai,'7  a  y  =  2.2  beam  with  R  =■  0.3  iC.51  cm  requires 
B0  —  5.5  121  kG  for  IS  800  A.  In  this  regime,  the  beam  is 
emittance  dominated  so  that  the  required  field  is  relatively 
insensitive  to  I.  In  principle,  smaller  radius  beams  could  be 
used  with  larger  Bo,  thereby  doubly  reducing  v4 .  However, 
experimentally  this  is  very  difficult  at  high  current  levels. 16 

To  analyze  the  effect  of  a  shielded  diode  on  propagation 
through  the  wiggler.  we  repeat  the  above  calculations  with 
an  initial  v,  corresponding  to  Peo  =*  0.  Note  that  I  *  1.75 
kA  in  this  case,  which  is  the  current  required  for  constant 
radius  propagation  in  only  the  solenoidal  field  with  these 
initial  conditions.  As  shown  in  Fig.  5,  the  electron  propa¬ 
gates  at  nearly  constant  radius  until  the  wiggler  amplitude 
becomes  large  enough  that  the  gradient  drift  begins  to  domi¬ 
nate.  Then  the  particle  is  lost  just  as  in  the  vt0  m  0  case. 

The  behavior  of  an  off-Axis  electron  in  a  linear  wiggler 
and  axial  guide  field  should  be  compared  to  that  when  a 
helical  wiggler  is  used. 13  We  approximate  the  helical  wiggler 
field  by6 


B* -b\z\BJ  -sin/kzfl  -r  j  *  J(3xr -r/|] 

-f  }fc:jrycoskz), 

(111 

B  J  =  b  |  cos  *z[  I  + 1  *  5(  jr  *  3/1] 

—  $  k  :xy  sin  kx  j . 

B*  =  -  kb\z)B„  [  1  -r  jk3!  jt  -y2)]!  x  cos/cr-y  sin  kz\. 

This  expansion  of  the  true  Bessel  function  expression  for  B " 
is  valid  for  fcr<  1.  Fnedland's  treatment13  of  electron  propa¬ 
gation  b  this  case  assumes  a  radially  uniform  wiggler,  and 
his  stability  condition  [Eq.  (2a)]  is  not  sringent  enough  when 
the  radial  variation  is  included.  For  example,  Friedland 
finds  stable  orbits  when  y  «  1.587,  B^  =4  cm,  i\0  =0, 
r0  =  0,  B0  =  1.26  kG.  and  B —  1.04  kG.  so  that  Eq.  I2al  is 
barely  satisfied,  and  we  can  duplicate  his  results  if  we  remove 
the  radial  variation  from  B  \  However,  with  the  B '  given  in 
Eq.  Ill),  we  find  that  the  wiggler  field  must  be  reduced  to 

—  625  G  to  obtain  stable  orbits. 

Although  the  radial  dependence  in  B  *  does  narrow  the 
allowable  range  of  operating  parameters,  stable  orbits  with 
(r)  s constant  are  achievable  with  a  helical  wiggler  in  a 
guide  field.  Electron  trajectories  in  a  helical  wiggler  for  the 
same  conditions  as  used  previously  with  a  linear  wiggler  are 
shown  in  Figs.  6  and  7.  Figure  6  is  to  be  compared  to  Fig.  2 
and  Fig.  7  to  Fig.  5.  In  both  cases,  the  electron  is  well  con¬ 
fined.  It  is  interesting  to  note  that  in  Fig.  6,  the  electron  born 
at  (  xjit  **  (0.4,0)  initially  E  X  B  drifts  in  the  —  9  direction, 
but  then  reverses  direction  as  the  wiggler  amplitude  in¬ 
creases.  Since  grad  B 6  is  radial,  the  grad  B  6  drift  is  in  the 

—  6  direction  and  does  not  lead  to  the  beam  expansion  ob¬ 
served  with  the  linear  wiggler. 

The  linear  wiggler  drift  described  here  imposes  addi¬ 
tional  constraints  on  the  parameters  of  an  FEL  experiment. 
Obviously,  the  beam  radius  must  be  kept  as  small  as  possible 


FIG.  5  Electron  uajecionei  in  linear  *mkr  wuh  •  0.  /  «  I  75  kA. 
r-2.2.  J„-2kG.  8,  -  l  kG.r,-0.6cm. 


FIG.  6.  Electron  trajenonw  with  helical  wiggler  tad  »  —  2.2, 8^  «*  2kG. 
8,  »  I  kG.  Xm  m  3  cm.  /  »  230  A,  r,  m  0  4  cm.  aod  —  0 
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FIG.  7.  Electron  trajectories  in  helical  wiggler  with  i°,0  -  0.  /  -  l  ? 3  kA 
y  m  xx  B%  m  2  kG.  3m  »  l  kG.  and  r,  —  0.6  an. 
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to  mimmiM  particle  loss  from  the  edge  of  the  beam.  Prelimi¬ 
nary  experimental  results  by  our  group  with  a  field  im¬ 
mersed.  apertured  source  indicate  that  particle  losses  can  be 
kept  acceptably  small  in  this  way.  Also,  if  y  is  large  enough 
that  B ,  >Bo>B,  can  be  satisfied  for  relatively  large  Bm ,  then 
the  drift  can  be  kept  small  while  achieving  acceptably  large 

Finally,  it  should  be  noted  that  the  drift  arises  from  the 
asymmetry  of  the  linear  wiggler  and  the  corresponding  ab¬ 
sence  of  focusing  forces  in  the  direction  perpendicular  to  the 
wiggler  field.  Therefore,  it  should  be  possible  to  stabilize  the 
drift  by  imposing  an  additional  focusing  force  in  that  direc¬ 
tion.  For  example,  preliminary  results  indicate  that  electron 
propagation  through  a  “square"  or  symmetrized  linear 
wiggler  is  very  stable.  Such  a  wiggler  has  an  additional  com¬ 
ponent  B  ‘f  —  cosh  ky  cos  kz  and  a  corresponding  addition 
to  B[  of  —  B„  smh  ky  sm  kz.  An  electron  trajectory 
through  such  a  wiggler  is  shown  in  Fig.  3. 

In  conclusion,  a  gradient  drift  has  been  shown  to  exist 
for  a  linear  wiggler  in  an  axial  guide  field.  The  drift  can  be 
substantial  with  small  or  large  beam  current  in  some  param¬ 
eter  ranges  for  a  wide  range  of  initial  conditions.  However, 
the  advantages  of  a  linear  wiggler  are  sufficient  in  many 
cases  to  either  limit  operation  to  a  “stable”  parameter  regime 
of  to  impose  additional  focusing  forces  to  stabilize  the  drift. 
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A  BROADBAND  HIGH  POWER  MILLIMETER  TO 
CENTIMETER  SPECTROMETER 


I.  INTRODUCTION 


Over  Che  last  several  years  there  has  been  considerable  interest  in 
generating  high  power  microwaves  from  relativistic  electron  beams. *  One  of 
the  most  important  parameters  to  be  determined  in  such  experiments  is  the 
power  spectrum.  Interferometers  have  been  used  for  spectral  measurements,^ 
but  they  are  only  useful  over  a  very  narrow  band  of  frequencies.  An  improved 

•5 

bandwidth  is  obtained  with  a  grating  spectrometer  but  noise  and  expense  then 
become  issues.  Dispersive  lines  have  also  been  used  but  they  are  impractical 
for  long  pulses  (>  1  us).  Also,  all  of  the  above  mentioned  devices  must  be 
calibrated  for  absolute  power.  In  the  beginning  of  a  microwave  experiment 
quick  but  less  precise  information  about  the  entire  power  spectrum  are  of 
primary  importance.  To  fit  this  purpose,  we  have  developed  a  simple  and 
inexpensive  spectrometer  which  is  insensitive  to  noise  and  can  dire"~y 
measure  selected  wavelengths  at  full  power. 


II.  SPECTROMETER  DESIGN 


Figure  I  shows  a  schematic  of  Che  spectrometer.  Radiation  is  collected  in 
the  horn  on  the  left.  It  then  travels  through  a  high  pass  filter  and  expands 
through  the  horn  on  the  right  to  a  collimating  or  weakly  focusing  lens.  The 
microwaves  are  reflected  from  a  metallic  boundary  which  is  located  inside  the 
gas  filled  chamber.  Gas  breakdown  occurs  at  cne  spatially  periodic  peaks  in 

Che  electric  field  which  result  from  interference  between  the  reflected  and 
Manuicript  approved  May  23,  1983. 
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incident  waves.  Gas  breakdown  occurs  at  about  half  wavelength  intervals 
apart.  The  wavelength  is  related  to  the  spot  spacing  d  and  the  angle  of 
incidence  (see  Fig.  1)  by 

X  •  2  d  cos  a. .  (1) 

i 

Thus,  by  measuring  from  a  photograph  of  the  breakdown  region  the  spot 
spacing  and  diameter,  and  knowing  the  gas  pressure,  then  the  electric  field, 
power  and  wavelength  can  be  determined.  A  relationship  between  breakdown 
electric  field  and  gas  pressure  will  be  given  later.  It  should  be  noted  that 
the  chamber  length  does  not  determine  the  wavelength  since  the  radiation  is 
reflected  from  only  one  highly  reflecting  boundary.  For  low  power  microwaves 
^  the  focusing  lens  increases  the  electric  field  to  facilitate  gas  breakdown, 
j  The  aluminum  horns  are  conical  with  an  inside  diameter  that  is  reduced 
from  12.5  cm  to  5  cm.  An  aluminum  taper  matches  the  end  of  the  horns  to  the 
high  pass  filter.  Focusing  is  accomplished  with  one  of  three  lucite  lenses 
s  having  focal  lengths  of  12,  24,  and  48  cm.  Focusing  inside  and  outside  the 
chamber  end  plate  has  been  tried;  however,  the  clearest  interference  patterns 
are  observed  when  focusing  is  outside  the  30  cm  long,  14  cm  i.d.  lucite 
chamber. 

Ill .  EXPERIMENTAL  RESULTS 

In  the  present  experiment  microwaves  are  generated  when  a  130  A,  550  keV, 
2  us  long  electron  beam  is  passed  through  a  "wiggler"  magnetic  field.  The 


rr 


wiggler  field  is  produced  from  Che  eddy  currents  induced  in  a  set  of  spatially 
periodic  conducting  rings  when  an  axial  magnetic  field  is  pulsed.  In  this 
particular  case  the  rings  are  made  of  aluminum  and  have  an  inside  diameter  of 
7.5  cm,  a  radial  thickness  of  0.6  cm  and  an  axial  extent  of  1  cm.  The  ring 
period  is  6  cm  and  the  composite  ring  structure  and  axial  guide  field  are  200 
cm  long. 

Figure  2  is  a  time  integrated  photograph  of  the  resulting  interference 
pattern.  Microwave  radiation  has  entered  from  the  left  of  the  photograph, 
where  the  lens  is  located  (but  cannot  be  seen),  and  is  reflected  off  a  copper 
plate  located  at  the  right  of  the  photograph  (arrow  indicates  copper  plate). 
The  white  spots  are  due  to  light  produced  from  the  gas  breakdown.  Type  57 
(ASA  3000)  Polaroid  film  was  used  with  a  Graphflex  camera  (f/4.5)  to  obtain 
these  results.  The  gas  density  was  selected  experimentally  so  that  enough 
light  was  available  to  be  recorded  photographically  but  not  high  enough  to 
allow  microwave  reflections  from  the  plasma. 

From  Fig,  2  the  spot  spacing  is  measured  to  be  2.25  cm  which  corresponds 
to  a  frequency  of  6.67  GHz.  The  breakdown  electric  field  is  estimated  from 
the  semi-empirical  formula 

Eb  -  AP{l  +  (a/PX)2}1/2,  (2) 

where  A  ■  3000  Vm~*Torr_1  and  a  »  0.9  Torr-m  for  air  or  nitrogen,  P  is  the  gas 
pressure  in  Torr,  and  X  is  the  wavelength  in  meters. 

Equation  (2)  has  been  derived  semi-empirically  from  McDonald's  precision 
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data4  and  agrees  in  che  worse  case  Co  within  302.  this  formula  is  valid  on 
the  high  pressure  side  of  the  breakdown  curve,  i.e.  for  PX  >  0.26  Torr-meters, 
pulse  lengths  >1  us,  low  repetition  rate  (<  100  pps ) ,  and  for  no  electron 
collisions  with  the  chamber  walls.  For  che  results  of  Fig.  2  the  measured 
pressure  is  25  Torr  and  the  wavelength  is  4.5  cm,  corresponding  to  an  electric 
field  of  about  1  kV/cm.  From  the  spot  diameter  of  2.5  cm  the  power  is 
estimated,  for  a  plane  wave,  to  be  8  kW. 

Wavelength  discrimination  is  accomplished  by  utilizing  the  pressure  and 
wavelength  dependence  of  Eq.  2.  That  is,  at  a  given  pressure  a  shorter 
wavelength  will  require  a  higher  electric  field  in  order  to  achieve 
breakdown.  As  s faced  before,  Eq.  2  is  valid  for  PX  >0.26  Torr~meter.  In  che 
ocher  limit,  i.e.,  for  PX  <  0.26,  it  will  be  shown  later  that  short 
wavelengths  cannot  be  resolved  as  a  result  of  diffusion. 

Discrimination  is  not  possible  for  all  wavelengths  in  certain  electric 
field  versus  wavelength  distributions.  If  any  two  spectral  lines  in  a 
discrete  distribution  satisfy  che  relation 


AE  -  Ei  (l- 


,  Aa  s2_.  ( 2-a)  1/2  , 

^  +  (X^°  a - -  f 


(3) 


( 1-a) 


then  discrimination  is  not  possible.  In  Eq.  (3),  AX  -  X 
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X2.  X1  >  x2, 


AF  ■  E^  -  E^,  E ^  <  E^  and  a  ■  AX/X^,  where  E^  i  is  the  electric  field  for 
wavelength  X^  Eq.  3  simply  states  chat  che  two  spectral  lines  satisfy  Eq.  2 


at  che  same  pressure.  Even  if  Eq.  (3)  is  satisfied  there  is  a  way  to 
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circumvent  this  difficulty,  provided  the  wavelengths  are  sufficiently 
separated.  In  such  a  case,  a  high  pass  filter  can  be  used  to  eliminate  the 
longer  wavelength. 

There  are  two  more  points  Co  consider  regarding  Eq.  3.  If  AE  is  greater 
than  the  right  hand  side  of  Eq.  3  then  the  pressure  can  be  adjusted  so  that  £2 
is  not  large  enough  to  breakdown  the  gas.  However,  if  X E  is  less  chan  the 
right  hand  side  of  Eq.  3  then  the  pressure  can  be  adjusted  to  prevent 
breakdown  by  Ep  To  resolve  every  wavelength  in  a  continuous  distribution 
E  (1),  it  follows  from  Eq.  3  chat 


d£ 

dX 


> 


(4) 


for  every  X. 

In  order  to  verify  that  Che  observed  breakdown  pattern  is  due  to  an 
interference  effect,  the  boundary  conditions  on  the  incident  wave  are  changed 
with  the  results  shown  in  Fig.  3.  In  each  case  microwaves  are  focused  by  a 
lens  located  at  the  left  of  the  photograph  to  a  region  in  front  of  the 
boundary  at  the  right.  By  comparing  Figs.  3a  and  3b,  it  is  clear  that  the 
total  electric  field  has  increased  by  increasing  the  surface  reflectivity.  In 
Fig.  3a  the  boundary  is  composed  primarily  of  graphite  whereas  in  Fig.  3b  the 
material  is  lucite  (index  *  1.63).  It  should  be  noted  that  the  thickness  of 
the  lucite  was  selected  not  to  satisfy  the  resonant  reflection  condition.  In 
Fig.  3c  an  absorbing  surface  is  used  to  eliminate  the  reflected  wave.  The 
interference  pattern  is  absent  but  a  single  region  of  breakdown  remains. 


I 

I 
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This  residual  breakdown  region  results  from  Che  method  used  to  absorb  the 
radiation.  The  absorber  was  constructed  from  Eccosorb  (Emmerson  and  Cuming, 
Inc.)  in  the  shape  of  a  "V"  and  oriented  so  chat  microwaves  enter  the  opening 
of  Che  ''V".  This  technique  results  in  substantial  reduction  of  the  reflected 
wave.  The  gas  breakdown  now  occurs  within  the  "V”  shaped  absorber  region. 

IV.  LIMITATIONS 

There  are  several  limitations  which  are  inherent  in  the  above  type  of 
spectrometer.  Practical  considerations  limit  both  the  wavelength  and  the 
electric  field  strength  that  can  be  measured.  Long  wavelength  ( > 1 0  cm) 
operation  is  impractical  while  trying  to  maintain  the  geometrical  optics 
condition.  Short  wavelengths  (<1  mm)  are  bounded  by  the  ability  to  resolve 
the  spot  spacing.  Resolution  is  maximized  by  operating  with  a  highly 
reflective  chamber  end  plate  and  a  high  gas  pressure,  which  limits  electron 
diffusion.  High  pressures,  however,  require  very  high  electric  fields  (>  100 
kV/cm)  for  submillimeter  resolution.  At  low  electric  fields  (<  100  V/cm)  a 
limit  is  set  by  the  ability  of  the  electrons  to  gain  sufficient  energy  to 
breakdown  Che  gas. 

In  order  to  make  the  wavelength  measurements  two  conditions  must  be 
satisfied:  the  breakdown  region  must  be  visible  and  the  wavelength  must  be 
resolvable.  The  dashed  curve  in  Fig.  4  shows  the  total  RMS  electric  field  of 
the  standing  wave  divided  by  the  peak  incident  electric  field.  When  the 


breakdown  electric  field  is  above  the  peak  total  RMS  electric  field  (top  line 
in  Fig.  4),  the  gas  does  not  breakdown.  Also  wavelength  resolution  is  not 
possible  if  the  breakdown  electric  field  is  below  the  minimum  total  RMS 
electric  field  (as  in  the  bottom  line),  since  breakdown  occurs  everywhere. 
Resolution  is  maximized  when  the  breakdown  field  (solid  line  in  Fig.  4)  is 
just  slightly  below  the  peak  total  RMS  electric  field,  so  that  the  length  of 
the  breakdown  region  is  considerably  smaller  than  the  wavelength.  The  length 
of  the  breakdown  region  can  be  determined  from  the  total  RMS  electric  field 
for  two  plane  waves  traveling  in  opposite  directions, 

E  *  {E,2/2  +  E  2/2  -  E  E  cos(2kx) >i/2,  (5) 

rms  i  r  l  r 

where  Ei  and  Er  are  the  incident  and  reflected  electric  field  amplitudes  for  a 
plane  polarized  wave  incident  normally  on  a  partially  reflecting  surface 
located  at  X  *  0,  and  k  »  2* A  is  the  wavenumber.  From  Eq.  5  the  length  of 
the  breakdown  region  Is  given  by 


where  E^  -  Er  <  2  Ej,  <  E^  +  Er  for  resolvable  breakdown.  An  important 

observation  here  is  that  for  a  poorly  reflecting  boundary,  the  minimum  E^g 
approaches  Che  maximum  E^g  and  breakdown  will  occur  nearly  everywhere  or  not 
at  all,  and  again  a  wavelength  measurement  cannot  be  made.  Figure  5  is  a  plot 
of  the  ratio  of  the  total  RMS  electric  field  to  the  peak  incident  electric 
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field  versus  position  for  three  different  values  of  E  / Ej.  It  is  clear  that  a 
highly  reflecting  boundary  is  desirable  for  wavelength  resolution.  It  should 
be  mentioned  that  a  shallow  focusing  angle  is  necessary  in  order  to  keep  the 
reflected  wave  amplitude  comparable  to  the  incident  amplitude.  This  accounts 
for  the  fact  that  clearer  interference  patterns  were  observed  when  the 
radiation  was  focused  beyond  the  end  plate. 

So  far,  only  the  initial  breakdown  region  has  been  considered.  However, 
electron  diffusion  will  increase  Che  length  of  this  initial  region,  and  hence 
further  restrict  the  ability  to  resolve  a  given  wavelength.  Diffusion 
contributes  a  length  given  by 


xd  -  (Dt)1/2,  (7) 

where  D  is  the  diffusion  coefficient  and  C  is  the  time.  In  order  to  relate 
the  diffusion  length  to  properties  of  the  gas  and  the  electrons,  the  diffusion 
coefficient  must  be  evaluated.  From  elementary  kinetic  theory  the  diffusion 
coefficient  for  a  Maxwellian  distribution  of  electrons  can  be  shown  to  be  in 
MRS  units 


D  -  2eU/(3mf  ), 
c 


(8) 


where  m  i3  the  electron  mass,  eU  is  Che  mean  electron  energy  (U  is  in  eV),  and 
f  is  Che  electron-neutral  collision  frequency.  From  Reference  4,  a  simple 
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relation  is  available  chat  connects  the  collision  frequency  to  the  pressure 
and  is  given  by, 

fc  -  bP,  (9) 

Q 

where  b  ■  5.3  x  10  /sec-Torr  for  air  or  nitrogen.  Also,  from  Reference  4  the 
measured  electron  energy  is  typically  a  couple  of  eV.  When  Eqs.  (7)-(9)  are 
combined  it  becomes  clear  that  for  long  pulses  and  short  wavelengths  a  high 
pressure  is  required  in  order  to  achieve  wavelength  resolution.  Higher 
pressures  will  then  require  high  electric  fields  to  obtain  gas  breakdown. 

Figure  6  demonstrates  the  combined  effects  of  diffusion  and  breakdown 
field,  on  wavelength  resolution.  In  Fig.  6a  the  image  is  cloudy  which  makes  a 
wavelength  determination  dubious.  By  raising  the  pressure  (from  2  to  20  torr) 
the  image  becomes  clear  as  shown  in  Fig.  6b.  A  lucite  wall  was  used  as  a 
reflector  in  Fig.  6a  and  6b.  At  20  torr  the  one  dimensional  diffusion  length 
is  from  Eq.  (7)  equal  to  0.26  cm.  Since  diffusion  occurs  in  both  directions 
from  the  center  of  the  breakdown,  the  calculated  spot  thickness  will  be  twice 
the  one  dimensional  diffusion  length,  namely  0.52  cm.  The  measured  spot 
thickness  (Fig.  6b)  is  0.5  cm  which  is  in  very  good  agreement  with  the 
calculated  value.  At  2  torr  the  two  directional  diffusion  length  is  1.6  cm 
which  accounts  for  smearing  of  discrete  spots.  However,  there  may  be  an 
additional  effect  which  clouds  the  image  of  Fig.  6a.  Plasma  motion  has  been 
observed  when  the  electric  field  is  substantially  above  the  breakdown  electric 
field. ^  The  conditions  for  Fig.  6a  are  such  that  the  electric  field  is  three 


times  Che  breakdown  field  which  is  sufficient  for  plasma  motion  to  occur,  and 
since  this  is  a  time  integrated  photograph  plasma  motion  could  account  for 
some  of  the  image  clouding. 

In  general,  Che  breakdown  field  versus  pressure  exhibits  a  minimum.  This 
fact  combined  with  diffusion  imposes  a  limit  on  the  minimum  resolvable 
wavelength.  This  electric  field  minimum  corresponds  to  a  maximum  transfer  of 
energy  to  the  electrons  from  the  electromagnetic  field,  which  occurs  when  the 
time  becween  collisions  is  some  fraction  of  half  the  wave  period.  That  is 

l/fc  -  T/( 21 ) ,  (10) 

where  T  is  the  wave  period  and  i  is  a  number  that  depends  on  the  electron- 
neutral  scattering  angle.  For  a  180°  angle,  i  »  1;  for  smaller  angles,  i 
becomes  larger  chan  1.  For  electrons  with  a  low  energy  (<  5  eV)  the 
scattering  angle  is  large,  thus  i  will  be  of  order  1.  By  combining  Eqs.  (9) 
and  (10)  a  relation  for  the  pressure  at  the  minimum  breakdown  field  can  be 
obtained.  When  this  is  done  and  the  results  are  compared  with  the  data  of 
Reference  4,  then  i  *  1.86  is  obtained  for  air. 

At  the  minimum  breakdown  field  the  shortest  resolvable  wavelength  can  now 
be  calculated.  By  setting  x^  ■  X/4  in  Eq.  (7)  and  solving  Eqs.  (7),  (8),  and 
(10)  for  Che  wavelength,  the  shortest  resolvable  wavelength  is  found  to  be 

\  *  16  eUt/(3  mic),  (11) 
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where  c  is  Che  speed  of  light.  This  result  represents  a  rather  fundamental 

relationship  for  any  gas  and  only  a  knowledge  of  the  mean  electron  energy  and 

scattering  angle  constant  is  required.  But  since  the  mean  electron  energy  for 

any  gas  is  less  chan  5  eV,  the  scattering  angle  will  be  large  and  i  will  be  of 

order  1.  For  Che  present  experiment,  t  ■  2  us,  i  *  1.86,  and  U  ■  1.5  eV,  so 

chat  X  •  0.5  cm. 
s 

If  it  is  desirable  to  resolve  shorter  wavelengths,  then  shorter  pulses 
and/or  higher  electric  fields  could  be  used  as  shown  in  Fig.  7.  The  dotted 
curve  is  a  plot  of  the  diffusion  relation  (Eq.  7)  with  x^  *  X/4,  (J  »  1.5  eV, 
and  c  *  2  us.  The  solid  lines  are  plots  of  the  higher  pressure  side  of  the 
breakdown  curve  (Eq.  2)  for  air  at  different  electric  field  strengths.  Since 
Xd  <  X/4  is  required  for  resolution  when  diffusion  is  considered,  only 
wavelengths  above  the  diffusion  curve  can  be  resolved.  Also,  for  a  particular 
electric  field,  X  must  lie  on  the  breakdown  curve  in  order  for  the  gas  to 
breakdown  only  at  Che  standing  wave  peaks.  Thus,  the  minimum  resolvable 
wavelength  at  any  particular  field  is  given  by  the  intersection  of  the 
breakdown  and  diffusion  curves.  When  the  curves  do  not  intersect,  diffusion 
is  not  a  limiting  factor  and  the  minimum  wavelength  is  just  determined  by  the 
field  strength.  It  should  be  noted  that  very  large  electric  fields  are 
required  as  one  approaches  submillimeter  wavelengths. 

An  expression  for  the  minimum  resolvable  wavelength  in  air  versus  electric 
field  can  be  found  by  eliminating  the  pressure  in  Eq.  (2),  (7),  (8)  and  (9) 
and  is  given  by 
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(12) 


(13) 


This  very  slow  dependence  on  electric  field  gives  a  practical  limit  to  the 
minimum  resolvable  wavelength. 

In  order  to  use  this  spectrometer  for  longer  wavelengths  a  larger  chamber 
and  lens  would  have  to  be  built  in  order  to  preserve  the  geometric  optics 
condition.  This  then  sets  a  practical  limit  on  the  utility  of  this 
spectrometer  for  long  wavelengths.  Effects  from  electrons  hitting  the  chamber 
wall  and  ion  motion  occur  at  such  low  frequencies  that  they  are  not  a  limiting 
factor  in  view  of  the  diffraction  condition. 


V.  CONCLUSION 


In  summary  then,  a  simple  spectrometer  has  been  described  which  is  capable 
of  both  high  power  and  broadband  wavelength  measurement.  Wavelength 
resolution  has  been  shown  to  depend  on  the  reflectivity  of  the  chamber  end 
plate  and  electron  diffusion.  Electron  diffusion  and  the  characteristic 
breakdown  curve  limits  the  resolvabilitv  of  submillimeter  wavelengths  by 
requiring  very  high  electric  fields  (>  iOO  kV/cm).  At  high  fields,  the 
minimum  resolvable  wavelength  is  proportional  to  the  reciprocal  square  root  of 
the  breakdown  field.  This  slow  dependence  on  electric  field  gives  a  practical 
limit  to  both  electric  field  and  submillimeter  wavelength  measurements. 


12 


Figure  I.  Schematic  of  the  spectrometer. 


rms 


Figure  5.  The  effect  of  Che  boundary  reflectivity  on  the  total  RMS  electric 


field.  Reflected  to  incident  wave  amplitude  is  indicated  in  the 
upper  right  hand  corner.  Low  reflectivity  greatly  diminishes  the 
ability  to  select  the  proper  breakdown  condition  for  wavelength 
measurement . 
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P  (Torr) 

Shortest  resolvable  wavelength  at  various  electric  field 
strengths.  The  dotted  curve  is  the  diffusion  length  equated 
to  X/4,  and  the  solid  curves  are  the  characteristic  breakdown 
curves  for  different  electric  field  strengths.  Minimum  resolvable 
wavelength  is  given  by  the  intersection  of  the  diffusion  and 
breakdown  curves ,  or  simply  by  the  bottom  of  the  breakdown  curves 
when  the  two  curves  do  not  intersect.  The  curves  are  for  air  with 
a  diffusion  time  of  2  usee  and  an  electron  energy  of  1.5  eV. 
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ELECTRON  DRIFT  IN  A  LINEAR  MAGNETIC  WIGGLER  WITH  AN 
AXIAL  GUIDE  FIELD 


Recently,  there  has  been  much  interest  in  free  electron  lasers  (FEL's) 
which  use  spatially  varying  magnetic  fields  to  modulate  a  relativistic 
electron  beam  (REB)  .1-^  Three  types  of  magnetic  wigglers  have  been  widely 
used:  helical  (usually  produced  by  a  bifilar  helical  current  winding^’^), 

Q 

linear  (produced  by  an  array  of  permanent  magnets  or  by  linearly  alternating 
windings  ) ,  and  radial  u  (produced  by  a  series  of  spaced  conducting  or 
ferromagnetic  washers11  immersed  in  an  axial  field  or  by  a  series  of 
alternating  coils  ) .  Of  these,  only  the  first  two  produce  a  perpendicular 
field  on  axis  and  are  therefore  suitable  for  use  with  a  small  diameter  beam  of 
solid  cross-section.  In  this  paper,  we  will  analyze  the  propagation  of  a 
solid,  high  current  REB  through  linear  and  helical  wigglers  with  a 
superimposed  axial  guide  field. 

FEL  experiments  fall  into  two  categories  depending  on  the  beam  current. 
Compton  regime  FEL's1’^’^  have  used  high  energy  (10' s  of  MeV)  ,  low  current 
(<  I  A)  beams  while  Raman  FEL's  have  used  lower  energy  (~  1  MeV) ,  high 
current  (>  1  kA)  beams.  Helical  wigglers  have  been  used  in  both  current 

Q 

regimes,  but  until  recently  linear  wigglers  were  used  only  in  the  Compton 
regime.  In  this  case,  the  linear  wigglers  consist  of  permanent  magnets  and 
there  is  no  guide  field.  However,  for  high  current  beams,  a  guide  field  is 
required  to  contain  the  beam.  Although  there  are  advantages  to  a  linear 
wiggler  from  the  standpoint  of  ease  of  assembly  and  versatility  (e.g., 
changing  the  periodicity  or  capering  the  period  and/or  field  amplitude  for 
efficiency  enhancement  )  we  will  show  that  there  is  no  equilibrium  for  off- 
axis  particles  when  a  beam  is  propagated  through  a  linear  wiggler  in  an  axial 
field.  The  particles  drift  out  of  the  wiggler  unless  additional  focusing 
forces  are  provided. 
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Various  authors  have  considered  particular  cases  of  electron  motion 
through  magnetic  wigglers.  These  have  typically  been  single  particle  (low 
current)  calculations.  Blewett  and  Chasman^  considered  motion  of  high  energy 
electrons  (~  24  MeV)  through  a  helical  wiggler  and  found  stable  helical  orbits 
with  superimposed  betatron  oscillations.  Friedland^  has  treated  the  case  of 
lower  energy  electrons  (  ~  300  keV)  in  an  idealized  radially  uniform  helical 
wiggler  with  a  superimposed  axial  guide  field.  He  shewed  that  various  stable 
trajectories  were  possible  and  derived  stability  criteria  relating  the  allowed 
wiggler  and  guide  field  strengths  to  bean  energy  and  wiggler  period  A^. 

These  "stable"  regions  are  given  by 


o  o 


(2a) 


and 


flo  >  ck’ 


(2b) 


where  BQ  is  the  axial  guide  field,  *  eBQ/'fm  is  the  corresponding  cyclotron 
frequency,  v  is  Che  electron  velocity,  and  k  -  Physically,  these 

conditions  3Cem  from  the  resonance  in  Che  perpend icular  velocity  of  an 
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electron  in  the  wiggler  and  guide  fields: 
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Freund  and  Drobot^ 


have  considered  this  case  further  and  also  conclude  chat 


stable  trajectories  with  nearly  constant  axial  velocities  and  relatively  large 
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wiggler  amplitudes  are  possible  when  Q  <<  ck.  This  is  consistent  with 
condition  (2a). 

The  present  analysis  employs  a  relatively  simple  computer  code  which 
solves  the  equations  of  motion  of  an  electron  in  any  electric  and  magnetic 
field  configuration  using  a  fourth-order  Runge-Kutta  method.  Self  electric 
and  magnetic  fields  are  calculated  by  assuming  that  the  electron  is  at  the 
edge  of  an  azimuthally  symmetric  beam  of  current  I,  so  that  the  self  fields 
can  be  written  as  (tnks) 
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-  I 

2ite  rv 
o  z 
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(4) 


where  r(t)  is  the  electron  radius.  The  axial  self  fields  are  neglected. 


For  the  particular  cases  considered  here,  the  external  magnetic  field 
consists  of  a  solenoidal  field  produced  by  a  15.3-cm-I.D.  x  2-m-long  solenoid 
together  with  a  wiggler  that  begins  in  the  uniform  portion  of  the  solenoidal 
field.  The  wiggler  field,  which  may  be  either  helical  or  linear,  rises 
adiabatically  over  ten  periods  and  then  oscillates  with  constant  amplitude. 
The  envelope  enclosing  the  wiggler  amplitude  is  given  by 


b(z) 


i  Kife-)2  +  0 '  *  ‘ 10 
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z  >  10  X. 


(5) 
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Here,  z  is  Che  axial  distance  from  Che  beginning  of  che  wiggler.  This 
variation  closely  fits  che  linear  wiggler  used  in  the  NRL  induction  linac 
FEL.  The  field  on  axis  from  this  wiggler  is  plotted  in  Fig.  1  along  with  the 
envelope  equation  (5). 

The  linear  wiggler  field  components  are 

l 

3  ■  b(z)B  cosh  kx  cos  kz 

X  w 

-  0  (5) 

■  -b(z)B  sinh  kx  sin  kz . 
z  w 

Q 

where  Bw  is  the  peak  wiggler  field  on  axis.  For  the  MRL  linear  wiggler, 

which  has  \  *  3  cm  and  winding  layers  spaced  by  3.2  cm,  these  expressions  are 
w 

valid  for  r  <  1  cm. 

Ue  will  be  primarily  concerned  with  relatively  low  energy  (~1  MeV) ,  high 
current  (~kA)  beams.  Of  particular  interest  are  the  effects  of  self-fields, 
guide  field  and  wiggler  amplitude,  and  initial  beam  conditions  corresponding 
to  a  field- immersed  (vgQ  ■  0)  or  a  shielded  source  (PQq  ■  0),  where  v^q  is 
the  intitial  azimuthal  electron  velocity  and  ?gQ  the  initial  canonical 
angular  momentum.  If  che  axial  self  magnetic  field  is  neglecced,  PgQ  -  0 
implies  vg  ■  erB2/(2Ym)  in  a  uniform  field  Bz. 

First,  we  consider  a  beam  with  Y  -  2.2,  I  *  250  A,  SQ  -  2  kG,  Bw  -  1  kG, 
and  Aw  ■  3  cm.  Note  chat  Eq .  (2a)  would  require  3W  <  3  kG  for  stability  in 
this  case  with  a  helical  wiggler.  The  x-y  trajectories  of  electrons  injected 
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envelope  b(z)  from  K .  (5) 


with  v,  *0  and  r  •  0.4  cm  into  a  linear  wiggler  are  shown  in  Fig.  2.  Each 
lo  o 

particle  initally  begins  to  E  x  3  drift  in  the  self  fields,  but  when  it 
reaches  a  region  of  large  wiggler  field  it  begins  drifting  in  the  in¬ 
direction.  In  each  case  the  electron  reaches  an  imaginary  wall,  located  at  r 
*  1.2  cm,  after  only  partially  traversing  the  wiggler.  It  should  be  noted 
that  the  assumptions  used  to  calculate  the  self  fields  become  invalid  as  the 
beam  distorts.  Consequently,  to  determine  if  these  self  fields  are 
responsible  for  this  drift,  the  calculations  are  repeated  with  1*0  (Fig. 

3) .  A  third  trajectory  is  also  plotted  for  Bw  *  0.5  kG.  The  electron  that 
originates  on  the  y-axis  is  well  confined,  but  electrons  off  the  y-axis  again 
drift  to  the  wall,  with  a  velocity  much  higher  for  Bw  =  1  kG  (<v^>  *  *047  c) 
than  for  Bw  *  .5  kG  ( <v^  >  -  .011  c) . 

This  behavior  can  be  explained  by  the  increasing  gradient  in  the  wiggler 
field,  and  consequently  the  emergence  of  a  significant  axial  field  component, 
as  the  distance  from  the  y  axis  increases.  Although  the  present  drift  arises 

-V 

from  the  gradient  in  Br,  and  is  in  the  direction  of  7B  x  B  ,  it  is 
a  v  wo 

quantitatively  different  from  the  usual  guiding  center  approximation  because 
the  field  variation  over  one  gyroperiod  is  so  extreme.  For  example,  when  kx  - 

l 

.8  (x  *  .4  cm  in  the  present  case) ,  3z  varies  from  +  Bw  to  -  Bw  over  one 

period.  However,  the  physical  mechanism  is  the  same  as  for  the  usual  gradient 

drift;  i.e.  ,  the  gyroradius  in  the  part  of  the  orbit  where  32  is  a  minimum,  or 
where  .1x11  is  a  minimum  in  Fig.  (3),  is  larger  than  where  Br  is  a  maximum.  The 
addition  of  self  fields  merely  imposes  an  additional  E  x  B  rotation  on  the 

obit,  so  that  an  electron  that  originates  on  the  y-axis  (and  is  therefore 

confined  when  1*0)  begins  to  drift  Into  a  region  of  increasing  B^ 

Therefore,  it  actually  is  "lost”  sooner  chan  an  electron  on  the  x-axis  which 
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x  (cm ) 


Fig.  2 


Ftg.  3  -  Electron  trajectories  in  linear  wiggler  with  *  3  cn,  3Q  ■  : 
Y  “  2.2,  rQ  -  .4  cm,  and  I  -  0.  3y  »  1  kG  for  particles  a  and  b,  and  Bw 
kG  for  particle  G. 


+  +  J £ 

initially  E  x  B  drifts  into  a  region  of  smaller  Bz 


An  approximate,  empirical  expression  for  the  drift  has  been  found  which  is 
in  quite  good  agreement  with  the  code  results  when  ^w.<  H  <  kvT  (below  the 
cyclotron  resonance).  In  terms  of  the  wiggler  gradient,  the  expresson  is 


-  1  7(BZ)  Bo 

Vl“  2  a-  2(bV  X  B- 


n 


which  has  the  same  form  as  the  usual  grad  B  drift  but  is  quantitatively 


different.  The  single  particle  equations  of  motion  are  v  =■  v  S»  , 

x  y  z 


-  v  SI  +  v  SI  ,  and  v 
2  X  X  2  2 


Assuming  v^  *  constant  and  that  the 


gyration  velocities  are  v  *  v  sin  kz  and  v 

y  iy  x 


v,  cos  kz ,  it  follows  that 
lx 


v  , 

iy 


V 

z 


kv  SI  coshkx 
z  w 


k2v 


v 


lx 


v 

z 


£1  SI  cosh  kx 
o  w 


,  2  2 
k  v 

z 


(3) 


l  l 

for  »  B^  .  These  assumptions  are  reasonably  valid  for  kx  <  .3  and  for 
B  sufficiently  far  from  resonance  and  Bu  small  enough  that  v  <<  v 
Furthermore,  to  insure  that  cosh  kx  *  constant  over  the  orbit,  we  restrict 
v,x  <<  -  vjL  (i.e.  ,  SJQ  <  kv2)  .  Then  using  Eqs.  (8)  and  (6)  in  Eq.  (7), 

we  obtain 
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- z - z — r— r-  cosh  kx  sinh  kx. 
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Table  1  compares  drift  velocities  for  various  cases  from  the  code  with  I  ■  0 
to  those  from  Eq.  (9).  In  general,  the  agreement  is  very  good. 

Clearly,  the  gradient  drift  is  more  severe  when  Bw  is  a  considerable 

fraction  of  BQ.  However,  if  Bw  is  held  constant  at  1  kG  and  BQ  increased  to  4 

kG,  v^(y  -  2.2)  is  not  significantly  reduced.  This  relative  insensitivity 
to  Bq  is  due  to  the  increased  v^  ,  and  hence  a  larger  v^,  as  BQ  approaches  the 

cyclotron  resonance  field  3r.  One  consequence  of  this  large  drift  near 

resonance  is  that  it  limits  the  degree  of  gain  enhancement  achievable  through 
the  magneto-resonance  effect. 

If  Bq  sufficiently  exceeds  Br,  the  drift  can  be  quite  small.  For 

Y  ■  2.2  and  *  3  cm,  the  resonant  magnetic  field  B^  *  7  kG.  With  Bq  »  8 

kG  and  Bw  *  1  kG,  the  particle  still  drifts  to  the  wall  as  shown  in  Fig.  4a. 

The  drift  is  now  in  the  opposite  direction  to  that  below  resonance 

since  v^  changes  sign  when  Bz  >  Br-  This  has  the  effect  of  changing  the  phase 

of  v  oscillations  with  resoect  to  those  of  3 i.e.  ,  v  is  positive  (for  x  > 
y  *  z  y 

0)  when  S2  is  a  minimum,  so  that  the  particle  drifts  in  the  +y  direction. 

When  BQ  is  increased  to  10  kG  for  this  case,  the  drift  is  small  enough  that 
the  electron  remains  confined  for  >  30  periods  as  shown  in  Fig.  (4b).  The 
confinement  remains  very  good  when  the  self  fields  of  a  500  A  beam  are  added. 

In  principle,  operation  of  FEL  experiments  with  BQ  >  Bf  is  possible  and 
has  been  demonstrated  with  a  high  current,  y  ■  3-5  beam  in  a  helical 
wiggler.4  However,  competing  processes  such  as  the  cyclotron  maser^ 
Interaction  can  produce  large  radiated  powers  at  frequencies  close  to  those  of 
Che  FEL  interaction  in  this  beam  energy  and  magnetic  field  regime. 


in  linear  wiggler  with  I  ■  0,  Y  ■  2.2 
.  a)  BQ  -  8  kG,  b)  Bq  -  10  kG. 
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Consequently ,  analysis  of  experimental  results  is  more  difficult.  Also, 

arbitrarily  large  guide  fields  are  not  possible  with  a  magnetically  shielded 

diode  simply  because  the  electrons  will  be  mirrored  by  the  field.  The 

equilibrium  beam  radius*''  in  such  a  case  depends  only  on  y,  I,  3Q,  and  the 

beam  emittance  e,  so  for  particular  beam  parameters  suitable  values  of  3Q  are 

1 8 

limited.  The  required  field  in  kG  for  a  matched  beam  radius  R  in  cm  can  be 
written 


3 


2 


1.36  I  + 
R23y 


11.56 


(10) 


where  I  is  in  kA  and  is  the  normalized  emittance  in  rad  cm.  For  example, 
if  ■  .14  ir  rad-cm  (about  the  lowest  value  expected  for  a  thermionic  cathode 
beam  with  I  »  750-1000  A)1’®,  ay-  2.2  beam  with  R  »  .3  (.5)  cm  requires 
3Q  »  5.5  (2)  kG  for  I  <  800  A.  In  this  regime,  the  beam  is  emittance 
dominated  so  that  the  required  field  is  relatively  insensitive  to  I.  In 
principle,  smaller  radius  beams  could  be  used  with  larger  B0,  thereby  doubly 
reducing  v^.  However,  experimentally  this  is  very  difficult  at  high  current 
levels. 

To  analyze  the  effect  of  a  shielded  diode  on  propagation  through  the 
wiggler,  we  repeat  the  above  calculations  with  an  initial  va  corresponding 
to  Pg0"  0.  Note  that  I  -  1.75  kA  in  this  case,  which  is  the  current  required 
for  constant  radius  propagation  in  only  the  solenoidal  field  with  these 
initial  conditions.  As  3hown  in  Fig.  5,  the  electron  propagates  at  nearly 
constant  radius  until  the  wiggler  amplitude  becomes  large  enough  that  the 
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gradient  drift  begins  to  dominate.  Then  the  particle  is  lost  just  as  in 


the  VgQ  -  0  case. 

The  behavior  of  an  off-axis  electron  in  a  linear  wiggler  and  axial  guide 

1  *2 

field  should  be  compared  to  that  when  a  helical  wiggler  is  used.  Me 
approximate  the  helical  wiggler  field  by^ 

Bh  *  b( z) B  {-sin  kz  [l  +  4  k2(3x“  +  y2)  1  +  4  x  y  cos  kz} 

(ID 

»  b(z)Bw(cos  kz  [l  +  -g-  k"(x2  +  3y2)]  -  j  k2  x  y  sin  kz} 

Bh  -  -kb( z) B  [ 1  +  4  k2(x2  +  y2) ] fx  cos  kz  +  y  sin  kz] . 

This  expansion  of  the  true  Bessel  function  expression  for  Bh  is  valid 

for  kr  <  1.  Friedland's  treatment1^  of  electron  propagation  in  this  case 

assunes  a  radially  uniform  wiggler,  and  his  stability  condition  (Eq.  2a)  is 

not  stringent  enough  when  the  radial  variation  is  included.  For  example, 

Friedland  finds  stable  orbits  when  y  ■  1.587,  A  *  4  cm,  v,  *  0,  r  »  0, 

w  io  o 

B0  ■  1.26  kG,  and  Bw  »  1.04  kG,  so  that  Eq.  2a  is  barely  satisfied,  and  we  can 
duplicate  his  results  if  we  remove  the  radial  variation  from  311.  However, 
with  the  B^  given  in  Eq .  11,  we  find  that  the  wiggler  field  must  be  reduced  to 
~625  G  to  obtain  stable  orbits. 

Although  the  radial  dependence  in  Bh  does  narrow  the  allowable  range  of 
operating  parameters,  stable  orbits  with  <r>  -  cons' t  are  achievable  with  a 

helical  wiggler  in  a  guide  field.  Electron  trajectories  in  a  helical  wiggler 
for  the  same  conditions  as  used  previously  with  a  linear  wiggler  are  shown  in 
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Figs.  6  and  7.  Figure  6  is  to  be  compared  to  Fig.  2  and  Fig.  7  to  Fig.  5.  In 
both  cases,  the  electron  is  well-conf ined.  It  is  interesting  to  note  that  in 
Fig.  6,  the  electron  born  at  (x,y)  »  (.4,0)  initially  E  x  3  drifts  in  the  + 

9  direction,  but  then  reverses  direction  as  the  wiggler  amplitude 
increases.  Since  grad  3^  is  radial,  the  grad  3^  drift  is  in  the  -  9 
direction  and  does  not  lead  to  the  beam  expansion  observed  with  the  linear 
wiggler. 

The  linear  wiggler  drift  described  here  imposes  additional  constraints  on 
the  parameters  of  an  FEL  experiment.  Obviously,  the  beam  radius  must  be  kept 
as  small  as  possible  to  minimize  particle  loss  from  the  edge  of  the  beam. 
Preliminary  experimental  results  by  our  group  with  a  field  immersed,  apertured 
source  indicate  that  particle  losses  can  be  kept  acceptably  small  in  this 
way.  Also,  if  y  is  large  enough  that  Br  >>  BQ  >>  Bw  can  be  satisfied  for 
relatively  large  Bw,  then  the  drift  can  be  kept  small  while  achieving 
acceptably  large  vw. 

Finally,  it  should  be  noted  that  the  drift  arises  from  the  asymmetry  of 
the  linear  wiggler  and  the  corresponding  absence  of  focusing  forces  in  the 
direction  perpendicular  to  the  wiggler  field.  Therefore,  it  should  be 
possible  to  stabilize  the  drift  by  imposing  an  additional  focusing  force  in 
that  direction.  For  example,  preliminary  results  indicate  that  electron 
propagation  through  a  "square"  or  symmetrized  linear  wiggler  is  very  stable. 

I 

Such  a  wiggler  has  an  additional  component  B^.  ■  cosh  ky  cos  kz  and  a 

l 

corresponding  addition  to  B  of  -3„  sinh  ky  sin  kz.  An  electron  trajectory 
through  such  a  wiggler  is  shown  in  Fig.  8. 
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30 

I  -  1.75  kA,  Y  *  2.2,  3Q  -  2  k G,  3W  -  1  kG,  and  rQ  -  .6  cm 
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In  conclusion,  a  gradient  drift  has  been  shewn  to  exist  for  a  linear 


wiggler  in  an  axial  guide  field.  The  drift  can  be  substantial  with  snail  o 
large  beam  current  in  some  parameter  ranges,  for  a  wide  range  of  initial 
conditions.  However,  the  advantages  of  a  linear  wiggler  are  sufficient  in 
many  cases  to  either  limit  operation  to  a  "stable"  parameter  regime  or  to 
impose  additional  focusing  forces  to  stabilize  the  drift. 

The  auchors  gratefully  acknowledge  useful  discussions  with  C.A. 
Kapetanakos,  H.  Freund,  and  C.M.  Tang. 
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NUMERICAL  SIMULATION  OF  HIGH  ENERGY  LIGHT 
ION  BEAM  FLOW  THROUGH  A  FOILLESS  DIODE 

I.  INTRODUCTION 

Interest  remains  high  in  the  quest  for  a  compact,  high  energy,  intense  light-ion  beam  accelerator. 
Some  of  the  areas  of  application  for  such  an  accelerator  include  directed  energy  systems,  inertial 
confinement  fusion  drivers  and  magnetic  fusion  fuel  injectors.  The  desired  operating  parameters  for 
such  an  accelerator  are  an  ion  energy  of  from  0.1  to  1  GeV,  a  beam  intensity  of  from  1013  to  1016 
ions/ cm2  per  puise.  and  a  pulse  length  of  from  less  than  one  nanosecond  to  a  millisecond.  To  be 
termed  "compact.”  the  device  should  not  measure  more  than  10  to  30  meters  in  any  dimension. 

The  ability  to  generate  intense  relativistic  electron  beams  (IREB’s)  in  foilless  diodes1  2  has 
brought  the  concept  of  a  compact  ion  accelerator  very  close  to  a  reality  in  the  form  of  collective  ion 
accelerators  (CIA's).  These  collective  accelerators  have  received  much  attention  over  the  years.3-6 
Their  operation  relies  upon  the  existence  of  localized,  collective  field  gradients  in  excess  of  1  MV/cm 
in  an  IREB.  In  comparison,  conventional  accelerators  provide  electric  fields  on  the  order  of  only  0.01 
MV/cm  which  is  unacceptable  for  a  "compact"  device.  Although  many  approaches  exist  for  achieving 
the  acceleration  of  positive  ions  via  the  collective  electric  field  produced  by  an  intense  relativistic  elec¬ 
tron  beam,  not  all  approaches  can  be  scaled-up  to  achieve  high  ion  energies.  Scaleabilitv  implies  a  sys¬ 
tematic  procedure  for  maintaining  synchronization  between  the  ion  being  accelerated  and  its  accelerat¬ 
ing  electrostatic  potential  well.  This  requirement  of  scaleability  may  be  generalized  to  any  type  of 
accelerator. 

The  "space  charge  wave  accelerator"  concept  has  been  proposed  at  NRL  to  achieve  such  scaleabil¬ 
ity.7  The  fundamental  idea  behind  the  concept  is  the  sustained  matching  of  ion  beam  velocity  with  the 
phase  velocity  of  the  accelerating  IREB  space  charge  wave  along  the  entire  length  of  the  CIA  (collective 
ion  accelerator)  The  most  formidable  challenge  to  this  concept  is  the  matching  of  injected  ion  beam 
velocity  with  the  minimum  electron  space  charge  wave  velocity  Technical  restraints  appear  to  place  a 
lower  limit  of  0. 2-0.3  c  on  the  phase  velocity  of  the  fundamental  mode  of  a  1  MeV  IREB  space  charge 
wave.  This  dictates  an  injected  proton  energy  of  16-30  MeV.  If  that  synchronization  of  miual  veloci¬ 
ties  can  be  met,  the  geometry  of  the  accelerator  cavity  can  be  designed  to  keep  v  loo  and  v  P(UK  in  step. 
The  specific  design  used  is  termed  the  Converging  Guide  Accelerator  or  CGA.8  The  phase  velocity  of 
the  wave  in  a  conducting  cylinder  depends  on  the  ratio  of  beam  to  wall  radii.  The  phase  velocity 
increases  inversely  with  this  ratio.  Thus,  the  wave  can  be  accelerated  by  accelerating  the  beam.  The 
beam  can  be  accelerated  by  passing  it  through  a  metal  tube  with  a  converging  radius. 

In  Fig.  1,  ions  are  injected  from  the  right  along  the  central  axis  of  the  electron  beam  which  is  also 
propagating  from  right  to  left.  Next,  a  space  charge  wave  is  excited  and  grows  around  the  ions.  The 
ions  loaded  into  the  wave  are  accelerated  by  it  as  it  too  accelerates  in  step.  The  electron  beam  is  then 
dumped  and  the  ions  extracted  from  the  far  left.  Two  meters  of  magnetic  field  at  a  strength  of  20  kG 
are  required  for  this  system.  Initially,  the  electron  beam  planned  for  use  with  this  system  will  have  an 
energy  of  0.5  MeV  and  a  current  of  1.4  kA  for  250  nsec.  Also,  the  electron  beam  generator  can  be 
operated  at  10  pulses  per  second.  The  trace  that  is  shown  in  Fig.  I  is  a  photograph  of  the  machine  vol¬ 
tage  versus  time.  The  high  voltage  machine  was  used  in  conjunction  with  a  resistive  load.  For  proof- 
of-principle  testing,  low-current.  30  MeV  proton  beams  could  be  injected  from  a  conventional  ion 
accelerator  such  as  the  cyclotron  at  NRL. 
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Such  a  low  ion  current  experiment  can  indeed  test  the  wave  accelerator  concept  but  there  remains 
another  question  to  be  answered  before  kiloamp-level,  infected  beam  devices  are  constructed.  What  are 
the  consequences  of  injecting  into  a  foilless  diode  a  30  MeV  proton  beam  of  comparable  amperage  to 
the  current  normally  carried  by  the  electrons  alone  ;n  the  diode.  How  will  the  structure  of  the  potential 
well  inside  the  hollow  electron  beam  be  effected0  Will  the  net  electron  current  flowing  through  the 
diode  increase  and.  if  so.  by  how  much0  These  are  the  matters  addressed  b>  the  computational 
research  reported  herein.  The  findings  demonstrate  that  an  injected  30  MeV  proton  beam  carrying  10% 
of  the  normal  diode  electron  current  will  not  seriously  effect  the  potential  well  structure  m  or  the  elec¬ 
tron  current  of  the  foilless  diode.  A  five- fold  increase  in  the  ion  beam  current,  on  the  other  hand, 
destroys  the  well  and  doubles  the  electron  current.  Thus,  clearly  "safe"  and  clearly  ’unsafe"  ion  current 
levels  have  been  established.  The  determination  of  an  exact  "critical"  ion  current  for  use  in  the  CGA 
device  is  not  important  at  this  time  and  has  been  left  for  future  research. 

In  the  section  that  follows,  the  actual  physical  device  to  be  simulated  is  described  in  detail.  Sec¬ 
tion  III  then  outlines  the  key  characteristics  of  the  D10DE2D  computer  code  which  was  used  to  con¬ 
duct  the  numerical  simulation  The  numerical  results  are  presented  in  Section  IV  including  such  infor¬ 
mation  as  radial  profiles  of  beam  current  density  and  plots  of  the  electrostatic  potential  well  structure  in 
the  diode.  In  order  to  give  added  emphasis  and  understanding  to  the  numerics.  Section  V  is  devoted  to 
a  simple  theoretical  analysis  of  the  physics  of  the  ion  beam  and  diode.  Good  agreement  is  found 
between  the  theoretical  predictions  and  the  numerical  findings.  Finally.  Section  VI  summarizes  the 
major  conclusions  that  can  be  drawn  from  this  work  and  suggests  material  for  future  research. 

II.  THE  FOILLESS  DIODE 

The  physical  diode  designed  for  the  experimental  testing  of  this  new  CGA  concept  closely  resem¬ 
bles  that  shown  in  Fig.  2.  It  corresponds  roughly  to  the  "wave  growth"  section  of  the  accelerator  shown 
in  Fig.  1.  The  coaxial,  vacuum,  power  feed  line  tapers  inward  radially  as  it  traverses  the  axial  length  of 
the  near-conical  cathode.  The  tip  of  the  cathode  coincides  with  the  entrance  plane  of  a  long  (>1 
meter)  cylindrical  drift  tube.  The  outer  radius  of  the  cathode  is  about  1.5  cm.  A  hole  0.5  cm  in  radius 
is  bored  through  the  center  of  the  cathode  to  serve  as  the  injection  port  for  the  high  energy  ion  beam. 
The  coaxial,  conical  A-K  gap  has  a  length  of  about  five  centimeters.  The  drift  tube  is  given  an  inner 
radius  of  about  1.7  cm.  The  entire  apparatus  is  immersed  in  a  uniform  axial  magnetic  field  of  about  10 
kilogauss. 

Ideally,  the  model  diode  to  be  numerically  simulated  would  correspond  exactly  to  that  depicted  in 
Fig.  2.  Unfortunately,  the  simulation  code,  DIODE2D,  which  was  used  in  these  studies  was  incapable 
of  efficiently  handling  that  specific  geometry.  Particularly  troublesome  are  the  sloped  walls  of  the 
cathode  and  of  the  coaxial  power  feed  line  as  well  as  the  overall  meter-long  length  of  the  physical  sys¬ 
tem.  The  existing  code  was  capable  of  treating  only  conducting  surfaces  and  boundaries  which  are 
strictly  radial  or  strictly  coaxial  (i.e.,  they  must  have  a  rectangular  R-Z cross-section) .  Thus,  a  cylindr¬ 
ical  cathode  shank  of  some  finite  thickness  must  be  substituted  for  the  conical  one.  Furthermore,  a 
combination  of  limited  resolution  and  finite  computer  data  storage  capabilities  prohibits  the  accurate 
treatment  of  a  system  less  than  2  cm  in  radius  but  over  100  cm  long.  The  total  axial  length  of  the 
model  system  will  be  shortened  instead  to  10  cm.  The  peculiarities  of  the  Poisson-solver  for  the  elec¬ 
tric  and  magnetic  potentials  in  that  system  also  required  that  the  tube  be  terminated  by  a  fixed- 
potential.  conducting  wall.9  The  final  approximation  to  the  original  svstem  took  the  form  shown  in  Fig. 
3. 

The  cathode  used  in  the  numerical  model  is  a  hollow  cylinder  1.2  cm  in  radius  with  a  shank  wall 
thickness  of  0.4  cm.  It  projects  axially  2.5  cm  into  the  10.0  cm  long  drift-tube  (anode).  The  anode 
tube  inner  radius  is  at  1.7  cm.  This  leaves  a  radial  A-K  (anode-cathode)  gap  of  0.5  cm  which  very 
comfortably  exceeds  the_  0.17  cm  gyroradius  of  a  0.5  MeV  electron  in  a  10  kG  magnetic  field.  Analysis 
using  existing  theory1041  leads  to  estimates  of  foilless  diode  electron  beam  currents  of  8-10  kA  or  diode 


power  levels  of  4-5  GW  Along  the  central  axis  of  this  device  will  be  injected  a  0.7  cm  radius  beam  of 
30  MeV  protons.  Three  separate  simulations  were  conducted  of  this  configuration.  The  first  was  a 
"benchmark"  run  to  determine  the  normal,  eiectrons-only  operating  characteristics  of  the  foilless  diode. 
Of  importance  is  not  only  the  gross  e-beam  current  but  also  the  beam  profile  and  the  structure  of  the 
electrostatic  potential  well  inside  the  beam.  The  final  two  simulations  were  conducted  with  the  high- 
energy  proton  beam  flowing  through  the  device.  Specific  ion  beam  currents  of  1  kA  and  5  kA  were 
chosen  for  these  two  test  cases. 

III.  NUMERICAL  IMPLEMENTATION 

A  modified,  electrostatic-magnetostatic  version  of  NRL's  D10DE2D  computer  code  was  used  to 
carry  out  the  simulations.  The  details  of  the  code  may  be  found  elsewhere. 1J  It  is  sufficient  here  to 
note  that  the  code  is  2-1/2-dimensional  in  that  it  explicitly  moves  particles  in  R  and  Z  'radially  and  axi¬ 
ally)  in  space  while  keeping  track  of  all  three  components  of  their  canonical  momentum.  The  treatment 
is  fully  relativistic  for  both  electrons  and  protons.  Unlike  a  full  electromagnetic  simulation  code. 
D10DE2D  is  designed  to  find  only  equilibrium  charge,  current,  and  field  configurations  in  the  device 
under  study. 

Particle  contributions  to  net  charges  and  currents  as  well  as  the  resulting  E  and  B  field  values  are 
stored  in  a  mesh  of  discrete  data  points  covering  the  entire  diode.  The  spacing  between  data  points  was 
fixed  at  R  —  0.05  cm  radially  and  Z  “  0.078  cm  axially.  This  allowed  for  34  data  cells  to  span  the  l.7 
cm  from  the  centerline  to  the  inner  radius  of  the  anode  tube  and  128  cells  from  the  back  of  the  cathode 
to  the  anode  endplate.  The  protruding  cathode  shank  measured  8  data  cells  thick  and  32  cells  long.  On 
its  front  face,  electron  emission  was  permitted  from  the  outer  7  of  the  8  possible  emission  points.  The 
eighth  cell  was  not  allowed  to  emit  in  order  to  obtain  additional  numerical  accuracy  in  the  field  values 
governing  electron  emission  from  the  inside  lip  of  the  cathode  shank.  Emission  was  permitted  all  along 
that  inner,  as  well  as  the  outer  shank  surface  except  for  the  rearmost  6  cells. 

The  magnitude  of  the  electron  emission  from  a  given  cathode  surface  is  determined  by  the 
amount  of  space  charge  necessary  to  zero  out  the  perpendicular  component  of  the  electric  field  at  the 
emission  point.  In  order  to  simulate  the  injection  of  the  30  MeV  proton  beam  of  predetermined 
amperage,  the  inner  14  cells  of  the  cathode  endplate  inside  the  shank  were  designated  as  proton  "emis¬ 
sion"  points.  The  emission  at  those  points,  however,  was  completely  independent  of  the  adjacent  elec¬ 
tric  field  values.  A  constant,  emitted  current  density.  J,.  of  0.65  kA/cm3  was  imposed  for  the  1  kA 
beam  case  and  3.25  kA/cm3  for  the  5  kA  case  at  each  of  the  14  points.  Particles  of  both  species  were 
perfectly  absorbed  upon  hitting  any  of  the  boundary  surfaces,  including  the  surfaces  of  the  protruding 
cathode  shank. 

All  quantities  of  interest  in  the  device  were  monitored  by  the  computer  code’s  extensive  diagnos¬ 
tics.  Electron  and  ion.  emitted  and  collected  current  densities  were  recorded  over  all  relevant  surfaces. 
Electrostatic  equipotential  contour  lines  and  magnetic  field  lines  throughout  the  device  were  plotted. 
Sample  values  of  all  three  components  of  the  magnetic  field  are  similarly  recorded.  Net  charges  and 
currents  were  periodically  listed.  Finally,  the  positions  of  statistical  samplings  of  electrons  and  protons 
were  plotted  at  equilibrium. 

IV.  RESULTS 

First,  a  "benchmark"  simulation  was  run  in  which  only  electrons  were  allowed  in  the  diode.  In 
agreement  with  theory  (see  Section  V),  a  net  electron  current  of  9.6  kA  was  measured.  Of  that  total 
amount,  2.8  kA  was  emitted  from  the  face  of  the  cathode  shank  tip  while  the  remaining  6.8  kA  ori¬ 
ginated  from  the  outside  surface  of  the  shank.  Emission  from  the  inner  shank  surface  was  negligible. 
A  net  electron  charge  of -1.86  x  104  staicouiombs  was  found  to  be  in  the  system.  Figure  4  presents  a 
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plot  of  sample  electron  positions  at  equilibrium.  Of  note  is  the  waviness  of  the  outer  envelope  of  the 
electron  beam.  Similar  beam  behavior  was  observed  in  foilless  diode  simulations  conducted  by  R  Jack- 
son14  15  and  others.  The  electrostatic  potential  contour  plot  for  the  device  modeled  in  simulation  is 
shown  in  Fig.  5.  The  beam  envelope  waviness  is  mirrored  in  these  equipotentials.  It  is  also  clear  from 
this  figure  that  the  electric  field  inside  the  hollow  cathode  is  negligible,  in  agreement  with  the  vanish¬ 
ingly  small  electron  emission  observed  there.  It  also  appears  that  electrons  leaving  the  cathode  are 
fairly  rapidly  accelerated  to  the  40%  equipotentiai  <0.2  Me  Vi  before  entering  a  virtual  "drift  region” 
about  3.5  cm  long  starting  about  2.0  cm  downstream  of  the  cathode  tip.  Electron  energy  and  velocity 
are  nearly  constant  in  this  region  as  would  be  the  case  in  an  actual  full-length  foilless  diode  drift  tube. 
It  is  therefore  that  region  where  beam  behavior  will  most  closely  mimic  a  realistic  device.  Finally,  the 
self-magnetic  field  lines  generated  by  the  e-beam  are  depicted  in  Fig.  6.  Imposed  over  this  self-field 
pattern  is  the  uniform  10  kG  background  field.  Up  to  about  2  “  8.5  cm.  the  beam  field  is  diamagnetic 
in  nature,  lowering  the  inside,  imposed  Bz  by  about  40-80  gauss  and  intensifying  the  field  near  the 
anode  tube  by  a  like  amount.  The  self-field  between  Z  “  8.5  and  10.0  cm.  where  electrons  receive  a 
final  0.3  MV  of  acceleration,  is  of  the  opposite  orientation.  However,  that  region  as  well  can  boast  of 
only  50-100  gauss  field  strengths. 

The  second  case  to  be  tested  was  that  involving  the  1  kA,  30  MeV  proton  beam  injected  in  the 
manner  described  in  the  previous  section  The  net  electron  current  in  the  new  equilibrium  increased  by 
2  kA  to  11.6  kA.  Of  that  total  current.  4.3  kA  was  being  emitted  from  the  tip  (compared  to  2.8  kA  for 
the  electrons-only  case),  0.4  kA  were  coming  from  the  inner  shank  surface,  while  the  outer  shank  sur¬ 
face  emission  remained  practically  unchanged  at  6.9  kA.  The  ion  beam  had  introduced  a  total  charge  of 
1-4  06  x  103  statcoulombs  into  the  diode  while  the  electron  charge  had  increased  to  -2.11  x  ]04 
statcoulombs.  Thus,  the  electron  beam  characteristics  of  the  diode  itself  had  changed  very  little.  In 
addition,  the  ion  beam  envelope  experienced  virtually  no  ballooning.  This  is  evident  from  Fig.  7  which 
presents  the  steady-state  ion  beam  envelope  superimposed  over  the  sample  electron  position  plot.  The 
qualitative  character  of  the  hollow  electron  beam  has  likewise  remained  essentially  unchanged.  The 
first  truly  major  change  manifests  itself  in  the  equipotentiai  plot  of  Fig.  8.  Comparison  with  Fig.  5 
shows  the  distortions  caused  by  the  concentration  of  positive  space  charge  along  the  axis.  The  contours 
are  sucked  into  the  hollow  cathode  up  to  the  40%  potential.  This  accounts  for  the  turn-on  of  electron 
emission  there.  Also,  electrons  reaching  the  previously-mentioned  "drift  region"  between  Z  -  4.5  to  Z 
“  8.0  cm  are  now  at  aboui  0.25  MeV  instead  of  the  0.20  MeV  of  the  electrons-only  case.  Finally.  Fig. 
9  shows  the  changes  made  to  the  self-magnetic  field  lines  generated  by  the  electron  beam.  The  magni¬ 
tudes  of  those  fields  in  the  drift  region  remain  unchanged  in  the  range  of  40  to  80  gauss,  but  the  orien¬ 
tation  of  the  field  has  reversed.  The  imposed  Bz  is  now  strengthened  inside  the  beam  and  weakened 
outside.  Although  this  effect  is  interesting,  the  small  field  strengths  make  it  unimportant  to  the  opera¬ 
tion  of  the  diode. 

For  the  third  and  final  test  case,  the  current  of  the  30  MeV  proton  beam  was  increased  to  5  kA. 
The  steady  state  electron  current  rose  to  19.7  kA  with  8.3  kA  from  the  shank  face.  4.7  kA  from  its 
inner  surface,  and  an  outer  shank  emission  current  almost  identical  to  the  previous  two  cases  at  6.7  kA. 
Of  course,  the  ion  charge  increased  by  the  same  factor  of  five  as  its  current  to  ■*■2.04  x  ]04 
statcoulombs.  In  partial  compensation,  the  net  electron  charge  had  grown  to  -3.32  x  10*  statcoulombs. 
Clearly  the  normal  operation  of  the  foilless  diode  has  been  seriously  disrupted.  Its  electron  current  has 
more  than  doubled.  Even  more  dramatic  is  the  destruction  of  the  ion  beam.  This  is  illustrated  by  the 
equilibrium  sample  particle  plot  of  Fig.  10.  The  proton  beam  envelope  has  expanded  halfway  through 
the  width  of  the  hollow  electron  beam  by  the  time  it  arrived  at  the  end  of  the  drift  tube.  Similarly  the 
electrostatic  potential  contours  shown  in  Fig.  1 1  no  longer  bear  much  resemblance  to  the  ion-free  case 
of  Fig.  5.  Electrons  now  experience  almost  80%  of  the  total  diode  potential  increase  before  entering 
the  "drift  region."  There  no  longer  exists  any  effective  "potential  well"  through  which  the  proton  beam 
can  travel.  Finally,  the  self-magnetic  field  of  the  electrons  (see  Fig.  12)  can  no  longer  be  ignored.  Its 
drift  region  magnitudes  now  range  from  about  200  to  700  gauss  which  is  an  appreciable  fraction  of  the 
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imposed  10  kG  field.  As  in  the  1  kA  ion  beam  case,  this  self-field  strengthens  the  Bz  along  the  axis 
and  weakens  it  outside  the  e-beam  radius. 

The  major  results  of  the  above  three  simulations  are  summarized  in  Table  1.  The  statement  can 
be  made  that  the  injection  of  a  1  kA.  30  MeV  proton  beam  along  the  axis  of  the  foilless  diode  being 
studied,  does  not  seriously  disrupt  the  normal  operating  characteristics  of  the  device.  The  injection  of  a 
5  kA,  30  MeV  ion  beam  does  cause  serious  disruption.  Correspondingly,  the  1  kA  beam  is  virtually 
uneffected  by  its  passage  through  the  diode,  while  the  5  kA  beam  undergoes  severe  ballooning.  In 
order  to  better  understand  these  macroscopic  results,  it  is  important  to  examine  the  modifications  made 
by  each  case  on  the  individual  current  density  profiles  of  each  case  Figure  13  depicts  the  radial  profiles 
of  the  electron  current  density.  Je.  emitted  from  the  up  face  of  the  cathode  shank.  The  profiles  all 
have  a  characteristic  double-pronged  shape.14  They  show  that  for  this  diode  the  doubling  of  the  emitted 
currents  on  this  surface  is  predominantly  occurring  at  the  points  of  lower  radius.  Those  points  are 
closer  to  the  ion  beam,  whose  space-charge  has  grossly  distorted  the  electric  field  values  there.  From 
Fig.  14.  however,  it  is  obvious  that  the  electron  emission  from  both  the  inner  and  the  outer  shank  sur¬ 
faces  are  virtually  identical  for  the  eiectrons-only  and  the  1  kA  beam  case.  Since  some  70%  of  the  total 
diode  current  originates  on  the  outer  surface,  this  agrees  with  the  overall  close  similarity  between  the 
two  cases.  The  5  kA  beam  case  is  quite  a  different  matter.  Although  the  net  outside  shank  emission  is 
about  the  same  as  the  previous  cases,  the  axial  profile  of  Jt  is  grossly  depressed  over  the  front  half  of 
the  shank  and  mildv  enhanced  over  the  rear  half  This  is  a  manifestation  of  the  2  1 /2-fold  increase  in 
the  value  of  B„  along  that  surface.  The  near-normal  amount  of  negative  charge  emitted  from  the  rear 
half  is  prevented  from  leaving  the  near  vicinity  of  the  shank  as  it  travels  toward  the  front  face.  This 
negative  space  charge  then  impedes  the  electron  emission  along  the  front  half  of  that  surface.  The 
most  dramatic  change,  however,  is  in  the  inner  shank  surface  emission.  For  the  5  kA  beam  case,  the 
emission  there  rivals  that  of  the  outside  surface.  For  the  other  two  cases  that  inner  surfaces  had  been 
effectively  non-emming.  The  combined  emission  changes  are  all  reflected  in  the  radial  profiles  of  elec¬ 
tron  current  densities  collected  at  the  anode  face  (Z  “  10  cm)  These  are  shown  in  Fig.  15.  In  addi¬ 
tion  to  the  overall  enhancement  of  collected  current  from  case  to  case,  a  distinct  pinching  of  the  beam 
toward  the  central  axis  can  be  seen.  This  is  to  be  expected  due  to  the  increased  self-magnetic  field 
strength  but  the  magnitude  of  the  radial  displacement  of  the  mean  current  density  is  limited  by  the 
imposed  10  kG  axial  5-field. 

The  final  diagnostic  measured  the  radial  profiles  of  the  30  MeV  proton  beam  in  the  plane  located 
at  Z  ”  7.5  cm.  well  inside  the  electron  beam  "drift  region."  Two  profiles  are  shown  for  the  1  kA  beam 
case  in  Fig.  16.  The  first  represents  the  initial,  imposed  beam  profile  a  mere  50  timesteps 
(At  —  2  x  10~a  sec)  after  the  beginning  of  the  simulation.  The  other  profile  shows  the  modest  spread 
of  the  beam  about  0.8  nano-seconds  later.  The  change  is  quite  negligible,  indicating  excellent 
confinement  and  stability  of  the  beam.  Once  again,  the  5  kA  case  is  dramatically  different.  The 
corresponding  profiles  for  that  high  current  beam  are  shown  in  Fig.  17.  Over  the  same  0.8  nanosecond 
period,  this  beam  has  expanded  from  its  initial  0.7  cm  radius  to  over  1.0  cm.  well  beyond  the  inner 
cathode  shank  radius  and  the  inner  radius  of  the  hollow  electron  beam.  That  configuration  is  clearly 
unstable. 

V.  THEORY  COMPARISON 

t 

In  this  section  analytic  expressions  are  derived  and  compared  to  the  simulation  results.  The  two 
major  results  are  the  radial  ballooning  of  the  ion  beam  and  the  increase  in  the  electron  limit  current  as 
the  ion  current  is  increased.  Agreement  between  theory  and  computation  is  quite  good. 

The  radial  equation  of  motion  for  the  ions  decouples  and  only  the  ion  beam's  self  fields  determine 
the  ion  motion.  This  is  true  for  times  that  are  short  in  comparison  to  the  axial  gyro-period.  That  is. 

r1  <  <  nr:,  where.  Cl.  -  — -  and  B.  is  the  axial  magnetic  field.  Also,  for  <  <  1/2.  where  J3  is 
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the  ion  velocity  divided  by  the  speed  of  light,  then  changes  in  the  relativistic  mass  factor  can  be 
neglected.  With  these  assumptions  the  equation  of  motion  for  the  tons  at  the  beam  surface  is  given  by 


q£,  -  qV-Bf, 
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where;  q  and  A/  are  the  ion  charge  and  mass.  £r  and  Be  are  the  seif  electric  and  magnetic  fields  in  the 
ion  beam,  y  is  the  relativistic  mass  factor.  V.  is  the  axial  ion  velocity,  and  r  is  the  radial  position  of  ions 
at  the  beam  surface.  Assuming  current  conservation  Eq.  (1)  reduces  to1' 
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Equation  (2)  can  be  solved  in  series  form18  to  give. 


is  the  dielectric  constant  of  free  space 
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where;  r0  is  the  initial  beam  radius  at  t  —  0  and  r  is  the  final  radius  at  time  t  In  Table  2,  a  comparison 
is  made  between  the  code  time,  u  to  reach  the  final  radius,  r,  and  the  time  calculated  by  Eq.  (3).  This 
comparison  is  made  for  30  MeV  protons. 


It  is  clear  from  the  comparison  that  the  self  fields  of  the  ion  beam  determine  its  motion.  Of  the 
self  fields,  the  electric-field  dominates  in  the  expansion.  This  result  does  not  limit  the  maximum  ion 
current.  Since,  the  ion  beam  expands  to  the  diameter  of  the  electron  beam,  it's  electric  field  is  reduced 
thus  allowing  for  more  ion  current  to  flow.  Increasing  the  ion  current  then  increases  the  electron  iimii 
current  which  in  turn  allows  for  an  increase  in  the  ion  current,  etc. 


In  addition  to  the  problem  of  ion  beam  ballooning,  the  simulation  also  addressed  the  major  ques¬ 
tion  of  changes  in  the  diode's  electron  current  due  to  the  presence  of  the  ion  beam.  The  maximum 
electron  current  is  determined  by  the  potential  distribution  along  the  electron-emitting  cathode  surfaces. 
When  ions  are  present  the  potential  is  raised  and  hence  the  limit  current  is  increased. 


In  Fig.  18  electrons  and  ions  are  moving  in  the  same  axial  direction.  A  large  external  magnetic 
field  is  applied  in  the  direction  of  particle  motion,  z  This  field  is  large  enough  to  magnetize  the  elec¬ 
trons  but  not  the  ions.  Furthermore,  the  assumptions  are  made  that  (a)  Je  and  J,  are  uniform,  (b) 
l£»,  I  and  \Bt,\  «  Ifi.l,  (c)  l£j  <  l£.ic,  and  (d)  Electron  velocity  shear  can  be  neglected  (where  Je 
and  7,  are  the  electron  and  ion  current  density;  B#,.  Btf  and  B.  are  the  ion.  electron  and  external,  mag¬ 
netic  fields;  Er,  is  the  radial  ion  electric  field;  c  is  the  vacuum  speed  of  light). 


From  Gauss'  law  the  potential  inside  the  electron  beam,  region  3.  (see  Fig.  20)  is  given  by  the 
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where  /,  and  /,  are  the  ion  and  electron  current,  Vi  and  K  are  the  ion  and  electron  velocity,  and  «0  is 
the  permitivity  of  the  free  space.  The  quantity.  tt> ,  may  be  eliminated  from  Eq.  (4)  by  using  conserva¬ 
tion  of  energy. 
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The  maximum  electron  current  can  then  be  found  bv  setting  — —  to  zero,  thus  yielding 
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mc'y o  is  the  initial  electron  energy,  m  and  e  are  the  electronic  mass  and  charge,  and  0  is  the  ion  velo¬ 
city  normalized  to  the  vacuum  speed  of  light. 

By  inspecting  the  first  term  in  Eq.  (8)  it  is  clear  that  the  presence  of  the  ion  current  raises  the 
final  electron  energy.  This  energy  increase  plays  only  a  small  role  in  the  limn  current  expression.  The 
biggest  ion  effect  arises  in  Eq  <4)  where  the  potential  is  directly  increased  by  an  increase  in  the  ion 
current. 

The  magnitude  of  the  limiting  electron  current  of  Eq.  (6)  varies  with  the  radial  position,  r.  In 
general  Eq  (6)  should  be  integrated  from  r,  to  r;  to  give  a  current  that  is  weighted  by  the  changing 
potential  across  the  beam.  However,  since  the  equation  itself  is  only  a  simple  approximation,  seeking 
such  extra  "precision”  would  be  pointless.  Instead.  Eq  <61  will  be  evaluated  at  different  radii  and  the 
effects  discussed-  An  upper  bound  to  the  limit  current  can  be  obtained  by  setting  r  **  r:  then  Eq.  (6> 
becomes 
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This  gives  an  upper  bound  to  the  limit  current  since  all  particles  are  experiencing  the  maximum  poten¬ 
tial  near  the  wall.  A  comparison  is  made  in  Table  3  ff  Eq.  ( 6 )  is  evaluated  at  r  —  r,  then  the  result  is. 
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Note  that  when  no  ions  are  present,  this  will  give  a  lower  bound  to  the  limit  current  since  all  of  the 
particles  are  exposed  to  a  potential  minimum.  Table  4  shows  a  comparison  of  Eq.  UOi  with  the  code 
results.  As  expected,  the  theoretical  prediction  is  wei’  below  the  code  findings  for  the  case  with  no  ion 
beam  present.  For  the  two  cases  with  ion  injection,  however,  the  addition  of  positive  space  charge 
inside  r  —  r,.  elevates  the  potential  there  to  such  a  degree  that  for  the  5  kA  case  it  actually  exceeds 
that  observed  in  Table  3  for  r  •  r2.  It  is  clear  from  both  Table  3  and  Table  4  that  the  comparison 
between  simulation  and  theory  is  quite  good,  considering  the  theoretical  approximations  that  were 
made.  Thus  a  scaling  law  for  the  electron  limit  current  dependence  on  the  ion  current  can  be  inferred. 
This  effect  can  have  a  major  impact  on  wave  phase  velocity  control  for  the  space  charge  wave  a'celera- 
tor.  For  example  if  the  injection  current  is  fixed  a  change  in  the  ion  current  will  shift  the  wave  phase 
velocity  thus  removing  ions  from  wave  synchronization.  However,  at  high  y0  and  for  modest  :on 
currents  the  effect  is  small.  For  example,  if  /,//,  -  0.1  and  yo  ”  TO. 
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Thus,  the  two  major  code  results  have  been  explained.  Analytic  expressions  have  been  derived 
for  predicting  these  results  and  good  agreement  was  found  between  theory  and  code  results.  The 
mechanism  for  ion  ballooning  has  been  shown  to  be  due  primarily  to  the  self  electric  field  of  the  ion 
beam.  Increases  in  the  electron  limit  current  are  shown  to  be  associated  with  an  increase  in  the  ion 
current.  As  discussed,  a  coupling  between  the  ion  and  electron  currents  can  negatively  impact  on  lon- 
wave  synchronization  for  the  space  charge  wave  accelerator 

VI.  CONCLUSIONS 

The  numeric  simulations  presented  here  give  the  first  documented  evidence  that  a  high  current, 
high  energy  light  ion  beam  can  be  injected  down  the  central  axis  of  a  conventional  foilless  diode 
without  either  seriously  disturbing  the  operating  characteristics  of  the  diode  or  causing  significant  disr¬ 
uption  of  the  ion  beam  itself.  For  the  4-5  GW  foilless  diode  studied  and  for  a  proton  energy  of  30 
MeV  a  beam  current  of  1  kA  can  easily  be  tolerated.  Increasing  the  ion  beam  current  to  5  kA  leads  to 
a  general  break-down  of  the  beam-diode  system.  The  precise  current  value  between  1  kA  and  5  kA  at 
which  system  stability  ceases  remains  to  be  found.  At  y  >  7  a  ten  percent  ion  current  will  cause  less 
than  a  one  percent  variation  in  beam  energy.  This  appears  acceptable  for  collective  wave  accelerators. 
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Table  1  —  Diode  Currents  and  Charges 
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Table  2  —  Comparison 
of  Code  and  Analytic  Calculated  Times 
to  Reach  Final  Radius 
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Table  3  —  Comparison  of 
Code  and  Analytic  Limit  Currents 
Potential  evaluated  at  r  —  r,. 
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Table  4  —  Comparison  of 
Code  and  Analytic  Limit  Currents. 
Potential  evaluated  at  r  -  r,. 


/, on  <kA> 

0 

1 

5 

1  electron  (kA)  from 

Eq.  (10) 

7.3 

11.01 

27.85 

/electron  <kA)  from 

Code 

9.6 

11.6 

19.7 

%  Relative  Deviation 
from  Eq.  (10) 

31.5 

5.45 

29 

10 


Wave  Acceleration 


wave  Growth 


I  =  1.4kA 

7  =  250  Rs 
RR  =  lOpps 


V 


100  ns/DIV 


Fig.  1  —  The  converging  guide  accelerator 


S.Ocm 


Fig.  2  —  The  experimental  diode  geometry 


Self-magnetic  fieid  for  diode  with  1  kA. 
30  MeV  proton  beam 


Fig.  11  —  d>(r,  i)  for  diode  with  S  kA, 
30  MeV  proton  beam 


Fig.  12  —  Self-magnetic  field  for  diode  with  5  kA. 
30  MeV  proton  beam 


16 


R(cm ) 

Fig.  13  —  J'U)  emitted  from  cathode  up 


ZICM) 


Fig.  14  -  Profiles  of  electron  shank  emission 


JUA/cni*) 


High  Current,  High  Voltage  Accelerators  as 
Free-Electron  Lasers  Drivers 

C.  W.  Roberson,  J.  A.  Pasour,  C.  A.  Kapetanakos,  Sprangle, 
0.  Golden,  F.  Mako  and  R.  Lucey 
Naval  Research  Laboratory,  Washington,  DC  20375 

Free-Electron  Generators  of  Coherent  Radiation 
June  22-25,  1981 


-186- 


Free-Electron  Generators  of 
Coherent  Radiation 


Based  on  Lectures  of  the  Office  of  \aval  Research,  Sponsored  Workshop. 
June  22-25.  1081.  Sun  Valley.  Idaho 


edited  by 

Stephen  F.  Jacobs,  Optica/  Sc.-mccs  Center.  University  of  Arizona,  Tucson 

Gerald  T.  Moore,  Institute  for  Modern  Optics ,  Department  of  Physics  and 
Astronomy,  U nrversity  of  Ve.v  Mexico 

Hcrschel  S.  Pilloff,  Office  of  Saval  Research,  Arlington,  V:'ginia 

Murray  Sarye.it  III,  Optical  Sciences  Center,  University  of  Arizona,  Tucson 

Marian  O.  Scully,  'institute  for  Modem  Optics.  Department  of  Physics  and 
Astronomy,  University  of  .Vow  Mexico 

Richard  Spit/cr,  Stanford  linear  Accelerator  Center,  Stanford  University, 
Stanford.  California 


A 

77 

1982 


Addison-Wesley  Publishing  Company 
Advanced  Book  Program/World  Science  Division 
Reading,  Massachusetts 

!  C»-.fo n  •  Am>terdjm  •  Don  MuU  OntAtiO  •  $.  .frev  •  *‘AyO 


high  current,  high  voltage  accelerators  as 
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P.  Sprangle,  J.  Golden,  F.  Mako$  and  R.  Lucey* 
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ABSTRACT 


Thera  are  several  approaches  to  generate  high  voltage 
m.ulti-kiloam.pere  electron  bears  under  Investigation  at  the 
Naval  F-ssearch  Laboratory.  Such  accelerators  could  provide 
compact  drivers  for  high  power  Free  Electron  Lasers.  Three 
of  these  potential  FEL  drivers  are  discussed  in  this 
paper.  .hey  are  ! 1)  the  Long  Pulse  Induction  L ; nee,  (2)  the 
Racetrack  Induction  Accelerator  ami  (3)  the  Modified 
Betatron. 

INTRODUCTION 


During  the  past  fifty  years  accelerators  have  developed 
primarily  along  two  diametr ica 1 ly  opposite  directions:  (  1) 
low  current,  high  voltage  and  (2)  high  current,  low  voltage 
devices. 


$  JAYCOR 

t  Pulse  Sciences,  Inc. 

727 


S 


72$ 


C.  ROBERSON  ET  AL 


33 


Figure  1  is  a  plot  of  bean  energy  vs  average  bean 
current  of  typical  accelerators  developed  during  this 

period.  The  history  of  the  traditional  accelerator 

1  2 

development  is  quite  well  known.  '  Marx  generators  with  a 
pulse  forming  line  have  been  developed  primarily  as  flash  x- 
ray  sources  or  inertial  fusion  drivers. 

The  free  electron  laser  requires  high  voltage  for  short 
wavelength  operation  and  high  current  for  efficiency  and 
high  gain  operation.^  Only  the  Induction  Linac  has  been  a 
serious  candidate  as  an  PEL  driver  in  this  parameter  range 
to  date. 
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This  paper  describes  three  accelerator  approaches  being 
pursued  at  NHL  that  could  be  used  as  drivers  for  free 
electron  lasers;  (1)  the  Long  Pulse  Induction  Linac  (2)  the 
Racetrack  Induction  Accelerator  and  (3)  the  Modified 
Betatron. 

THE  LONS  PULSE  INDUCTION  LINAC 

Figure  2  is  a  schematic  of  the  Long  pulse  Induction 
Linac.  Ibis  device  was  built  at  the  National  Bureau  of 
Standards4  and  is  currently  in  operation  at  the  Naval 
Reserach  Laboratory.  'Rie  accelerator  consists  of  two  major 
components;  (1)  an  injector  and  (2)  and  induction 
accelerator  module. 

The  electron  gun  in  the  injector  has  a  16.5  cm  diameter 
tungsten  dispenser  thermionic  cathode.  Electrons  are 
accelerated  in  the  gun  through  a  series  of  22  annular 
electrodes,  spaced  by  ceramic  insulator  rings.  The  last 
electrode  has  a  95%  transmission  screen  of  tungsten  and  is 
at  ground  potential. 

The  electron  gun  is  immersed  in  an  oil  filled  tank.  The 
gun  voltage  is  fed  from  a  pulse  line  driving  a  12:1  step  up 
transformer.  The  injector  typically  produces  a  0.8  kA 
current  beam  pulse  of  400  keV  energy.  The  electron  beam  is 
transported  to  the  induction  accelerator  module  by  a  series 
of  focusing  coils. 

The  induction  accelerator  module  consists  of  two  core 
sets.  One  of  the  core  sets  gives  a  4  to  1  step  up  voltage 
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and  the  other  •  5  to  1.  The  cores  ere  wound  with  0.001  inch 
mild  steel  foil,  separated  by  0.0002S  inch  mylar  sheets. 

Each  accelerating  gap  produces  about  200  kV  of  acceleration. 

Typically,  the  output  energy  of  the  electron  beam 
generated  by  the  linac  is  0.8  MeV,  the  current  approximately 
0.8  kA  and  the  pulse  length  2  microseconds.  The  temporal 
variation  in  the  voltage  is  less  than  3  percent  over  1.6 
microseconds.  This  is  the  longest  pulse  induction  linac  in 
existence.  Shorter  pulse  induction  linacs  of  40  nsec®'®  and 
300  nsec*  have  been  built.  The  pulse  length  becomes  an 
important  consideration  for  free  electron  laser  experiments 
where  one  wishes  to  study  the  nonlinear  dynamics  of  the 
beam,  or  efficiency  enhancement  schemes.  Applications  that 
require  a  significant  amount  of  energy  in  the  radiation 
field  also  require  long  pulses  to  avoid  the  problems 
encountered  with  excessive  electric  field  strengths  at  short 
pusles. 

He  have  designed  and  constructed  a  free  electron  laser 
experiment  that  will  use  as  a  driver  the  beam  from  the 
induction  linac,  which  will  be  focused  to  0.5  -  0.75  cm.  A 
pulsed,  120  cm  long,  3  cm  wavelength  linear  wiggler  will 
provide  the  field  modulation.  The  wiggler  amplitude  rises 
adiabatically  in  30  cm,  has  a  uniform  straight  section  of  60 
cm  and  decays  adiabatically  in  the  last  30  cm.  There  is  a 
uniform  axial  field  over  the  length  of  the  wiggler  that  can 
vary  from  1  to  5  kG.  The  output  radiation  is  expected  to  be 
in  the  3-4  mm  range  and  the  theoretical  efficiency  is  about 
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THE  RACETRACX  INDUCTION  ACCELERATOR 

The  Racetrack  Induction  Accelerator  is  a  new  high 

Q 

current-high  voltage  cyclic  accelerator  concept.  The 
accelerator  is  designed  to  take  advantage  of  long  pulse 
induction  accelerator  modules  to  obtain  a  high  voltage,  high 
current  electron  beam  in  a  geometry  that  is  suitable  for 
free  electron  laser  operation.  The  geometry  is  an  oval 

racetrack  shaped  toroid,  similar  to  a  model  C 

g 

stellerator.  The  geometry  is  shown  in  Figure  3.  In  one 
leg  of  the  racetrack  a  long  pulse  induction  accelerator 
module  is  inserted.  ttie  voltage  gain  of  a  particle  as  it 
goes  around  the  racetrack  is  approximately 


where  V  is  the  voltage  gain  per  pass  through  the  induction 
module,  T  is  the  time  the  induction  accelerator  module  is 
on,  T  is  the  time  it  takes  a  particle  to  go  around  the 
racetrack,  and  is  the  final  voltage.  In  the  first  stage 
of  operation  we  plan  to  construct  a  5  MeV,  1  kA  racetrack 
accelerator.  The  induction  accelerator  module  from  the 
existing  linac  is  capable  of  accelerating  a  beam  to  20  MeV 
in  the  geometry  shown  in  Figure  3. 

The  high  current  operation  of  the  racetrack  requires  a 
magnetic  field  in  the  direction  of  particle  motion.  This 
prevents  the  beam  from  blowing  up  due  to  space  charge. 

The  toroidal  magnetic  field  can  cause  the  beam  to  drift 
out  of  the  accelerator  as  a  result  of  the  field  curvature  in 
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Figure  3 


2  METERS 
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the  bends,  this  drift  can  be  averaged  out  by  applying 
aa  *  2  helical  winding  in  the  bends.  For  moderate 
energies,  the  magnetic  fields  of  such  an  accelerator  can  be 
steady  state.  An  advantage  of  such  a  fixed  field 
accelerator  is  that  the  shape  of  the  accelerating  pulse  is 
not  critical,  the  final  energy  of  accelerated  particles  is 
simply  the  volt-secs  of  the  induction  module,  divided  by  the 
transit  time  and  is  independent  of  pulse  shape,  this  is 
true  provided  the  final  energy  does  not  exceed  a  maximum 
determined  by  magnetic  field  and  geometry.  With  a  1  meter 
radius  of  curvature  bend,  computer  orbit  calculations1*1 

indicate  that  a  stellerator  magnetic  field  of  1  kG/MeV  is 

v 

required  to  contain  single  particles. 

High  energy  operation  requires  a  time  dependent  vertical 
field  to  reduce  the  curvature  drift  in  the  bends.  'Hie 
stellerator  fields  then  serve  to  provide  a  stable  energy 
bandwidth  of  1  MeV/kG  about  the  equilibrium  orbit. 

For  FEL  applications,  a  wigglar  can  be  inserted  and  an 
optical  cavity  formed  in  one  leg  of  the  accelerator  provided 
the  gain  is  not  too  low  as  a  result  of  energy  spread  in  the 
beam.  One  can  achieve  variable  wavelength  operation  over  a 
wide  range  during  the  acceleration  phase  or  long  pulse 
narrow  band,  operation  after  acceleration. 

As  a  result  of  Inserting  the  long  pulse  induction 
accelerator  module  in  a  racetrack  geometry  with  stellerator 
fields  we  have  the  possibility  of  a  factor  of  40  gain  in  the 
energy.  This  could  result  in  a  high  gain  free  electron  •: 
laser  in  the  infrared  region. 
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THE  MODIFIED  BETATRON 

The  betatron  is  a  compact  cyclic  accelerator  with  a 

well-known  technological  background.  Betatrons  have  long 

been  used  to  accelerate  electrons  to  high  energies  300 

MeV).  However,  for  a  bean  of  energy  me  y, the  current  in 

these  devices  is  quite  low  because  of  the  space  charge  limit 

(I.  *  V*)  and  because  of  various  beam  instabilities  (which 

1  in 

typically  scale  as  Ij^  *  y)  .  Consequently,  a  conventional 
betatron  is  not  suitable  to  accelerate  large  current 
(>  several  kA)  electron  beams  that  are  injected  at  low 
energy  (y  <  10)  . 

The  modified  betatron  accelerator  overcomes  the  current 
limitation  inherent  in  the  conventional  betatron  by  the 
addition  of  a  strong  toroidal  magnetic  field  component 
parallel  to  the  particle  orbits  (Figure  4).  It  can  be  shown 
from  space  charge  considerations  alone  that  the  maximum 
number  of  electrons  that  can  be  confined  in  a  modified 
betatron  is  larger  than  the  corresponding  number  for  a 
conventional  betatron  by  a  factor  (Bg/3^)  , 

where  B.  and  B  are  the  toroidal  and  vertical  or  betatron 
magnetic  fields,  respectively.  The  maximum  current  that  can 
be  confined  in  a  modified  betatron  then  is 

r.  2  B  2 

Ima*-2-1(r)  T  to  [k*5' 

o  * 

where  rb  and  r^  are  the  minor  and  major  beam  radii.  Thus, 
for  rb  *  1  cm,  rQ  «  1  m,  y  *  5,  Bx  ■  100  G,  and  B^  «  5  kG, 
the  maximum  current  is  about  65  kA.  Although  this  value 
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represents  the  limit  due  only  to  space  charge  effects,  it 
can  be  shown  that  the  constraints  imposed  by  stability 
considerations  are  also  substantially  reduced  by  the 
imposition  of  a  toroidal  field. 

In  Figure  4,  a  possible  injection  and  ejection  scheme  is 
depicted.  Electrons  from  an  external  source  are  injected 
along  open  field  lines  into  the  torus.  Conditions  can  be 
chosen  so  that  the  forces  on  the  beam  arising  from  both  the 
external  and  self  fields  cause  the  beam  to  drift  inward  on 
an  open  orbit.  When  it  reaches  the  center  of  the  torus, 
an  external  parameter  is  changed  so  that  the  beam  is  on  a 
closed  orbit  and  thus  is  trapped.  One  relatively  simple  way 
to  achieve  this  trapping  is  to  change  the  external  field 
index  n,  which  determines  the  variation  of  with  r» 
i.e.  B^Ir)  «  r  n  .  After  trapping,  the  beam  is  accelerated 
just  as  in  a  conventional  betatron  by  increasing  the 
vertical  field.  After  acceleration,  the  beam  can  be 
extracted  by  suddenly  energizing  the  field  annihilation 
coils,  thereby  allowing  the  beam  to  leave  the  torus  along  a 
tangent  to  the  orbit. 

A  proof  of  principle  experiment  is  now  being  designed  at 
NRL  in  which  a  modified  betatron  will  be  used  to  accelerate 

a  10  kA  electron  beam  from  ~3  MeV  to  a  final  energy  of  ~  SO 

1 2 

MeV.  Theoretical  analyses  and  computer  simulations  have 
indicated  that  such  parameters  are  feasible  for  aim  radius 
beam.  The  initial  and  final  betatron  fields  would  then 
be  ~  170  G  and  -*1.7  kG,  respectively,  with  a  toroidal 
field  of  ~  10  kG. 
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The  applicability  of  the  modified  betatron  to  FEL 
experiments  will  depend  to  a  large  extent  on  the  quality  of 
the  extracted  beam.  Although  little  work  has  been  done  on 
this  aspect  of  the  problem,  the  current  and  energy  levels 
that  are  theoretically  achievable  could  make  possible  a 
high-gain,  infrared  FEL. 

SUMMARY 


We  have  briefly  discussed  three  approaches  to  high 
current  beam  accelerators  presently  being  investigated  at 
the  Naval  Research  Laboratory.  These  accelerators  are 
scalable  to  energies  in  excess  of  those  obtainable  with  a 
Marx-pulse  line  beam  generator.  The  induction  linac  at  NRL 
is  unique  because  of  its  long  pulse  length.  Although  the 
linac  is  scalable  to  high  energies,  the  cost  encourages  one 
to  examine  the  possibility  of  high  current  cyclic 
accelerators.  Two  high  current  cyclic  accelerator  concepts 
that  are  scalable  to  high  energies  have  been  discussed. 
These  devices  are  presently  in  the  conceptual  design  phase 
at  NRL.  The  suitability  of  these  accelerators  as  FEL 
drivers  must  await  their  development. 

A  number  of  FEL  experiments  have  been  planned  for  the 
long  pulse  induction  linac.  The  design  parameters  of  the 
first  such  experiment  is  discussed  in  this  paper.  This 
experiment  is  expected  to  be  in  operation  in  the  near 
future. 
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